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impact;  b)  the  effect  of  fabric  finishes;  and  c)  the  effect  of  weave.  A 
simple  model  of  impact  together  with  the  experimental  data  showed  that 
inter-layer  effects  are  of  second  order.  High  friction  finishes  and  bonding 
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agents  vara  found  to  be  detrimental  to  the  performance  of  Kevlar.  Triaxial 
fabrics  were  found  to  be  inferior  to  biaxial  fabrics.  Part  II  describes 
the  development  and  applications  of  computer  models  of  impact  using  both 
characteristics  and  finite  element  methods. 
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GENERAL  INTRODUCTION 


This  report i  presented  in  two  parts,  describes  the  investigation 
carried  out  in  UMIST  in  the  period  1979-1991  into  the  performance  of 
textile  structures  under  ballistic  impact.  The  investigation 
follows  two  complementary  routes,  namely  an  exhaustive  program  of 
experimentation  centred  around  ballistic  impact  tests  and  secondly  a 
detailed  survey  and  development  of  theoretical  and  computer  models  for 
simulation  of  the  textile  and  projectile  dynamics.  The  design  and 
evaluation  of  textile  structures  for  protection  against  high  speed 
normal  impact  must  hinge  upon  the  many  material,  surface  and  structure 
parameters  intrinsic  in  the  structure  design.  -  The  projectile, 
impacting  in  the  range  150  m/s  -  550  m/s,  and  with  a  nominal  mass  of  1 
gm,  encounters  the  textile  structure  both  at  the  fibre  level  and  at  the 
weave  level.  In  the  first  instance,  one  must  consider  firstly  for  the 
arrest  of  the  projectile,  resistance  to  cutting  by  the  projectile 
edges,  and  opening  of  the  weave  structures.  These  are  not  important 
considerations  in  the  present  configuration,  but  would  obviously  be 
considered  for  very  small  projectiles  where  the  projectile  diameter  is 
of  the  order  of  the  weave  pitch.  A  more  important  consideration  is  in 
the  indentation  jf  fabric  surface,  leading  to  high  fibre  strains,  and 
when  these  strains  exceed  the  fibre  fracture  strain,  the  structure  will 
fail.  The  structure  is  also  considered  inadequate  if  the  total 
indentation  is  sufficiently  large  as  to  present  a  hazard  to  the  sensitive 
zone  behind  the  fabric,  even  through  the  fibre  fracture  strain  is  never 
exceeded.  Important  to  these  two  criteria  are  firstly  fibre  material 
properties  quantified  by  fracture  strain,  modulus  of  elasticity  and 
fibre  weight  (denier).  Then  the  method  of  assembly  into  the  textile 


material  is  important,  quantified  by  the  warp  and  weft  pitch  leading  to 
the  area  density  (textile  weight)  and  crisp  parameters,  by  the  fibre 
friction  which  impede  the  relative  motions  of  the  fibres  in 
accommodating  the  impact  energy  and  finally  by  the  assembly  of  textile 
layers  into  a  structure  considering  layer  mixes  and  layer  fixing  by 
stitching  or  by  resin  bonding.  The  important  projectile  impact 
parameters  are  projectile  mass,  and  impact  velocity  or  energy. 

In  the  first  part  of  this  report,  the  results  of  an  experimental 
research  programme  are  presented.  The  main  aspects  of  this  work  were 
investigations  of  the  effects  of  the  addition  of  finishes  designed  to 
alter  the  degree  of  yarn  to  yarn  cohesion  within  a  fabric  and  also  the 
compare  the  multi-ply  response  to  impact  to  that  of  the  single  layer. 

In  the  second  part  of  this  report,  various  theoretical 
considerations  are  examined  with  the  motives  of  generating  procedures 
for  predicting  and  stimulating  textile  structure  performances.  The 
first  considerations  examine  the  mode  of  propagation  to  stress  signals 
through  the  structure,  leading  to  the  ideas  in  fibre  matching  at  the 
nodes  (intersection  of  warp  and  wefts)  and  then  describes  the 
variational  model  applied  to  triaxlal  fabrics.  This  model  has 
previously  been  developed  and  applied  to  regular  orthoganal  weaves 
(Journal  of  Textile  Institute  Vol  70,  No.  Ill,  page  469  et  seq,  1972) 
and  is  applicable  in  the  very  approximate  sense  to  single  and  then 
multilayer  structures  with  no  slip  either  between  fibres  or  layers. 
Then,  the  method  of  characteristics  is  examined,  this  leading  to  a 
computer  simulation  which  is  very  applicable  to  single  layer  structures 
In  which  the  fibre  materials  and  fibre  orientation  are  significant.  The 
high  speed  impact  results  are  quite  accurate  since  this  model  is  stress 
wave  based.  The  next  set  of  models  are  finite  element  models, 
focussing  on  various  details  in  the  structure.  The  first  modal,  the 


membrane  model  Incorporates  the  crimp  effects  into  the  element  and  is 
appropriate  to  single  layer  structures.  The  second  looks  at  the  detail 
of  the  fibres  in  the  node  area.  The  third  and  fourth  programs  consider 
the  multilayer  structure  the  former  considering  the  layers  bonded  by  a 
resin  or  glue  effective  only  in  shear,  and  the  latter  the  more  important 
of  the  two,  allowing  layers  to  slide  relative  to  each  other. 
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1.1  Introduction 


The  aims  of  the  present  study  were  to  construct  a  ballistic  range 
capable  of  impacting  1  gm  steel  projectiles  onto  fabric  samples  and  to 
also  provide  apparatus  for  the  measurement  of  impact  velocity  and  the 
energy  absorbed  by  penetration.  Apparatus  was  also  developed  in  order 
to  (1)  provide  multi-flash  photographs  of  fabrics  under  impact  (2) 
measure  the  arrival  of  the  stress  wave  at  the  fabric  boundary.  In 
addition  work  was  undertaken  to: 

1 )  establish  the  effects  on  the  ballistic  performance  of  the  addition 

\ 

to  fabrics  of  chemical  finishes  designed  to  alter  the  degree  of 
yarn  to  yarn  cohesion  within  a  fabric; 

2)  compare  the  multi-ply  response  to  impact  to  that  of  the  single 
layer  over  a  range  of  impact  velocities. 

These  two  factors  were  known  to  be  of  importance  because: 

1)  Laible1  has  found  that  the  addition  of  high  friction  finishes 

significantly  improved  the  ballistic  performance  of  polypropylene 

2 

fabrics;  alternatively  Morrison  has  shown  that  when  Kevlar  is 
used  in  either  a  polyester  resin  or  a  silicone  rubber  composite, 
the  ballistic  resistance  is  reduced  when  compared  to  the  untreated 
fabric ; 

2)  most  theoretical  work  has  only  been  developed  to  the  level  of 
single  layer  impact. 

In  addition  to  the  ballistic  performances  of  the  above,  two 
trlaxlal  weave  fabrics  and  a  lightweight  nylon  fabric  were  compared  to  a 
heavyweight  biaxial  fabric  more  representative  of  those  normally  used 
for  ballistic  protection. 


1 .2  Development  of  the  Ballistic  Range 


During  the  previous  research  programme  supported  by  the  Ministry 
of  Defence,  SCRDE,  facilities  were  developed  to  allow  transverse 
ballistic  impact  onto  fabrics.  Lead  airgun  pellets  were  fired  from  a 
cartridge  gun  at  velocities  within  the  range  100  m  s_1  to  550  m  s-1. 

The  velocities  of  the  projectile,  before  impact  and  after 

penetration  of  the  fabric  (VQut) ,  were  measured  in  the  following 

manner.  Four  pencil  leads  were  aligned  along  the  path  of  the 

projectile;  two  in  front  of  the  fabric,  and  two  behind.  The  pencil 

leads  formed  part  of  a  simple  resistive  circuit  uith  an  input  to  a 

storage  oscilliscope.  As  each  pencil  lead  in  turn  was  broken  by  the 

projectile,  the  voltage  input  to  the  oscilloscope  was  consequently 

reduced.  From  the  stored  dljplay  the  time  between  contact  breaks  am 

thus  to  T.  and  V  t  velocities,  were  obtained, 
in  out 

Photographs  of  the  deformation  process  after  impact  but  before 
penetration  were  obtained  by  illuminating  the  event  with  a  single  flash 
of  light  of  1.2  u  s  duration.  These  experiments  were  performed  in  a 
darkened  room  with  the  camera  lens  open.  The  time  of  the  flash,  in 
relation  to  a  contact  break  a  short  distance  before  impact,  was  varied 
using  a  simple  delay  circuit.  By  varying  this  delay  for  each  trial, 
plots  of  pellet  position  against  time  were  obtained.  This  method, 
however,  gave  a  large  scatter  of  results  as  both  a  constant  impact 
velocity  for  each  trial,  and  identical  modes  of  penetration  for  each 
sample,  from  the  same  fabric,  had  to  be  assumed. 

This  apparatus  was  modified  for  the  present  research  in  order  to 
minimise  the  energy  absorption  due  to  projectile  deformation.  Steel, 
rather  than  lead  projectiles  were  required.  Cylindrical  steel 
projectiles  of  diameter  5.5  mm,  length  5.5  mm  ar.d  mass  1.004  (+  0.0008)g. 


were  obtained.  The  projectiles  are  located  in  plastic  sabots  and 
fired  from  a  0.303  inch  (.769  nan)  rifle  barrel.  So  far,  muzzle 
velocities*  of  between  262  ms”^  and  550  ms”^  have  been  achieved. 
Multiple  flash  photographs  show  that  before  impact  the  projectiles  yaw 
slightly  but  do  not  spin.  As  there  was  a  danger  that  the  steel 
projectiles  might  richochet,  the  whole  apparatus  was  enclosed  in  a 
steel  box.  Polycarbonate  windows  were  included  in  the  structure  to 
allow  the  fabric  sample  to  be  photographed.  The  present  apparatus  is 
shown  in  Figure  1. 

In  order  to  trigger  a  aeries  of  light  flashes  to  obtain  multiple 
image  photographs  of  the  projectile-fabric  interaction;  a  'Bowen'  ten 
channel  delay  was  acquired. 

This  unit  produces  electronic  delays,  in  either  parallel  or 
series,  down  to  intervals  of1us(±0»01f)»  A  projectile  moving  at  a 
velocity  of  say  550  m  s”1  moves  .55  mm  in  1  us,  thus  six  stages  of  fabric 
deformation  can,  in  theory  be  observed  before  the  pyramidal  deformation 
has  reached  3-3  nm.  Therefore  the  delay  is  more  than  adequate  for  the 
present  research. 

Previously  a  single  flash  of  a  stroboscope  was  used  as  the  light 
source  for  photography.  Although  this  flash  was  of  a  relatively  short 
duration  and  of  reasonable  brightness,  image  quality  was  impaired  due 
to  either;  (a)  blur  caused  by  projectile  movement  during  the  time  of 
the  flash,  b)  a  small  depth  of  focus  due  to  the  large  lens  apertures 
required  for  the  low  light  level  or  (c)  high  graininess  due  to  the  fast 
film  which  was  required.  These  effects  were  exaggerated  when  negatives 
were  greatly  enlarged  for  analysis.  In  order  to  improve  image  quality, 
three  'Pulse'  double  flash  argon  arc  light  sources  have  been  obtained. 
The  manufacturer  quotes  a  flash  duration  of  0.3  tis. 


Two  still  cameras,  for  exclusive  use  on  the  present  research 
programme,  were  purchased  to  enable  simultaneous  'side  on'  and  'end  on' 
photographs  of  the  deformation  to  be  obtained. 

Previously  there  was  a  relatively  large  error  in  the  measured 
values  of  Vin  and  It  is  particularly  important  to  reduce  the 

errors  of  these  measurements,  as  the  energy  loss  of  the  projectile  is 
obtained  from  the  difference  of  the  squares  of  V^n  and  V 
Previously,  time  intervals  could  only  be  measured  to  an  accuracy  of  ± 
2%.  In  order  to  improve  on  this  two  digital  timers  were  purchased. 
These  can  each  measure  time  intervals  down  to  0.1  u  s.  The  separation 
of  the  velocity  measuring  stations  are  now  of  the  order  of  80  cm. 
Therefore,  timing  errors  now  contribute  to  only  0.07%  of  the  error  in 
velocity  measurement  for  a  velocity  of  550  ms”1. 

The  above  facilities  were  used  to  obtain  micrographs  of  fabric 
systems  under  Impact. 

Figure  2  shows  three  such  single  flash  photographs  of  the 
deformation  of  a  double- layered  nylon  fabric.  The  impact  velocity  is 
approximately  350  ms”1  in  each  case.  A  general  view  of  the  whole 
fabric  is  shown  in  figure  2A.  Lines  were  drawn  on  this  fabric  at  0.5  cm 
separation,  parallel  to  both  the  warp  and  weft  directions.  From  this 
photograph  the  pyramidal  form  of  the  transverse  deformation  is  clearly 
seen.  A  line  drawn  on  the  print  parallel  to  one  of  the  grid  lines  helps 
to  show  the  in-plane  displacement  of  the  fabric  towards  the  point  of 
impact.  Figures  2B  and  2C  show  the  region  of  the  fabric  near  the 
impact  zone  at  approximately  20  us  and  60  us  after  impact.  In  figure 
2B  the  projectile,  which  is  moving  in  a  direction  from  bottom  to  top  of 
the  photograph,  has  not  yet  penetrated  the  fabric  while  in  2C  the  whole 
of  the  cylindrical  projectile  is  visible. 


f  **: 


(figures  5B  and  5C)  as  much  original  detail  is  lost  during  reproduction. 
In  figure  5A  the  time  between  flashes  is  10  ys  while  in  figure  SB  this 
interval  is  20  us.  Figure  5A  shows  the  fabric  at  times  before 
penetration.  In  figure  5B,  however,  one  image  shows  the  fabric  before 
and  two  images  show  the  fabric  after  penetration.  The  multiple  flash 
technique  has  been  used  by  previous  researchers  in  the  same  field. 
However  much  previous  work  using  multiple  flash  methods  has  been  at  low 
impact  velocities  in  order  to  be  able  to  observe  large  deformations 
before  penetration.  In  figure  5A,  the  pyramid  height  is  only  0.8  cm  at 
the  time  of  failure.  (Approximately  50  us  after  impact). 

Strain  Measurement 

It  would  be  of  enormous  theoretical  interest  to  be  able  to  measure 
the  spacial  and  temporal  distribution  of  statin  during  impact.  The 
measurement  of  the  level  and  distribution  of  stress  and  strain 
throughout  a  fabric  during  impact,  does  however,  present  serious 
experimental  problems,  as  the  measuring  instruments  could  easily 
themselves  influence  the  strain  distribution.  A  partial  solution  is  to 
restrict  measurements  to  the  fabric  boundary.  For  this  purpose  a 
commercial  pressure  transducer  has  been  mounted  in  such  a  manner  as  to 
enable  tensile  forces  to  be  measured.  One  yarn  from  the  fabric  was 
clamped  to  the  transducer.  The  sample  holder  was  then  adjusted  in  order 
to  ensure  that  the  yams  intersected  the  impact  zone. 

The  force  transducer  alone  has  a  resonant  frequency  of  200  kHz. 
The  output  is  therefore  filtered  by  a  180  kHz  low  pass  filter.  Stress 
wave  reflections  at  the  boundary  would  be  expected  to  give  rise  to 
frequencies  (f)  in  the  region  of 


where  v  a  sonic  velocity 

d  a  distance  from  impact  zone  to  fabric  edge. 


orlop 


This  ia  only  an  approximation,  aa  no  consideration  was  made  of 


For  Kevlar 

f 


4  x  103 
0.075x2 


26.6  kHz 


Thus  the  apparatus  was  adequate  for  the  present  research. 

The  output  from  the  transducer  was  fed  into  a  transient  recorded 
via  a  change  amplifier.  The  resulting  data  was  displayed  and  then 
stored  on  floppy  disc. 

To  allow  the  force  measurements  to  be  synchronised  with 
micrographs  of  the  fabric  during  impact,  the  output  from  a  photodiode 
was  displayed  on  one  channel  of  the  transient  recorder.  A  typical 
trace  is  shown  in  figure  6.  Ore  yarn  from  the  last  layer  of  a 
multilayered  sample  was  attached  to  the  force  transducer.  The 
uppermost  trace  is  the  output  from  the  force  transducer.  Below,  the 
photodiode  output  is  shown.  The  large  decay  time  of  the  photodiode 
output  does  not  reflect  the  light  flash  duration,  but  is  merely 
attributable  to  the  electrical  system.  The  peak  output  occurs  as  the 
projectile  is  arrested. 

Figure  7  shows  the  synchronised  output  from  another  trial.  Here, 
although  the  cone  deformation  has  already  occurred  and  the  photograph 
shows  the  fabric  at  a  time  of  at  least  20  us  after  impact,  no  signal  is 
recorded  by  the  transducer  until  22  u  s  later.  It  can  be  seen  from  the 
micrograph  that  the  fabric  is  not  completely  flat  in  the  region  near  the 
transducer.  It  seems  that  this  has  resulted  in  a  delay  of  the  stress 
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Only  preliminary  work  has  so  far  been  undertaken  using  the 
multiple  flash  and  the  boundary  transducer  but  results  will  be  obtained 
on  the  deceleration  of  the  projectile  and  the  resulting  force  at  the 


boundary  in  the  following  research  programme. 

hi  The  Ballistic  Performance  of  a  Lightweight  Fabric 

During  the  course  of  the  current  project  the  author  was  asked  to 

compare  a  lightweight  and  a  heavyweight  nylon  fabric  on  the  ballistic 

range  as  a  part  of  preliminary  work  for  a  proposed  project.  Details  of 

4 

the  two  fabrics  are  shown  in  Table  1.  Figucia  reports  that  the 
ballistic  performance  of  a  multilayered  fabric  of  given  weight 
increases  "as  the  area  density  of  the  individual  layers  is  reduced. 

Fabric  Area  Tarns  Tex  of  Mean  Mean 

Density  per  cm.  yarns  breaking  load  }  extension 

em-2  <**« 

Warp  Weft  Warp  Weft  Warp  Weft  Warp  Weft 

42.8  38.7  39.7  5.40  5.19  0.137  0.132  17.0  16.1 

378.0  9.75  9.75  197.1  193.6  13.3  13-6  23.8  23.6 

Table  1 


The  fabrics  were  tested  on  the  ballistic  range.  The  results  are 
presented  as  a  plot  of  the  energy  loss  against  the  combined  area  density 
of  the  multilayered  samples,  (fig.  8).  It  is  clear  that  weight  for 
weight,  less  energy  is  absorbed  by  the  lightweight  fabric.  It  is  of 
interest  to  compare  the  V^0  energies.  The  data  from  tests  on  28  layers 
(s  1198  gm  )  of  the  lightweight  fabric  are  shown  as  a  plot  of  absorbed 
energy  against  impact  velocity  (fig.  8b).  The  V5Q  energy  was  found  to 
From  the  Ve0  energy  curve  of  the  heavyweight  fabric,  the 


be  27.3  J 


Combined  Area  Density  of  Multilayer  Sample- kg  m* 

FIG.  8  Energy 
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V^q  energy  at  1198  gm  was  calculated  to  be  40.3  J.  Thi3  difference 

could  be  a  reflection  of  the  differences  in  the  quasi-static  work  of 

rupture  of  the  constituent  yarns.  Assuming  that  the  tensile  3train 

energy  is  the  only  mode  of  energy  absorption  and  that  only  yarns  which 

make  contact  with  the  projectile  are  strained  then  the  expected 

projectile  energy  loss  per  unit  area  density  was  calculated  to  be  - 

2 

heavy  duty  nylon  =  5.29  J  per  Kg/m 

2 

lightweight  nylon  =  3*39  J  per  Kg/m 
(Typical  quasi-static  stress  strain  curves  are  shown  in  figure  9). 
These  results  are  only  one  sixth  of  the  V^Q  energies  found  by 
experiment,  however,  the  ratios  of  the  experimental  md  calculated 
energies  are  similar: 

5  29 

Ratio  of  specific  work  of  rupture  of  principal  yarns  -  ^  •  *  =  1.56 

40 . 3 

Ratio  of  VgQ  energies  -  273  3  1 • 48 

However,  the  ballistic  performance  of  the  two  nylon  fabrics  might 
have  been  expected  to  be  more  similar  as  the  work  of  rupture  is  not  the 
only  mode  of  energy  absorption.  For  example,  the  transverse  kinetic 
energy  of  the  fabric  accounts  for  some  of  the  energy  absorbed  from  the 
projectile. 

One  possible  cause  of  the  inferior  ballistic  performance  of  the 
lightweight  fabric  is  a  difference  in  the  mode  of  penetration.  Many 
photgraphs  of  the  projectile  in  flight  after  penetration  revealed  a 
small  object  travelling  ahead  of  the  projectile.  This  was  identified 
as  a  fused  plug  of  nylon.  It  is  known^’^  that  fibre  melting  is  a 
factor  in  the  failure  mechanism  but  no  plugging  has  been  observed 
before.  The  nylon  plug  travels  ahead  of  the  impact  projectile  but  was 
slowed  due  to  air  friction.  It  must  be  concluded  that  the  plug  was 
ejected  before  the  final  disengagement  of  the  projectile  from  the 
fabric.  The  exact  mechanism  is  however  unknown. 


Typical  quais-static  stress  strain  curve; 
for  yams  extracted  from  (A)  a  heavyweipi 
and  (B)  a  lightweight  fabric. 


1.4  The  Ballistic  Performance  of  Two  Triaxlal  Fabrics 


In  normal  (biaxial)  fabrics  two  seta  of  parallel  threads  intersect 
at  90°.  In  triaxlal  fabrics,  three  sets  of  parallel  threads  intersect 
at  60°  intervals.  There  are  reasons  for  believing  that  the  ballistic 
performances  of  otherwise  identical  fabrics  would  be  different  for 
biaxial  and  triaxlal  fabrics,  although  as  yet  the  superiority  of  either 
cannot  be  confidently  demonstrated  due  to  lack  of  knowledge  of  the 
importance  of  various  effects  of  weave.  For  example  a  triaxlal  weave 
might  confer  added  ballistic  performance  due  to  the  spreading  of  load 
more  evenly  around  the  impact  zone  or  alternatively  might  confer  a 
reduced  ballistic  performance  due  to  the  inherent  openness  of  the 
triaxlal  weave. 

Two  triaxlal  fabrics  (A  and  B)  were  therefore  compared  on  the 
ballistic  range.  First  consider  fabric  A  which  became  available  early 
on  in  the  research  program.  As  the  ballistic  rig  was  only  partially 
operational  at  the  time,  and  the  fabric  quantity  was  limited  the 
triaxlal  fabric  was  compared  to  a  biaxial  fabric  using  a  simple  ad  hoc 
procedure.  The  fabric  parameters  are  shown  in  Table  2.  The  ballistic 
performance  of  each  fabric  was  assessed  by  finding  the  number  of  layers 
required  to  stop  the  projectile  at  a  known  velocity. 

The  results  are  shown  in  Table  3*  It  is  clear  that  the  total  area 

•2 

density  of  the  biaxial  fabric  (2400  gm  ),  which  is  required  to  stop  the 
projectile,  is  much  less  than  that  required  in  the  case  of  the  triaxlal 
fabric  (3094  ms"1). 


’  y 


Area 

density 

,_-2x 


Tarns  per  ca  Tex  of  yarns 


Breaking  load 
of  yarn 
(Newtons) 


Breaking  Extension 

(*) 


WARP 

WEFT 

WARP 

WEFT 

WARP 

WEFT 

WARP  WEFT 

-iaxial 

§ylon 

240 

11.4 

12.3 

103 

101 

67.1 

68.5 

24.9  25.5 

A 

B 

C 

A 

B  C 

A  B 

C 

ABC 

riaxial 

jprloa 

221 

3.8 

3.9 

3.7 

200 

205  201 

109  128 

130 

20.8  22.2  23 .2 

Table  2 


Although  the  fabrics  had  similar  area  densities,  it  cannot  be 
concluded  that  the  difference  in  construction  was  the  sole  reason  for 
the  difference  in  ballistic  resistance,  as  the  triaxlal  was  of  a  very 
open  structure.  A  triaxlal  fabrio  of  a  much  closer  weave  was  therefore 
tested  balllstlcally. 


Fabric  No.  of  layers  Area  Density  Impact  velocity  Was  sample 

of  sample  ,  -K  penetrated  *'  i 

VDL9  ) 


(ga-2) 


Biaxial 

nylon 


■■u  ... 

t  T 


Triaxlal 

nylon 


Table 
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The  area  dena.J  /  of  three  layered 
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sample  was  1281  go  .  The  VgQ  energy  was  found  to  be  38.3  J.  The 
equivalent  VgQ  limit  of  the  heavy  duty  nylon  described  earlier  can  be 
found  from  the  plot  of  VgQ  energy  against  area  density  (figure  8).  By 
this  method  the  V^q  limit  for  the  biaxial  nylon  of  combined  area  density 
1281  gm"2  is  42.5  J. 

The  energy  of  the  triaxial  fabric  was  thus  marginally  lower 
than  the  biaxial  fabric ,  but  this  is  not  statistically  significant.  It 
is  concluded  that  the  performance  of  the  biaxial  and  triaxial  fabrics 
were  similar. 

In  part  II  of  this  report  it  is  demonstrated  that  a  computer  model 
of  the  impact  onto  triaxial  structures  predicts  that  the  triaxial  weave 
is  ballistlcally  inferior  to  the  biaxial  weave.  Conversely  it  might  be 
expeoted  that  the  triaxial  fabrio  might  be  more  efficient  by 
considering  differences  in  the  transverse  wave  front. 

Adeyefa  predicts  that  for  triaxial  weaves,  the  transverse  wave 
front  would  be  hexagonal.  A  photograph  of  the  reverse  side  of  the 
triaxial  fabric  during  impact  is  shown  in  figure  10.  The  wave  front  is 
confirmed  to  be  hexagonal.  Figure  10b shows  a  biaxial  nylon  and  a 
biaxial  Kevlar  by  way  of  comparison.  In  these  cases  the  transverse 
wave  front  is  rhomboidal. 

By  geometrical  considerations  it  is  found  that  the  area  bounded  by 
the  transverse  wave  front  is  greater  for  the  triaxial  wave  over  the 
biaxial  weave  by  1  :  1.3.  It  is  expected  that  the  kinetic  energy  of  the 
transversely  moving  region  of  the  fabric  would  be  greater  for  triaxial 
than  biaxial  weaves.  This  effect  would  make  the  triaxial  weave  a 
better  absorbant  of  projectile  energy. 


FIG.  10b 


A  Biaxial  nylon  (top)  and  a  Biaxial  Kevlar  (bottom)  jM 

under  ballistic  impact. 
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1.5  Main  Research  Programme 


Experimental  Procedure 
For  the  purposes  of  the  present  research  two  methods  of 
quantifying  the  ballistic  performance  of  textile  materials  have  been 
adopted. 

First,  the  variation  of  energy  loss  against  Impact  velocity  Is 


established,  and  from  this  the  V^q  limit  Is  obtained.  The  variation  of 


Impact  velocity  provides  an  Independent  Input  parameter  for  models. 
The  experimental  and  theoretical  results  can  then  be  compared.  This 
method  also  provides  a  means  of  comparing  the  ballistic  properties  of 
different  fabrics.  The  V^q  limit  can  be  found  for  samples  composed  of 


! 


different  number  of  layers.  The  variation  of  the  limit  with 


the 


combined  area  density  of  the  multilayered  sample  enables  the  ballistic 

i 

performance  to  be  quantified.  This  method  has  the  advantage  that  the 

l 

Vg0  limit  is  a  widely  accepted  way  of  assessing  ballistic  resistance.  A 

.  ■ 

great  disadvantage  Is  the  large  number  of  tests  required  to  establish 

i 

the  V^q  limit  even  for  a  single  layered  sample.  Difficulties  also 

arise  due  to  the  large  scatter  of  results  at  velocities  near  the  iV__ 

i  50 


limit. 


The  second  technique  is  to  Impact  a  series  of  multilayered  samples 
at  one  impact  velocity  (approx.  500  m  s"1).  In  this  procedure  fabric 
performance  is  established  by  the  variation  of  energy  absorbed  with  the 
combined  area  density  of  the  multilayer  sample.  Figucia  (ref.  4) 
reports  that,  for  Kevlar,  this  relation  is  linear.  Figucia  also 
reports  that  the  projectile  energy  loss  was  found  to  be  independent  of 
impact  velocity  in  the  range  tested.  Figucia  only  considers  results 
where  less  than  50-6051  of  the  available  energy  is  absorbed.  The 


author’s  technique  was  to  include  all  results  for  velocities  where 
penetration  occurred,  for  it  has  been  found  that  although  the  absolute 


value  of  the  experimental  scatter  was  high,  the  fractional  variation  in 
absorbed  energy  was  acceptable.  The  second  method  has  proved  very 
successful  in  quantifying  fabric  differences. 

Fabric  Finishes 

Many  computer  models  have  been  developed  to  predict  behaviour  of 
woven  textiles  when  subjected  to  impact  by  a  free  flying  projectile.  A 
common  approach  is  to  treat  the  woven  fabric  as  a  pin  jointed  net. 
However ,  previous  research  by  the  author  has  produced  evidence  that 
longitudinal  movements  of  the  principal  yarns  can  occur  relative  to  the 
crossing  yarns.  The  evidence  originates  from  observations  of  residual 
displacements,  and  residual  damage  of  the  principal  yarns  (ref.  (?).  It 
is  expected  that  changing  the  ability  of  yarns  to  pass  freely  over  each- 
other  might  change  the  energy  absorbing  characteristics  of  the  fabric. 
Three  possible  effects  are  postulated.  First,  an  increase  in  crossover 
stability  will  cause  an  increase  in  the  reflected  part  of  the 
longitudinal  wave  at  each  crossover.  This  leads  to  increased  :  train  at 
the  impact  zone,  and  early  failure.  Second,  an  increase  in  crossover 
stability  increases  the  energy  transfer  to  parallel  yarns.  Third,  a 
decrease  in  crossover  stability  would  increase  the  ability  of  yarns  to 
move  laterally  at  the  impact  zone  and  allow  the  projectile  to  pass 

without  yarn  failure. 

2 

Morrison  has  shown  that  the  addition  of  a  rubber  matrix  to  a 
Kevlar  fabric  seriously  reduces  the  ballistic  performance.  Morrison 
concludes  that  for  the  materials  tested,  the  second  effect  i.e.  the 
Increase  of  nodal  reflections  must  be  predominant. 

Alternatively,  it  is  known  that  wet  samples  of  Kevlar  are 
seriously  weaker  than  dry  samples  (ref“?).  One  possible  explanation  is 
that  water  has  seriously  reduced  the  yarn  to  yarn  friction.  Alterations 
to  the  stress-strain  behaviour  of  wetted  Kevlar  were  not  found  to  be 
large  enough  to  explain  this. 


Laible1  has  shown  experimentally  that  the  ballistic  performance  of 
polypropylene  fabrics  is  improved  by  the  addition  of  high  friction 
finishes. 

Lateral  yarn  movements  have  been  observed  by  the  author  (ref  £). 
In  some  cases  single  layer  fabrics  were  penetrated  without  any  yarns 
being  broken. 

The  above  observations  indicate  that  there  is  an  optimum  yarn  to 
yarn  cohesion. 

In  order  to  investigate  whether  yarn  to  yarn  friction  affected  the 
ballistic  properties  of  fabrics  and  in  order  to  observe  whether  any  such 
effects  are  dependent  on  impact  velocity,  multilayer  fabrics  of  nylon 
and  Kevlar  were  tested  on  the  ballistic  range.  The  fabrics  had 
previously  been  treated  with  chemical  treatments  designed  to  either 
decrease  or  increase  the  degree  of  yarn  to  yarn  friction  at  crossovers 
within  the  fabric.  Ideally  it  would  have  been  desirable  to  find  a 
treatment  which  would  have  reduced  the  friction  in  Kevlar  fabrics,  but 
none  was  found.  Even  a  finish  (Slligen  E)  which  is  designed  to  reduce 
friction  was  found  to  Increase  the  level  of  friction  for  Kevlar.  The 
results  are  shown  in  figures  11  to  16  as  plots  of  the  energy  loss  of  the 
projectile  against  the  impact  velocity.  The  scatter  of  data  is  a 
consequence  of  real  differences  in  absorbed  energy  for  each  trial.  As 
no  measurements  -of  rebound  velocity  were  measured,  the  energy  loss  is 
strictly  the  loss  of  kinetic  energy  of  forward  momentum.  Thus,  for 
tests  where  the  fabric  was  not  penetrated,  the  data  points  lie  on  the 
curve  of  formula 

E  =  irnV2 

where  m  is  the  mass  of  the  projectile.  For  cases  where  the  fabric  was 
penetrated  the  curve  was  fitted  statistically  to  fit  the  exponential 
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The  choice  of  this  form  of  curve  is  somewhat  arbitrary,  however  the  aim 
was  to  be  able  to  apply  a  systematic  routine  to  evaluate  differences 
between  the  two  fabric  finishes.  When  the  constants  a  and  b  have  been 
evaluated,  the  V^q  limit  is  predicted  by  solving 

aebv  -  ImV2  *  0. 

The  energy  can  then  be  calculated.  The  relative  performance 
of  these  fabrics  can  be  found  from  the  plot  of  the  VgQ  energy  against 
the  combined  area  density  of  the  multilayered  fabric.  The  ballistic 
performance  for  these  tests  were  qualified  by  then  calculating  the  7^Q 
energy  loss  at  an  area  density  of  1000  gm”  .  This  value  was  chosen 
arbitrarily,  but  is  well  in  the  range  of  the  area  densities  tested  in 
the  current  research.  The  results  are  shown  in  table  4.  The 
following  conclusions  can  be  drawn: 

1)  the  high  friction  finish  conferred  an  improved  ballistic 
performance  on  the  nylon  fabric  compared  to  the  addition  of  a  low 
friction  finish; 

2)  both  the  high  friction  treatments  produced  a  large  reduction 
in  ballistic  performance  of  the  Kevlar  fabric; 

3)  the  water  treatment  improved  the  ballistic  performance  over 
the  'as  received'  fabric  possibly  because  residual  finishes  were 
removed ; 

4)  none  of  the  treatments  conferred  any  strong  velocity 
dependent  effects. 

In  view  of  this  last  conclusion,  it  was  decided  to  adopt  Figuoia's 
method  for  any  further  quantification  of  ballistic  performance.  In  the 
next  phase  of  research  a  more  open  weave  fabric  was  tested  as  it  was 
felt  that  the  inherent  freedom  of  yarns  would  provide  plenty  of  scope 
for  modification  of  yarn  to  yarn  cohesion.  Samples  of  these  fabrics 


were  immersed  in  aqueous  solution  of  the  fabric  finish  and  then  cured,  at 
120°  for  15  minutes.  As  a  control,  a  set  of  samples  was  immersed  in 
water  and  dried  at  the  same  temperature.  One  of  two  finishes  were 
appliedi- 

i)  LURAPRET  +  E30.  This  confers  a  high  yarn  to  yarn  friction. 
Chemically  it  is  an  aqueous  solution  of  silic  acid. 

(ii)  PARAPRET  HVN  +  LURAPRET  B30.  50/50  mixture.  Parapret  is 

an  aqueous  dispersion  of  an  anionic  polyacrylate. 

The  effects  of  these  treatments  were  quantified  using  the  yarn 
pull  test  described  in  the  previous  report.  A  typical  force-extension 
curve  obtained  from  this  test  is  shown  in  figure  17.  For  each  test, 
the  first  two  maxima  (A,B)  were  recorded.  The  results  for  the  Kevlar 
fabrics  treated  with  finishes  are  shown  in  table  5.  As  expected  the 
oaxlmum  forces  were  greatest  for  the  Lurapret-Parapret  (bonded)  sampler 
and  least  for  the  control  samples  (water  treatment).  For  the  '.-jaded 
samples  there  was  a  66.8}  reduction  in  force  between  the  first  and 
second  peaks,  whereas  for  the  other  samples  this  figure  was  lr.  the 
region  of  15}.  It  is  concluded  that  the  yarn  tc  yarn  bonds  are  broken 
at  the  first  peak  and  that  the  level  of  the  second  peaic  is  determined  by 
frictional  forces.  There  appears  to  be  very  little  yarn  to  yarn  bonding 
for  the  fabrics  which  were  treated  with  the  higher  friction  finish. 

I 

Reslilts  of  the  ballistic  tests  on  treated  fabrics 

Weight  for  weight  the  water  treated  samples  we  found  to  be 

significantly  superior  in  ballistic  performance  than  either  the  ’high 

friction'  or  bonded  samples.  The  latter  two  fabrics  were  found  to  be  of 

equal  ballistic  performance,  even  though  there  was  a  large  difference 

between  the  levels  of  force  in  the  yarn  pull  test.  The  BPI’s (Ball i sti 
performance  Indicator),  the  slope  of  the  energy  -  area  density  relation  is 

^Registered  trade  made  of  BASF,  England. 


listed  In  table  5.  The  results  are  presented  graphically  In  figure  18. 
These  results  hold  even  If  the  add  on  weight  of  the  finishes  is  not 
included  when  calculating  the  B.P. I,  although  the  differences  are 
lower.  The  reason  for  this,  of  course,  is  that  the  fabric  treatment 
adds  weight  but  not  strength.  These  renormalised  results  are  shown 
graphically  in  figure  19. 

Conclusions  on  the  effect  of  increased  yarn  to  yarn  friction 

Although  it  is  possible  that  similarity  in  ballistic  performance 
between  the  bonded  and  the  high  friction  finish  is  coincidental,  and  the 

some  intermediary  level  of  cohesion  would  produce  different  results;  it 

f 

|  seems  probable  that  there  exists  a  relatively  low  level  of  yarn  to  yarn 
cohesion  which  significantly  reduces  the  resistance  of  Kevlar  fabric. 
Above  this  level  the  ballistic  resistance  becomes  insensitive  to  any 
further  increase  in  friction  or  bonding,  alternatively  the  add-on 

'weight  of  the  finishes  could  have  modified  the  response  to  impact. 

| 

;  Fig. 8. 29a  and. 5.29b  show  strain  histories,  at  the  point  of  impact 
\  of  an  identical  fabric,  as  predicted  by  the  methods  described  in  part 
;  IX.  Fig?.  5.29a  shows  the  case  where  slipping  is  allowed  at  yarn 
crossovers  and  fig. 5. 29b  models  the  case  where  the  yarns  are  rigidly 
connected.  Although  the  time  to  projectile  arrest  is  less  for  the 
latter  case,  the  strains  are  generally  higher.  Thus  a  strain  failure 
criterion  is  applied  to  both  cases,  it  is  predicted  that  the  fabric  with 
rigidly  connected  yarns  would  fail  at  a  lower  impact  velocity  than  the 
fabric  where  yarn  slip  is  allowed.  This  is  consistent  with  the 
experimental  results. 

Ballistic  Resistance  of  Kevlar  29  and  Kevlar  49 

Kevlar  29  exhibits  twice  the  tensile  breaking  strain  but  almost 
half  the  tensile  modulus  of  Kevlar  49.  Roylance(ref .  3)  has  developed 
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a  finite  difference  model  for  the  ballistic  impact  onto  textile 
structures.  Using  these  quasi -static  fibre  parameters,  Roy lance 
predicts  that  Kevlar  29  is  superior  to  Kevlar  49.  However,  Roylance 
refers  to  experimental  results  which  seem  to  indicate  a  similarity  in 
ballistic  resistance . 

In  order  to  provide  our  own  data  a  Kevlar  29  and  a  Kevlar  49  fabric 
of  similar  construction  were  compared  on  the  ballistic  test  range.  Both 
fabrics  were  treated  with  the  high  friction  finish  as  the  results  were 
also  used  to  provide  data  on  the  effect  of  finishers. 

Results  of  Tests  on  Kevltr  49  and  Kevlar  29 

No  significant  difference  was  found  between  the  ballistic 
resistance  of  Kevlar  29  and  Kevlar  49.  Figures  20  and  21  are  plots  of 
projectile  energy  loss  against  combined  area  density  of  the  multilayer 
fabric  for  these  two  fabrics.  In  Figure  21  the  results  have  been 

_p 

renormalised  to  reflect  the  untreated  weight  of  the  fabrics  (175  gm  ). 

The  B.P.I's  are  shown  in  Table  5  . 

It  is  of  interest  to  compare  the  author's  results  with  those  of 

2 

Flgucia  (ref.  4)  who  found  B.P.I's  in  the  range  of  22.6  to  J1.5  J/kg/m  . 
The  two  untreated  Kevlar  samples  tested  by  the  author  were  found  to  have 
B.P.I's  of  29.6  and  33.6  J/kg/m  .  Thus  the  results  for  Kevlar  are  in 
good  agreement  even  though  the  authors  results  were  obtained  using  a 
blunt  cylindrical  projectile  instead  of  the  fragment  simulator  and  a  6" 
rather  an  a  1"  diameter  fabric  sample  size. 

Model  of  Multilayer  Impact 

A  simple  computer  model  was  used  to  calculate  and  display 
graphically  the  multilayer  performance  from  single  layer  data.  The 
model  assumed  no  physical  interaction  between  layers.  The  exit 
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AREA  DENSITY  (GRAMS  PER  SQUARE  METRE)  -  Based  on  untreated  weight. 

Projectile  energy  loss  (E)  against  the  combined  area  density  (A)  of  multilayered 
samples.  Comparison  between  the  ballistic  performance  of  Kevlar  20  and  Kevlar  49 
Both  fabrics  had  been  treated  with  a  high  friction  finish.  „  . 
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velocity  (v#)  from  each  layer  is  calculated  from  the  single  layer 
relation  (Esf(vj))  of  energy  loss  against  impact  velocity  (v^).  Thus 

2 

ve  s  i  oVj  -  f  (Vj) 

where  m  *  projectile  mass. 


This  velocity  is  then  taken  to  be  the  impact  velocity  onto  the  next 
layer.  This  iteration  is  continued  until  penetration  of  all  the 
layers,  or  until  the  projectile  is  arrested.  This  routine  is  repeated 
for  initial  impact  velocities  in  the  range  0  to  550  os*1  in  order  to 
obtain  the  multilayer  related  (E  *  f'(v^)}. 

This  routine  was  applied  to  single  layer  data  which  had  earlier 
been  obtained  for  nylon  and  Kevlar  treated  with  high  and  low  fabric 
finishes.  Experimental  data  from  the  multilayer  tests  was  then 
compared  to  that  predicted  by  calculated.  The  original  data  points, 
the  corresponding  least  squares  curve  fits,  together  with  the  predicted 
curves  of  energy  loss  against  impact  velocity  are  shown  in  figures  (11) 
to  (16). 

It  can  be  seen  that  generally  the  experimental  results  show  that 
the  maximum  energy  absorption  occurs  at  the  limit.  In  addition  the 
calculated  energy  losses  (shown  by  a  broken  line)  are  generally  higher 
than  those  found  experimentally  at  the  high  velocities.  The  one 
exception  to  this  is  the  Kevlar  treated  with  Siligen  E,  where  the  three 
layer  curves  are  similar. 

The  above  results  are  replotted  as  graphs  of  projectile  energy 
loss  at  the  V^q  limit  against  the  combined  area  density  of  the 
multilayered  sample.  These  are  presented  in  fig,  22  to  27.  It  is 
concluded  that  in  the  range  of  these  tests,  the  calculated  V^q  limit 
agrees  with  the  experimental  results. 
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General  Results  and  Conclusions 


A  ballistic  rig  has  been  designed  and  constructed  by  the  author  in 
order  to: 

1)  impact  1  gm  cylindrical  steel  projectiles  onto  a  fabric 
samples; 

2)  measure  the  impact  velocity  electronically; 

3)  measure  the  exit  velocity  photographically; 

4)  provide  single  or  multiflash  photographs  of  the  fabric  under 
impact; 

5)  measure  the  arrival  of  stress  waves  at  the  fabric  boundary. 
Test s  on  fabrics  have  produced  the  following  results: 

1)  the  shape  of  the  transverse  wave  front  is  octagonal  for  a  two 
ply  fabric,  if  the  yarn  directions  of  each  layer  are  set  at  45°  to 
each  other; 

2)  the  shape  of  the  transverse  wave  front  is  hexagonal  for 
triaxial  fabrics ; 

3)  two  triaxial  fabrics  have  been  shown  not  to  be  ballistically 
superior  to  biaxial  fabrics. 

4)  a  Kevlar  29  and  a  Kevlar  49  fabric  of  similar  weight  and 
construction  ’•'ere  found  to  have  similar  ballistic  resistances. 

The  following  conclusions  can  be  drawn: 

1)  the  octagonal  form  of  the  transverse  wave  in  the  bilayer 
impact  indicates  that  there  is  some  degree  of  interlayer 
interaction  and  that  any  computer  model  must  not  assume  a 
rhomboidal  form  as  a  starting  condition  unless  the  fabric  in  each 
layer  is  similarly  orientated; 

2)  the  hexagonal  shape  of  the  transverse  wave  for  a  triaxial 
fabric  indicates  that  more  energy  should  be  absorbed  as  kinetic 


energy  of  transverse  motion,  however  as  empirical  data  shows  that 
triaxial  fabrics  are  in  fact  inferior  to  biaxial  weaves  other 
effects  such  as  the  parting  of  yarns  at  the  impact  zone  must  be 
predominant; 

3)  the  high  friction  finish  probably  improved  the  ballistic 
performance  of  the  nylon  fabric  by  a  reduction  in  lateral  mobility 
of  yarns  at  the  impact  zone  or  due  to  the  influence  of  the  degree 
of  yarn  to  yarn  friction  on  the  strain  distribution  during  impact; 

4)  a  lightweight  nylon  fabric  was  shown  to  be  ballistically 
inferior  to  a  heavyweight  nylon  fabric  and  although  these  results 
could  be  partly  attributable  to  the  stress-strain  relations  of  the 
constituent  yarns,  it  is  postulated  that  the  mode  of  penetration 
(plugging  in  the  case  of  the  lightweight  close  weave  fabric)  is 
also  a  contributory  factor; 

5)  a  high  friction  finish  conferred  a  greater  ballistic 
performance  on  a  nylon  fabric  than  the  same  fabric  treated  with  a 
low  friction  finish; 

6)  both  the  high  friction  ana  the  bonding  agent  proved  seriously 
detrimental  to  the  ballistic  performance  of  the  Kevlar  fabrics 
tested; 

7)  a  simple  energy  model,  designed  to  simulate  Impact  onto 
separate  layers,  agrees  well  with  empirical  results  of  the 
energy  where  the  number  of  layers  is  low,  but  overestimates  the 
ballistic  performacna  when  penetration  occurs  at  high  velocities; 

8)  the  reduction  in  ballistic  performance  of  Kevlar  due  to  the 
addition  of  high  friction  or  bonding  agents  could  be  because  the 
level  of  friction  is  already  idealised  in  the  natural  state, 
unfortunately  no  finish  was  found  to  confer  a  lower  friction  on 


Kevlar.  Another  possible  effect  is  that  the  add-on  weight  could 
influence  the  mechanics  of  impact  and  penetration  and  therefore 
mask  any  effects  of  friction; 

9)  as  the  energy  models  for  multiple  impact  do  not  consider 
interlayer  interaction,  yet  generally  agree  with  empirical  data 
of  V^0  energy,  it  seems  that  interlayer  interactions  are  second 
order  effects,  where  the  number  of  layers  is  small. 

It  is  proposed  that 

1)  no  major  work  on  the  ballistic  properties  of  triaxial  fabrics 
alone  should  be  undertaken.  However  the  use  of  triaxial  fabrics  in 
hard  composites  should  not  be  dismissed  as  in  such  materials,  yarn 
mobility  will  be  reduced  and  penetration  due  to  the  pushing  aside 
of  yarns  will  be  precluded; 

2)  the  multilayer  impact  models  should  be  compared  with  data  from 
multilayer  impact  where  the  number  of  layers  is  greater,  in  order 
to  determine  whether  there  is  a  need  for  more  sophisticated 
models. 


CALCULATED  FROM 
SINGLE  LAYER 
DATA 
(Joules) 
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Siligen  E 
(high  friction) 


Lurapret 
(high  friction) 


Nylon  Siligen  E 
(low  friction 


EXPERIMENTAL 

RESULT 

(Joules) 

62.8 


Lurapret 

3' 

(high  friction) 

Table  4.  VKn  Energy  at  an  Area  Density  of 
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Introduction 


Textile  materials  in  the  forms  of  fabrics,  webbings  and  felts  are 
widely  used  in  energy  absorbing  systems  such  as  automobile  seat  belts,* 
parachute  lines,  and  personnel  protection.  The  particular  problem 
reported  here  .  .  is  an  investigation  of  the  behaviour  of  single¬ 
layer  and  multiple-layer  assembles  of  woven  fabric  used  in  protective 
clothing  worn  in  the  vicinity  of  medium  velocity  (500  ms"1)  particle 
hazards. 

These  protective  clothings  are  made  by  sawing  together  many  layers 
of  fabrics,  each  c-f  which  may  have  been  subjected  to  various  surface 
treatments.  Usually  for  easy  mobility,  the  fabrics  are  not  bonded 
together  and  can  therefore  slip  relative  to  one  another.  The  materials 
most  commonly  used  in  the  manufacture  of  these  fabrics  are  Nylon  and 
Kevlar.  The  yarns  are  bundles  of  thin,  long,  fibers,  held  together  by 
friction. 

The  mechanics  of  the  behaviour  of  the  clothing  is  complex  and 
various  simp 11 cations  are  made  to  render  the  problem  tractable.  The 
projectile  is  taken  to  be  rigid,  blunt,  and  travelling  perpendicularly 
to  the  fabric  assembly  without  spinning.  Each  yarn  is  assumed  to 
completely  lie  in  a  plane,  initially.  Under  these  conditions,  it  i3  no 
longer  possible  for  the  projectile  to  penetrate  the  fabric  simply  by 
pushing  aside  the  yarns  in  its  path  without  breaking  them.  Tarn  rupture 
occurs  when  a  given  maxiumum  strain  is  exceeded  and  visco-elastic 
effects  are  not  investigated. 

Different  mathematical  methods  are  used  to  investigate  the 
mechanics  of  the  fabric  assemblies  with  the  aim  of  developing  numerical 
methods  for  evaluating  the  effectiveness  of  various  fabric  parameters. 


Before  the  details  of  these  methods  are  given,  a  brief  review  of 
publications  by  previous  investigators  of  this  and  related  problems,  is 
given  in  Chapter  1. 

In  Chapter  2,  a  variational  model  is  described.  In  this  model, 
only  the  out-of-plane  transverse  motion  of  the  fabric(s)  is  considered. 
The  shape  of  the  indentation  is  first  predicted  from  the  given  weave 
pattern  and  the  equation  of  motion  is  then  derived.  Next,  the  effect  of 
yarn  crossovers  (nodes)  on  the  transmission  and  reflection  of  the 
strain  waves  which  are  initiated  at,  and  propagate  away  from,  the  impact 
area  are  investigated.  The  well-known  equation  for  the  propagation  of 
waves  along  a  string  forms  the  basis  of  this  enquiry.  Assuming  rigid 
connections  at  the  nodes,  equations  are  derived  for  the  impedance  of  a 
node  for  any  fabric  geometry.  Numerical  examples  are  obtained  to 
compare  various  geometries. 

In  Chapter  3,  a  description  is  given  of  an  application  of  the 
method  of  characteristics  to  the  investigation.  The  characteristic 
equations  of  the  three  types  of  waves  travelling  along  the  yarns  (no 
twist)  are  obtained.  The  motion  of  the  fabric  assembly  is  obtained  by 
following  the  characteristics  entering  each  of  the  nodes.  A  quasi- 
linear  model  is  also  developed  in  which  only  the  out-of-plane  motion  is 
considered:  an  attempt  to  reduce  the  amount  of  computation  effort.  In 
each  of  the  moaels  based  on  the  method  of  characteristics,  it  is 
necessary  to  give  the  fabric  a  non-zero  prestrain  to  prevent  numerical 
instability.  \ 

The  application  of  the  finite  element  method  to  the  development  of 
different  models  of  a  fabric  is  described  in  Chapter  4.  Attempts  were 
made  to  investigate  the  significance  to  overall  fabric  structural 
behaviour  ofj  crimp,  inter-yarn  slip,  and  inter-fabric  slip. 

In  Chapter  5,  the  numerical  results  are  discussed  and  compared 
with  the  available  experimental  results. 
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Chapter  6  is  the  conclusion  and  some  suggestions  are  also  included 
for  those  interested  in  pursuing  the  investigation. 

The  impacting  particle  is  always  a  blunt  cylinder,  mass  1.0003g,  * 

diameter  5.56  mm. 


r- 
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CHAPTER  1  LITERATURE  SURVEY 


The  geometrical  structure  of  the  arrangements  of  yarns  in  a 
fabrics  significantly  affects  some  of  the  fabric  properties.  This  fact 
makes  a  brief  review  of  publications  on  the  geometry  of  fabrics  of 
relevance  here.  Pierced)  made  an  Important  contribution  to  study  of 
the  geometry  of  the  plain  orthogonal  weave  consisting  of  two  families  of 
yarns  interlaced-  at  right  angles.  The  yarns,  in  that  investigation, 
were  considered  as  thin  rods  of  circular  cross  section.  Allowance  was 
made  for  the  effect  of  the  finite  radius  of  the  yarns,  at  the  cross 
overs,  that  is,  crimp  was  included.  The  deformation  characteristics  of 
the  model  was  not  considered  in  his  analysis.  Since,  geometric- 
mechanical  models  have  been  developed  for  various  fabrics  08),  09). 
01ofson(2)  developed  a  model  and  analysed  the  tensile  deformation  of  a 
fabric.  Kawabata  et  al,  in  a  series  of  papers  (3),  (4)  and  (5), 
developed  the  * stereo*  model  and  applied  it  to  the  study  of  uniaxial 
extension,  biaxial  extension  and  of  a  plain  weave  shear.  They  showed 
that  the  compressive  strength  of  the  yarn  significantly  affects  its 
mechanical  properties  under  these  loading  cases.  In  the  case  of 
uniaxial  extension,  localised  bending  at  the  yarn  cross  overs  was 
introduced  to  prevent  straightening  out  of  the  yarn.  This  bending 
resistance  was  produced  by  the  frictional  resistance  against  inter- 
fiber  slip.  In  the  analysis  of  biaxial  extension,  perfectly  flexible 
yarns  were  assumed.  The  resistance  of  the  fabric  to  shear  was  predicted 
to  be  due  to  the  torsional  resistance  to  the  rotation  of  the  yarns 
relative  to  one  another  at  the  cross-overs.  The  magnitude  of  the  shear 
stiffness  was  therefore  dependent  on  the  reaction  between  crossing 


yarns 


In  another  publication,  01ofsson(6)  investigated  the  rheology  of  a 
frictional-elasto-plastie  model  and  derived  expressions  for  the 
extensional,  shearing,  bending  and  buckling,  and  creasing  modulil  in 
terms  of  the  yarn  bending  stiffness,  its  crimp,  its  extensional 
stiffness  and  the  effect  of  sliding  friction  between  the  yarns.  There 
are  many  publications  of  studies  into  the  mechanics  of  fabric  under 
static  conditions  by  the  researchers  associated  with  Leeds  University, 
Grosberg  and  Kedia(7)  reported  that  the  initial  extension  of  a  fabric 
can  only  be  analysed  successfully  by  including  the  bending  resistance 
of  the  yarns,  and  that  the  assumption  of  perfect  flexibility  leads  to 
large  discrepancies  between  theoretical  and  experimental  results. 
Mac Roy,  McGraith  and  McNamara (8 ) ,  and  Bearle  and  Grosberg(9)  have 
analysed  the  mechanical  behaviour  cf  knitted  fabrics.  Skelton  and 
FreestonC 10)  reported  that  the  effect  of  a  decrease  of  the  pitch  of  a 
fabric,  that  is  an  increase  of  the  pick  per  length,  is  greater  interyarn 
force  and  poorer  translation  of  yarn  strength  to  fabric  strength.  The 
efficiency  of  translation,  however,  remains  above  90$  in  most  cases  at 
both  low  and  high  strain  rates  for  all  conditions  of  finish.  Shanahan 
and  Postle(ll),  and  Hearle(12)  have  also  reported  analyses  of  the 
mechanics  of  knitted  fabrics. 

The  mechanics  of  trlaxially  woven  fabrics  was  investigated  by 
Skelton(13).  The  results  showed  that  the  use  of  this  weave  pattern 
improves  the  shear  strength  of  a  fabric  over  that  of  a  biaxially  woven 

i 

fabric  of  the  same  area  density.  It  has  subsequently  been  reported(l4) 
that  a  triaxial  weave  is  superior  to  a  biaxial  not  only  in  resistance  to 
shear  but  also  in  resistance  to  uniaxial  and  biaxial  extensions. 

Huang(15)  analysed  the  finite  biaxial  and  uniaxial  extension  of  a 
completely  set  plain-woven  fabric  in  which  the  yarns  have  a  non-zero 
flexural  stiffness.  The  flexural  resistance  sharply  decreases  when  the 
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change  of  curvature  exceeds  a  fixed,  limiting,  value.  This  reduction  in 
flexural  resistance  is  due  to  slipping  of  the  fibers  in  the  yarn.  The 
resrlts  show  that  the  stress-strain  relation  for  the  yarn  is  non-linear 
even  for  a  Hookean  or  linear  fiber  material. 

Genensky  and  Rivlin(l6)  developed  a  theory  for  the  deformation  of 
*  a  network  of  orthogonal  cords,  assuming  a  stress-free  shear  deformation 
model.  They  analysed  different  deformations.  Recently 

Christo ffer sen (177  published  an  analysis  of  the  mechanics  of  the  In¬ 
plane  deformation  of  fabrics.  The  fabric  was  treated  as  an  orthotropic 
plane  which  behaves  elastically  in  stretching  along  two  orthotropic 
directions  and  is  capable  of  stress-free  deformation  in  shear.  No 
slippage  of  the  yams  was  allowed.  Results  were  published  for  stress 
concentration  at  a  crack  in  the  fabric. 

Another  line  of  investigation,  different  from  all  those  reported 
above,  was  adopted  by  Leech (20)  in  the  analysis  of  nets.  Fourier 
representation  of  the  thickness  changes  was  used  to  obtain  expressions 
fort  the  natural  frequency  of  the  net,  the  conditions  for  travelling 
waves,  and,  the  ray  and  front  theories  for  propagation. 

There  have  been  few  publications  of  the  analysis  of  the  dynamic 
response  of  uncoated  fabrics  and  dense  nets  to  impact  loading.  Most 
publications  on  the  mechanical  behaviour  have  been  for  statical  loading 
while  there  have  been  reports  on  the  dynamics  of  the  response  of 
Individual  yarns  to  impact.  The  paucity  of  publications  in  this  field 
is  reflected  in  the  fact  that  only  2  of  the  278  references  cited  in  the 
review  of  the  mechanics  of  penetration,  by  Jonas  and  Zukas  (21),  are  on 
the  penetration  of  uncoated  fabrics.  Manse 11(22)  compiled  an  up-to- 
date  list  of  publications  on  the  response  of  single  layers  of  uncoated 
fabrics  and  has  reported  some  new  experimental  work  and  finite  element 
analysis . 


In  the  analysis  of  the  mechanics  of  multilayer  of  fabrics ,  new 
factors,  either  negligible  in,  or  unconnected  with,  the  analysis  of 
single  layers  now  have  to  be  considered.  These  include  the  Influence  of 
slip  between  layers,  the  frictional  properties  of  the  fabric  surfaces, 
the  pattern  of  laying  the  layers  and  the  significance  of  the  compressive 
stiffnesses.  A  relevant  investigation  is  the  study  by  Marom  and 
Bodner(23)  in  which  results  show  that  a  multilay  assembly  of  thin 
aluminium  plates  has  a  higher  resistance  to  ballistic  impact  than 
either  a  monolithic  plate  of  the  same  total  thickness  or  seperated  thin 
plates.  The  different  modes  of  shear  and  compression  are  cited  as 
causes  of  the  different  resistances. 


The  publications  on  the  response  of  composite  plates  to  Impact  are 
also  relevant.  Daniel  and  Liber (24)  reported  that  the  primary 
deformation  appeared  to  be  a  flexural  wave.  Crlstescu  et  al(25),  on  the 
other  hand,  concluded  from  experimental  results,  that,  even  for 
composite  plates,  the  projectile  energy  is  absorbed  mostly  by 
extensional  deformation.  While  it  mostly  occurs  after  he^vy 
delamination,  extension  is  the  major  energy  absorption  mechanism,  but 
that  the  delamination  serves  to  spread  the  deformations,  thus  involving 


more  of  panel  in  the  energy  absorption.  The  sequence  of  delamination  is 


also  reported.! 


I  FINITE  ELEMENT  ANALYSIS  OF  FABRICS 
A  bibliography  of  most  of  the  published  finite  element  analyses  of 
fabrics  are  listed  in  Ref  (22).  Henghold  and  Russel(26)  developed  a 


formulation  for  Ja  cable  element  which  is  directly  applicable  to  the 


modelling  of  individual  yams.  0zdemir(27)  modified  the  strain 


formulation  in  the  last-mentioned  model.  In  the  new  formulation,  the 


cable  length  at  any  time  is  specified  as  a  function  of  the  lengths  from 


o 
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the  origin  to  the  nodes.  The  formulation  avoids  the  inconsistencies  in 
some  formulations  which  lead  to  non-zero  strains  when  the  shape  of  the 
cable  distorts  without  a  change  of  length.  The  equations  of  motion 
obtained  by  this  formulation  were  shown  to  coincide  with  those  obtained 
by  the  simple  rod  elements,  for  two-node  elements.  The  3-node,  and  4- 
node  elements  tend  to  be  overstiff  compared  with  other  published 
results,  de  Lynch (28)  analysed  the  dynamics  of  both  yarns/fibers  and 
fabrics.  Rod  elements  were  used  to  represent  the  yarns  while  semi¬ 
annuli  membrane  elements  were  used  for  fabric  analysis.  In  the  fabric 
analysis,  circumferential  buckling  was  reported  to  occur  during  the 
propagation  of  waves  prior  to  reflection  at  the  boundaries.  The 
accuracy  of  the  modelling  of  buckling  was  said  to  need  further  study. 

Shanks  and  Leech (29)  published  results  on  the  Influence  of  weave 
pattern  on  the  transient  response  of  coarse  nets  and  cloths  to  Impact. 
Zero  crimp  was  assumed  in  the  model,  an  assumption  which,  while 
realistic  enough  for  coarse  nets  may  lead  to  inaccuracies  in  the 
analysis  of  dense  fabrics.  Stubbs  and  Fluss  (18)  developed  a  space- 
truss  model  for  a  plain-weave,  coated  fabric.  The  model  incorporates 
crimp.  Each  element  has  6  nodes  and  is  a  combination  of  7  straignt 
rods.  The  formulation  allows  fcr  ?arge  deformation  and  non-line  .r 
material  properties  to  be  investigated.  While  results  were  not 
published  for  dynamic  loading,  this  model  is  considered  here  as  it  is 
easily  used  in  such  an  analysis.  A  similar  model  developed  by  Torbe(19) 
Involves  shear  resistance  and  has  the  nodal  variables  at  the  yam 
crossovers  condensed  cut,  leaving  only  4  nodes  for  element.  This 
crucifix  element  was  used  in  the  analysis  of  coated  airballons.  It 
should  make  a  more  efficient  use  of  computer  storage  than  the  6-node 
element  of  Ref  18.  However,  the  latter  model  allowed  for  yarn 
flattening,  or  yarn  compression,  which  the  former  neglected. 


MEMBRANE  MODELS 


Leech,  Hearle  and  Mansell (30)  developed  a  membrane  model  to 
analyse  the  arrest  of  projectiles  by  pretensioned  cloths  and  nets.  All 
in-plane  motions  were  neglected  while  the  shape  of  the  front  of  the  out- 
of-plane  transverse  motion  was  based  on  earlier  work  by  Leech  and 
Mansell(31)  on  the  prediction  of  wave  fronts  in  orthogonal  fabrics.  It 
was  shown  the  shape  assumed  for  the  cross  section  of  the  indent  does  not 
significantly  affect  the  numerical  results.  Already  mentioned  is  the 
use  of  semi -annulli  membrane  elements  in  the  finite  element  of  fabrics 
by  de  Lynch.  There  was  no  comparison  given  between  the  numerical 
results  and  experimental  results  and  as  such  it  is  difficult  to  specify 
the  extent  to  which  the  wave  front  could  be  predicted  using  annulii 
elements.  Annulli  elements  may  be  most  appropriate  in  the  analysis  of 
the  dynamics  of  knitted  fabrics  for  which  experimental  results  show  a 
circular  front  for  out-of-plane  transverse  motion.  Oden,  Kay  and 
Fost(32)  published  results  on  an  analysis  of  the  non-linear  response  of 
an  Incompressible  membranes  using  constant  strain,  triangular,  elements 
with  a  Green’s  deformation  tensor  formulation  for  strain.  The  loading 
was  provided  by  a  constant  force  which  was  applied  for  a  short  time  and 
then  removed.  Benzley  and  Kay (33)  obtained  numerical  results  for  the 
vibration  of  a  pretensioned  membrane  and  reported  that  these  results 
compare  well  with  analytical  results.  The  nonlinearities  caused  by 
large  deformation  were  shown  while  the  derivation  of  the  equations  of 
motion  of  the  nodes  of  the  elements  involve  the  determination  of  Cauchy 
stresses.  The  formulation  allows  for  the  Initial  orientation  of  yarns 
to  be  specified  and  allows  for  their  rotation  during  motion. 
Computationally,  the  disadvantage  of  this  formulation  results  from  the 
fact  that  the  directional  cosine  of  the  yarns  are  defined  relative  to 
the  initial  shape  of  each  element,  and  therefore  has  to  be  recalculated 
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for  each  element  shape.  This  is  except  when  a  membrane  is  divided  into 
square,  rectangular,  or  similar  triangular,  elements.  Another  method 
of  accounting  for  the  rotation  of  yarns  was  adopted  by  Tabaddor  and 
StaffordO^)  in  their  study  of  the  vibration  of  a  cord-reinforced  tyre. 
The  constitutive  equations  were  defined  to  incorporate  the  relation 
between  yarn  rotation  and  strains. 

Leonard  and  Verma(35)  used  a  double  curved  element  to  analyse  the 
properties  of  a  Mooney-Rivlin  membrane.  Coons  geometry  was  used  to 
obtain  accurate  representation  of  points  on  the  curved  surface.  The 
paper  contains  a  useful  bibliography. 

The  only  study,  known  to  this  author,  of  the  use  of  three- 
dimensional  finite  element  analysis  to  model  a  single  layer  of  a  fabric 
is  that  by  Lloyd (36).  In  this  analysis,  20-Node,  three  dimensional 
parallelpiped  elements  were  used  and  the  deformed  shape  of  a  knitted 
fabric  was  successfully  predicted.  However  the  model  fails  at  large 
strains. 


MULTILAYER  SYSTEM  OF  UNCOATED  FABRICS 
The  only  publication  encountered  in  this  literature  survey,  on  the 
response  of  layered,  uncoated  fabrics  to  impact,  by  Flaherty  (37),  is 
based  on  experimental  vc.k.  The  results  show  that  the  order  of 
arrangement  cf  the  fabrics  affects  their  combined  strength,  at  least, 
for  the  2-layer  system  reported.  There  are,  however,  reports  of 
numerical,  and  even  analytical,  studies  of  the  responses  of  layered 
beams  or  plates  made  from  aluminium  and  wood.  Some  of  these  studies  are 
briefly  mentioned  later. 

Some  of  the  phenomena  such  as  slipping,  which  may  significantly 
affect  the  overall  structural  behaviour  of  a  multilayer  system  of 
fabrics  have  already  been  investigated  in  the  study  of  metal 


■  t  - 

i 

»  1  :■ 

■;  ‘  \ 
v 

*  V 

/  w 

i  V 

\  • 


\ 


A  V 


■  u< 


deformation  of  Wilkins (38).  The  method  used  in  that  study  to  model 
separation  is  also  relevant. 

Thompson,  Goodman  and  Vanderbilt (39)  studied  the  effect  of 
interlayer  slip  on  the  deflections  and  stressed  in  a  statically  loaded 
layered  system  of  wooden  beams.  The  slip  was  related  to  the  shear 
across  the  thickness  of  each  layer.  With  the  assumption  of  negligible 
friction,  resistance  against  interlayer  slip  only  occurs  when  the 
layers  are  glued  together.  Suzuki  and  Chang (40)  showed  that  interlayer 
slip  may  significantly  affect  the  overall  structural  behaviour  of 
layered  wooden  structures,  in  a  study  of  bonding  failure  during  the 
bending  of  a  laminated  cantilever.  Analytical  expressions  were  derived 
to  relate  deflection  and  loading  force,  and  the  movement  of  the  sllp/no- 
slip  boundary,  with  Increasing  load,  was  demonstrated.  It  was 
concluded  that;  interlayer  slip  due  to  debonding,  if  it  occurs,  will 
reduce  the  overall  stiffness  of  the  cantilever.  This  conclusion  can  be 
compared  with  that  of  Refs  (24)  and  (25)  on  the  significance  of  the 
delamination  process  in  spreading  an  impact  load  over  a  large  area  of 
composite  sheet. 

There  are  other  publications,  which  are  not  on  the  analysis  of 
fabrics,  but  which  are  relevant  to  this  study  because  they  suggest 
possible  approaches  to  the  finite  element  analysis  of  multilayers (41 ) , 
(42),  (43). 

The  study  of  fluid  motion  in  containers  also  involve  the  modelling 
of  slip.  An  outline  of  the  problems  and  prospects  of  the  methods 
currently  used  in  that  analysis  are  given  in  a  book  on  structural 
mechanics  which  is  edited  by  Donea  (44).  One  of  the  most  important 
sections  in  the  book  is  the  description  of  an  extension  of  Wilkins's 
work  on  the  concept  of  'slave*  line,  'slave*  nodes  and  'master'  elements 
which  are  defined  at  the  interface  between  two  materials  which  can  slide 


relative  to  one  other  during  motion.  Other  important  contributions 
have  also  been  made  to  the  modelling  of  slip  in  finite  element 
analyses (45)  -  (51). 


CHAPTER  2  SIMPLE  MODELS 


2.0  Introduction 

In  this  chapter,  two  aspects  are  considered;  firstly  the 
variational  model  used  by  Leech,  Hearle  and  Mansell  (30)  in  analysing  the 
response  of  an  orthogonally  woven  cloth  to  transverse  impact  is 
extended  to  the  case  of  an  impacted  triaxially  woven  fabric.  The 
equations  for  the  two  fabrics  are  then  used  to  predict  their  relative 
effectiveness  in  retarding  a  projectile. 

The  second  aspect  considered  here  relates  to  nodal  impedances;  the 
relationship  between  yam  angular  arrangement  and  the  rates  at  which 
different  incidents  disturbances  are  transmitted. 

2.1  Wave  front  Variational  Model 

In  this  model  only  the  out-of-plane  motion  is  considered,  with  all 
in-plane  yam  motions  neglected. 

Firstly,  the  shape  of  the  transverse  wave  front  is  predicted  and 
the  size  of  the  disturbed  area  expressed  as  a  function  of  time  and  a 
similarity  coordinate,  ^  ,  which  measures  the  perpendicular  distance 
behind  the  wave  front,  the  magnitude  of  at  the  wave-front  being  1. 
The  deflection  at  any  point  in  the  disturbed  area,  W(  ^  ,t),  is 
represented  by  a  trial  function  which  satisfied  the  kinematic 
conditions  at  the  point  of  impact  and  at  the  boundary  between  the 
excited  zone  and  the  quiet  zone. 

The  total  energies,  (kinetic  and  strain)  are  calculated  and 
Hamilton's  principle  introduced  to  obtain  the  dynamics  of  the  system. 
Kamiltion'3  principle  states  that  between  times  t1  and 
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(T-U)dt  ♦  \  1  Wdt  =0  2.1 

t, 

for  any  dynamic  system,  where  the  kinetic  variables  vanish  at 
and  t2*  T  is  the  Kinetic  co-energy  of  the  system  and  W  is  the 
virtual  work  associated  with  slippage  and  internal  friction.  In 
this  case,  an  ordinary  2nd  order  differential  equation  results'  and 
can  be  solved  numerically  using  Runge-Kutta  methods. 

2.1.2  Application  to  a  Triaxially  woven  Fabric 

Leech  and  Mansell  (3D  showed  that  the  wave  front  in  a  transversely 

impacted  orthogonal  fabric  is  rhomboical.  This  prediction  is  confirmed 

by  experiments.  The  equation  of  motion  they  obtained  is  given  by 

t»x  *  y  ,  where  t  is  time;  Cx  and  Cy  are  respectively  the 
C  C 

*  y 

wave  velocities  in  the  x-  and  y-directions.  The  triaxially  woven  fabric 
consists  of  three  families  of  strands  each  of  which  makes  a  60  degree 
angle  with  the  adjacent  yarn,  Fig  2.*]  (o)  •  In  the  first  instance,  the 
wave  fronts  separating  the  ’quiet*  region  from  the  transversely 
displaced  region  is  necessarily  deduced. 

Consider  the  fabric  shown  in  Fig  2.1  and  let  it  be  transversely 
Impacted  at  afly  point,  I. 

The  axes,  si,  s2,  s3,  are  as  shown  in  Fig  2.1  along  the  strand 
directions.  The  fabric  is  then  divided  into  six  sectors.  Now  consider 
one  of  these  sectors. 

On  impact,  disturbances  travel  along  la'  and  lb'  respectively.  At 
b',  the  disturbance  along  lb'  is  partly  reflected,  partly  diverted 
along  b*a'  and  partly  transmitted  along  b’b  ,  similarly  for  the 
disturbance  along  la'  when  it  reaches  a*. 

A  strand  parallel  to  s3,  such  as  a  *b*  is  disturbed  at  both  ends. 
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After  the  two  waves  have  passed  through  it,  It  would  become  horizontal 
and  unstrained  for  the  rest  of  the  motion.  Considering  only  the 
tranalssion  of  out-of-plane  transverse  waves ,  It  is  clear  that,  at  time 
t  after  signal  initiation,  the  stress  front  in  this  section  is  given  by: 
t  *  si  s2  -  - 

El  ♦  U 

similarly,  the  wave  front  equations  in  the  5  other  sectors  are: 


t 

9 

-  si 

• 

s2 

Cl 

C2 

2.2b 

t 

9 

s2 

♦ 

s3 

C2 

C3 

2.2c 

t 

* 

s3 

si 

• 

C3 

Cl 

2. 2d 

t 

9 

-32 

- 

s3 

C2 

C3 

2.2e 

.  t 

9 

-a3 

♦ 

si  ) 

C3 

Cl 

2.2f 

where  Cl, 

C2,  C3  are  the  transverse  wave  velocities  in  the  directions  of 

si,  s2,  and  s3 

respectively.  The  fronts  whose  equations 

are  given  in 

equations 

(2.2 

,a-f) 

form 

the  sides  of  an  hexagon. 

The 

predicted  indentation  is  shown  in  Fig.  10 a  md  is  confirmed  by  the 
experimental  photograph,  Fig.  \  O'* 


2.1.3  Equation  of  Motion  for  small  deformation 

In  Appendix  1,  the  derivation  of  the  equation  of  motion  for  an 
Impacted  trlaxial  fabric  is  given,  using  the  wave  front  predicted  in 
section  2.1.2.  The  equation  relating  the  displacement  at  the  impact 
point,  Wo  and  the  non-dimensionalised  time is  given  by: 
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Also i  equation  2.3  could  be  written  as: 


(1+r2)Wo"+2rWo’+  (|l3+U-I2)wo  s  0  2.5 

15 

and  equation  2.4  as 

( 1+r2)Hol,+2rWo *  +  I3+I1-I2  Wo  =0  2.6 

U 

Equation  2.5  is  only  different  from  2.6  in  that  in  equation  2.5  the 
term  13  Is  multiplied  by  2/3,  while  in  equation  2.6  it  is  multiplied  by 
1.  The  last  expression  of  the  left  hand  side  of  equation  2.5  Cor 
equation  2.6)  was  obtained  from  strain  energy  considerations  and  the 
factor  of  2/3  in  equation  2.5  represents  a  reduction  in  strain  energy 

I 

per  unit  area  in  the  triaxial  weave  case  due  to  the  fact  that  in  each  of 
the  six  slanted  faces  of  the  pyramid,  one  family  of  strands  remain 

I 

horizontal  and  unstrained .  The  material  utilisation  factor  for  the 

! 

triaxial  weave,  relative  to  the  biaxial,  in  absorbing  projectile  energy 

I 

by  strand  straining,  is  2/3* 

The  ratio  of  the  displaced  areas  is  1.3:1,  to  the  advantage  of  the 
trlaxially  woven  fabric,  which  is  less  than  1:2/3  or  (1.5:1),  the  ratio 
of  material  utilisation  factors.  It  is  reasonable,  therefore,  to 
conclude  that  more  strain  energy  may  be  absorbed,  in  a  given  time,  by  an 
orthogonally  woven  fabric  than  by  a  trlaxially  woven  fabric  of  the  same 
area  density. 

The  ratio  of  the  time  factors,  afc/  *1.30, 
since  ab  =  Vlad  and  =®si 


This  ratio  results  from  the  fact  that  in  a  given  time,  the 


transversely  displaced  area  is  larger  in  the  fabric  with  0-60°  -120° 
arranged  strands  than  in  another  fabric  of  equal  area  density  in  which 
the  yarn  directions  are  perpendicular.  It  is  reasonable  to  suppose  that 
the  kinetic  energy  of  the  triaxial  fabric  is  high  r  at  any  time  after 
impact  than  that  of  the  biaxial  fabric.  This  combination  cf  these  two 
generally  produces  a  difference  between  the  performances  of  the  two 
fabrics  of  equal  area  density,  one  with  orthogonally  woven  yarns  and  the 
second  with  triaxial  weave. 

2.1.4  Equation  of  Motion  for  Large  Deformation 

Equation  2.3  and  2.4,  above,  are  only  applicable  for  small 
strains.  For  large  strains  the  full,  non-linear  strain  expressions  are 
needed. 

It  was  shown  in  Appendix  1  that  the  equation  of  motion  for  the 
triaxial  fabric  then  becomes: 

(Ur2)Wo"-»2rWo,-»6tfo(k1b1(uEi  -  1  )(1  ♦£,) 

f?U(WQa/C1r)2) 


+  k2b2(1-g.2),  ~  C 
£»  (  +c2  - 


- 2 -  ) 

/Tw(Woa/C2r)2) 


k3b3(1i£3)(u£3  -  _ 1 

£3  /fu<Woa/C3r)2) 


1/3)*  0, 


2.7 


where  kl 


C1C3+C1C2 
*  C1C3+C1C2+C2C3 


9 


k2  -  C1C3+C1C2 

"  C1C2+C1C3+C2C3 
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(1+  r2)Wo"+2rWo'+6Wo(-~l- 

i 


2.10 


*  «*r 


Comparing  equation  2.9  with  equation  2.10,  it  is  clear  that, 
again,  the  triaxial  arrangement  of  the  strands  leads  to  a  larger  time 
factor,  at,  and,  at  the  same  time,  causes  an  underutilisation  of  a  third 
of  the  fabric  as  far  as  the  absorption  of  energy  by  strain  is  concerned. 

The  above  observations  suggest  that  in  a  multilayer  ply,  there 
might  be  some  improvement  in  performance  if  the  direction  of  the  yams 
in  two  orthogonally  woven  fabrics  are  set  at  45°  to  one  another  so  as  to 
gain  the  advantages  of  the  two  different  weaves. 

This  possibility  may  be  examined  in  details  by  more  accurate 


methods 
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2.2  NODAL  IMPEDANCE  METHOD 
2.2.0  Introduction 

The  absorption  of  impact  by  fabrics  is  accomplished  by  the 
propagation  of  the  input  energy  away  from  the  directly  impacted  zone. 
Energy  is  dispersed  into  connecting  members  in  the  fabric  via  the  node 
connections.  In  general  three  types  of  signals  are  generated  in  the 
yarns  of  a  fabric  on  impact.  These  are,  in  the  plane  of  the  yarns;  the 
transverse  signal  and  the  extensional  signal,  and,  perpendicular  to  the 
plans,  the  out-of-plane  transverse  signal.  At  the  nodes,  each  type  of 
signal  may  generate  one  or  more  of  the  three  types  of  signals  in  the 
connecting  yarns,  depending  on  the  fabric  weave. 

The  yarns  normally  undergo  large  deformations  and  consequently  the 
complete  non-linear  string  equations  would  be  employed  to  model  the 
response  of  such  systems.  The  string  equations  have  been  developed  by 
various  Investigators;  for  example  see  Ref  30.  They  form  a  non-linear 
hyperbolic  system  with  two  distinct  propagation  speeds,  namely,  the 
propagation  speed  of  extensional  (fast)  and  transverse  (slow) 
disturbances.  The  nonlinearities  arise  from  considerations  of  gross 
deformation  of  the  string  (yarn)  and  from  the  admission  of  nonlinear 
constitutive  relations. 

In  this  section,  the  propagation  of  disturbances  through  fabrics 
of  various  weaves  is  considered,  the  linearised  decoupled  string 
equations  are  employed  and  coupling  between  extensional  and  transverse 
waves  is  initiated  when  nodes  are  encountered. 

The  effectiveness  cf  the  weaves  is  assessed  by  comparing  the 
percentages  of  the  incident  disturbed  signals  which  are  reflected  at  a 
node. 


I 


22 


r 

where  A  is  the  string  cross-sectional  area,  T  is  the  pre-tension  and 
?(X  ft)  is  the  transverse  displacement  of  the  string.  The  equation  of 
notion  for  transverse  disturbances,  admits  characteristic  solutions  * 

V  *  Ft(X  -Cfct)  and  V  =  G(x  +  Cfet^ 

9 

where  Ct(a  /  (T/^>  A)),  is  the  propagation  speed  of  transverse  signals. 

The  material  v train  induced  by  the  transverse  motion  is,  for  small 
motions: 


2.3. 


Although  this  is  a  second  order  quantity,  it  will  propagate  ahead  of  the 
transverse  signal  with  the  extensional  propagation  speed,  C  . 


2.2.2  MODE  EFFECTS  NEAR  IMPACT  ZONE 

When  a  disturbance  propagating  through  \  yarn  encounters  a  node, 
only  part  of  the  signal  will  pass  through  the  node,  while  another  part 
of  the  signal  is  reflected.  The  attenuation  will  be  dependent  upon  the 
type  of  disturbance  (extensional,  in-plane  transverse  or  out-of-plane 
transverse),  upon  the  material  characteristics  of  the  yarns  meeting  at 
the  node,  and  upon  the  weave  pattern.  This  node,  as  well  as 
transmitting  part  of  the  signal,  will  cause  a  reflection  back  towards 
the  signal  source  and  a  diversion  along  the  branch  yarns  and  for  in¬ 
plane  motions,  a  change  of  mode  along  the  connecting  yarns. 

In  an  impact  zone,  energy  is  continuously  provided  by  the 
projectile  as  it  is  retarded.  Three  types  of  signals  are  propagated 
along  the  impacted  yarns.  The  signals  are  consequently  partly 
transmitted,  partly  diverted  to  other  yarns  beyond  the  impact  zone,  and 


partly  reflected  back  Into  the  zone.  At  a  rate  of  energy  input,  the 

t 

chances  of  yarn  failui*e  increases  as  the  percentage  reflection  of 
signals  back  into  the  directly  impacted  yarns.  The  influence  of  yarn 
arrangement  on  this  percentage  is  an  indicator  of  weave  efficiency. 

In  the  first  instance,  the  passage  of  an  out-of-plane  transverse 

* 

Input  disturbance  through  a  node  is  considered.  Only  two  weave 
patterns,  the  orthogonal  and  the  (0  -  60°  -  120°)  triaxial  are 
considered.  The  general  equations  for  the  passage  of  in-plane  signals 
are  then  derived  and  expressions  obtained  for  a  few  weave  patterns. 


2.J.1  Out-of-plane  Transverse  Signals 

The  discussion  refers  to  Fig.  2.4(a),  2.4(b),  2.5(a)  and  2.5(b). 

In  an  orthogonal  weave,  a  disturbance  along  IM,  the  incident  at  I, 
creates  signals  in  the  three  ether  yarns  connected  to  I,  that  is  in  IK, 
IL  and  U.  The  input  signal  in  IM,  f^(  £  -  Cxfct),  gives  rise  to  a 
transmitted  signal  fj(  £  -  Ctxt),  a  reflected  signal  g(  ^  ♦  Cfcxt)  and 
signals  in  the  cross  elements  h(  Cfcgt)  where  is  a  coordinate 
along  a  yarn. 

Referring  to  Fig  2.5,  it  can  be  shown  that 


!' 


and 

J 


1 

1  » 


^  e1  h1 

J  f1»  b1 


i'1 


1 


1 

In 


2.3.5 


2.3.6 


where  (  )*  ■  *2)(  )/^£  »  th®  al°Pe* 

The  sum  of  the  amptitudes  of  the  incident  signal,  f^,  and  the 
reflected  signal,  g^, 

*  *  jjfj  from  equation  2.3.5  t 

•  Leech  and  Mansell  (3D  showed  that  the  corresponding  equations  for 


a  biaxial  fabric  arei 
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^BD 

a 


.11,,  11  i  ,i  wi  i  , 

h  *  iV,  «  i  ■  i  fn  h  ■  irj,  a.3.7 

where  f2,  g^  and  h  are  the  transmitted,  reflected  and  diverted  signals 
respectively. 

Therefore,  for,  a  biaxial  weave,  the  sum  of  the  incident  and 
reflected  signals  is  given  by 

fj  ♦  g^  *  2.3.8 

Since  the  sum  is  less  for  a  biaxially  woven  fabric,  failure  due  to 

overstraining  might  be  expected  to  occur  at  a  lower  impact  velocity  in  a 

trlaxlally  woven  fabrlo.  A  higher  proportion  of  the  incident  signal  is 

reflected  in  a  triaxial  fabric  because  the  greater  number  of  yams 

poses  a  stiffen  constraint  to  out-of  plane  motions. 

The  partitioning  of  the  total  input  energy  (kinetic  and  strain), 

II,  into  that  transmitted,  It,  that  reflected,  Ir,  and  that  diverted 

into  the  branch  yams,  1^,  and  Ib2,  can  be  shown  to  be  given  by 

It  t  1  i  Ir  i  Ibl  »  Ib2  «  1  ... 

li  9  li  9  Ii  Ii  9 

Hence  the  energy  in  the  yarn  carrying  the  incident  wave,  assuming 

continuous  input  of  energy  *  Ii.  (Ip  ♦  Ij). 

Under  similar  conditions  for  an  orthogonal  fabric,  the  energy  in 
e 

yarn  *  Ii. 


2.3.2  HOPAL  IMPEDANCE  FOR  IN-PLANE  SIGNALS;  GENERAL  EQUATIONS 

The  general  equations  for  the  Influence  of  a  node  subjected  to  an 
ln-plane  incident  wave  along  one  of  its  connecting  member  will  be 
derived.  The  percentage  of  the  incident  wave  which  is  reflected  along 
the  disturbed  yarn  is  of  paramount  importance,  since,  as  mentioned 
earlier,  a  yam  arrangement  which  causes  an  high  percentage  of  the 
Incident  wave  to  be  reflected  would  produce  yarn  failur**  ' 
signal  level. 
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Consider  the  node  0  in  Fig  2.5  at  which  n  ♦  1  yarns  meet.  The 
pretensions  in  the  yarns  are  enough  to  prevent  compression  (or 
buckling)  from  oceuring. 

The  yarn  labelled  1  is  subjected  to  an  extensional  wave  0(x-C#t) 
and  an  in-plane  transverse  disturbance  V(X  -  C^t),  where  the  wave 
speeds  are  given  by:  C#*  Jjj-  and  Ct  A. 

The  reflected  extensional  wave  is  taken  as  £U(x  ♦  C#t),  while  the 
reflected  in-plane  transverse  wave  is  dV  (x  ♦  Cfct).  The 

fractions  and  d  are  dependent  on  the  yarn  arrangement.  The 
extensional,  and  the  transverse  disturbance  generated  in  a  typical  yarn 
are  respectively  taken  as  0^(x  -C^t)  and  7^(x  -  C^t)  where  the  wave 
speeds  are  given  by 


Compatibility  at  the  node  requires  that  (Fig  2.6) 

•  •  * 

.  0(1- p  )  cos<<L  ♦  V(1-d)  sin  «<L,  2.3.10 

resolving  along  the  yarn,  and  perpendicular  to  the  yarn: 

\  «  -0(1  -p)sin<*L«.  V(1-d)cos«iL  ,  2.3.11 

where  the  dot  denotes  differentiation  with  respect  to  time  and 
primes  denotes  derivatives  with  respect  to  x. 

The  equations  can  be  rewritten  in  terms  of  derivatives  with 
respect  to  axial  length,  x,  as  ,  / 

CeI.0L  )  cos* L*GtV 1  ( 1 -d )  sin*L 


Ct  VL  *  “Ce°,(1“£  )  3in*L*Ctv'(1-d)  cos«<L, 
or 

*  Ce0'(1-  &  )  cos®^.*  Ct  v'(l-d)  sin*  _ 
“CeL  r  L  “CeL  L* 


•  Ct  V'(l-d)  cos  «<  L  - 
CtL 


2.3.12 


Ce  o'(1-  P  )  sin  °<  . 
CtL 


2.3.13 


The  equations  2.3. ’0  and  2.3*11  holds  for  each  yarn  connected  to 
node  0. 

Proa  equilibrium  considerations: 
along  disturbed  yarn, 

EA0'(1  )  a  Z<slAlUl  003 “  TLVL  sln#<L)  2.3.14 

Substituting  for  and  from  equations  2.3.12  and  2.3. 13> 
equation  2.3 J4  becomes 

KAO’d  ♦  f  )  •  ZtB^cos  *  L(§JL0'd-  P  Jain  01  L+^v’d-dJsin^) 


-  TLsin<*  l(§|l  V'(l-d)  cos*  L-§|L  o'd-  p  Jsin*  L)J 

i  • 

The  terms  involving  the  fractions  0  and  d  are  all  put  on  one  side  and 

! 

the  final  equation  is: 


°“So<L  *  \ 

***  tL  j 

- ‘-V’ JjiBjjSik  *  fJjL>  | 

.  -EM  ■  .  t  .  Tl  §l  sin®  <L) 


-  r'  Z'Vl 


Ct  Ct  .  sin2  o<  L 

CeL  lL  CtL;  2 


2.3.15 


Similarly ,  consideration  of  equilibrium  in  a  direction  perpendicular  to 
the  disturbed  yarn  gives  the  following: 


TV'd+d)  *  ^  siipf  ♦  Tl7j[  C0S<K  l) 


2.3.16 
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Upon  substitution  for  and  V^, 


Tv’d-Kl)  ■£elV1»  ^  L<ui  eti.  ('-f  >““V 


♦  TLC03  <  L  (»’  gL(t-4)  cos  <<  t  -  gL  D'  <1-  f  >  >in  <*  L>  * 


Ce  „1 


Ct  V'  (1-d)  sin  ^ 
CeL 


The  final  form  of  this  expression  is 


>/  [T  *  §L  S1”2<K  L  *  TL  §L  «¥»*•<  L 


-  f  u1  (EL^ 


Ce 

CeL  “  TL  CtL 


Ct  sin  2«<L 


*  -  TV1 


*  a'Z(E 


.  Ce  Ce  sln2  ^  L 
LaL  CeL  +  xLCtL  2 


-  r’  Z(h  §El  0032 *  L  Vl  %  Sl”  V  J-3'17 

'  The  suimation  is  over  all  xnfe  yarns  connected  to  the  node,  apart 
from  the  initially  disturbed  yarn. 

Given  the  number  of  the  yarns  connected  to  a  node  and  the  angle 
with  a  fixed  direction,  equations  2.3.15  and  2.3.17  would  give  the 
values  of  the  fractions  of  the  proportions  of  an  in-plane  disturbance 
which  are  reflected  along  the  disturbed  yarn.  The  expressions  for  the 
general  solutions  of  the  two  equations  are  long  and  were  not  very  useful 
for  numerical  work.  It  is  bette  for  any  given  fabric,  to  generate  the 
values  of  the  summed  expressions  from  the  given  values  oto(  T^, 
and  A^,  and  then  to  solve  the  two  resulting  simultaneous  equations. 

Once  the  values  of  j2  and  d  are  evaluated,  the  level  of  the  signal 
(IT.  and  V.)  in  the  other  yarns  can  be  calculated  using  equations  2.3.12 

U  Id 


2.3.3  IN-PLANE  TRANSVERSE  DISTURBANCES 

The  equations  connected  with  this  case  are  obtained  by  putting  U*  * 
0  in  equations  2.3.15  and  2.3.17.  They  are: 

.tl>  ■ 

and, 

»'•  1  [*  *  it  *  TL§L“>J<*J 

.  -  IV'  .  V'  I  ttjj,  =o.2  *l  .  Vl  ^2*  j 

The  expression  in  equation  2.3.10  shows  that  an  incident  in-plane 
transverse  disturbance  generates,  in  a  branch  yarn,  an  extensional 
wave,  if 

sino(^;o,  that  is  if  the  yarn  in  question  is  not  parallel  to  the 
L 

disturbed  yarn. 

In  the  case  when  only  two  yarns  are  connected  together  at  the  mode 
(  °^  ^  *  180°),  Fig  2.7(a),  the  equations  reduce  to 
0  a  0  and 

T'-d[T  *  Tl§J  2.3.20 

If  the  yarns  are  made  of  the  same  material,  the  right  hand  side  of 
equation  2.3.20  is  zero  since  the  tensions  are  the  same.  This  is  as 
could  be  expected,  and  the  wave  would  pass  through  the  yarn  without 
reflection  since  the  two  yarns  are  equivalent  to  a  single  continuous 
yarn.  If  yarn  1  is  denser  than  the  disturbed  yarn,  that  is  C.  ^  r  L 

W  t  * 
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some  reflection  occurs.  The  value  of  d  increased  with  the  ratio  Of  the 
densities  until  the  case  of  a  rigid  (yam)  (Ct^  C^)  when  the  whole 
signal  is  reflected. 

On  the  other  haM,  a  wave  of  the  opposite  sign  is  reflected  in  the 
case  when  yarn  1  is  lighter  than  the  disturbed  yarn  (Ct^  Cfc).  In  the 
terminal  case  with  the  second  yarn  removed,  d«-l.  This  is  as  could  be 
expected  since  the  free  end  must  be  stress-free. 


2.3.4  EXTENSIONAL  DISTURBANCES 
The  governing  equations  are: 


Cw 

-EAO  +  O'  lOi^cos2*  JSSl  ♦  T,  §L  sin2*  L) 

£.0*  [eA  +  ^(B^A^cos2 L  Ct_  ♦  TL  Ce_  sin2o< L>j|  7  2.3.21 


P  •®l£<It*Lc2  "  TL  CtL5 


Ce  .  sin  2  *  L 


-  "'I'Viifc  *  2is4ii  . 


2.3.22 


Por  the  case  ^  i  *  180,  Fib  2. (a),  equation  2.3.22  is  trivial, 
while  equation  2.3.21  gives 


t  ■(-EA.Et*,gL)/(EA  ♦  E,A,  gj  . 


2.3.23 


For  the  same  material,  £  a  0,  no  reflection.  For  a  free 
end,  CeL^  Ce,  fj  *  -1.  At  the  other  extreme,  when  the  yarn 
carrying,  the  disturbance  is  connected  to  an  extremely  dense  yarn,  = 
1,  complete  reflection  occurs  and  any  incident  strain  is  doubled. 


'  / 
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2*3.5  NODAL  IMPEDANCE  OF  DIFFERENT  YARN  ARRANGEMENTS 

The  different  patterns  of  yam  arrangement  which  are  discussed  are 
shown  in  Figs  2.7  (a  -  e) 


(a)  ORTHOGONAL  PATTERN 
Joint  data 


*1  *  90°,  <*2  a  180°,  0(3  *  270°  . 

(Assume  equal  tensions  and  yam  materials) 
^aae  1>  V  s  1,  o*  =  0  (in-plane  transverse) 


T+2EA—+T  a  -  T  ♦  2EA  ~ 

Ce  Ce  ♦  T 


T  ♦  EACt 

Ce 


EACt 

TCe 

1  +  EACt 
TCe 


_  P  .  EACt 

*  TTT  »  wher®  P  a  T§r 


2.3*24 


^a38  2,  7  *  0|  D1  *  1  (Extensional  disturbance) 
^.EA  ♦  £  (EA  +  2T  §r)»  -  EA  +  EA  +  2T  £2. 


P  =  T. 


L  Ct 

EA  ♦  T_Ce  ’ 


_ 1 

P  +  1 


2.3.25 


/ 

/  / 
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s 


1  +  p 
1  +  2P 
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2EA  *  Tl§| 

«* 

*  O  ♦  P)2  .  JL 

1  +  2P  1+P 


2.3.28 


Case  8 

Incident  in-plane  transverse  signal,  V1  s  1,  O1  :  0 


d  ^T  +  EA  (i+i+O+i+J)  +  T  (i+l+1+l>i)J 


f  (EA  -  T  S? 


ct)  Cl-l+i-D 


a  —  T  ♦  T  (i+i+1+i+i)  *  EA  g—  (l+i+O+J+J) 


|jtT  ♦  2EA  gjj 


d  |  4T  ♦  2EA  ~  |  =  2T  +  2EA  ~ 


d  a 


2T  ♦  2EA— 
_ Ce 

4T  ♦  2EA^ 


a  1  ♦  P 
2  ♦  P 


2  ♦  P 


1 

1  ♦  P 


2.3.29 


(d)  LIMIT  PATTERN  WITH  LARGE  NUMBER  OF  YARNS 

In  the  limit  when  the  number  of  yarns  are  large  and  uniformly 
arranged,  the  summation  signs  of  equations  2.3*15  and  2.3.17  are 
replaced  by  integration  signs.  This  is  how  the  result  below  is 
obtained.  This  limit  case  is  practical  with  the  regular  ’fan’ 
arrangement  of  yarns  which  are  then  held  together  in  a  matrix. 


Case  9  Incident  extensional  disturbance,  0  *  1,  V  =  0 


+  (EA  ♦  Ti 


+  d.o 


*  .  EA  .  Tt[eA  .  I  jjf] 


EA(  IT  -t)  ♦  ri 
EA(  If  ♦  1)  >TTgf 


0.6819P  +  1 
1.31bP  +  1 


(IT  -1)  P  -t-TT  (1  +  p)  _ 1_ 

(■jf  ♦  1)  P  +tt  1  ♦  P 


2.3.3C 


0.6819  P  •*•  1 
1.318P  +  1 


Cl  ♦  P)  . 


1 

TupT 


Case  10  Incident  in-plane  transverse  disturbance  V  a  1,  0  s  0 


d  [t  ♦  ♦  t)  a  -  T  +  (T  +  EA  |J)7r 


T  (f  -  1)  +1TEA  g 
T  (1T+  1)  ♦'n'EA  jg 


r 

IT  f..,l  +  p 

IT 
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(0.6819  ♦  P)  (1  ♦  P) 

1.3ia  ♦  P 


2.3.6  ENERGY  PARTITIONING 


In  order  to  compare  the  4  different  yarn  arrangement  patterns,  the 
rates  at  which  energy  is  reflected  are  calculated. 

The  energy  dissipated  in  time  t  are 


Incident  energy  a 
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p  * 

Reflected  energy  a  Ctt  T  d^V1) 

t  j— 


1  i  2 

CtEAf  (O1) 


2  2 

The  percentage  energy  reflected  is  proportional  toj&  and  d  . 
2  2 

The  values  ofj)  ,  d  are  compared  for  the  different  yarn  patterns. 

2  2 

In  Fig  2.8,  the  values  of  and  d  are  each  plotted  against 
values  of  Vp.  The  range  of  values  of  ^  used  is  0  -  1.  This  range 
corresponds  to  values  of  p  from  <»°  to  1.0  and  the  range  of  strain  from 
0.0  to  1.0. 

Note  that  for  a  yarn  made  of  Hooken  material: 


2.3.32 


2.3.33 


It  is  the  region  of  low  values  of  prestrain,  £,  ,  which  are  of 
practical  interest.  This  region,  as  the  graphs  show,  is  also  the  region 
where  the  percentac >  of  incident  signal  which  is  reflected  is  highest. 

The  most  "efficient"  yarn  arrangement  is  the  orthogonal,  since 
this  arrangement  produces  the  lowest  percentage  of  reflected  energy 
both  in-plane  transverse  disturbance  and  extensional  disturbance.  At 
very  low  prestrains,  the  percentage  of  input  in-plane  transverse  signal 
which  is  reflected  in  case  8  is  less  than  that  in  case  (1). 

The  graph  3  for  the  different  yarn  arrangements  are  very  close  to 
one  another  at  low  prestrains,  for  in-plane  transverse  disturbances; 
which  means  that  the  arrangement  of  tbe  yarns  would  not  appreciably 
affect  the  performance  of  the  fabrics.  At  zero  pretension,  in  Ref 
(61),  a  very  low  reflection  coefficient  (approximately  0.01)  is 
estimated  for  typical  fabrics. 


The  optimum  condition  occurs  at  p  *  1,  or  whan  the  prestrain  equals 
unity.  Most  yarns  would  faij.  long  before  this  strain  is  attained  so 
this  condition  cannot  be  realised.  Even  if  such  prestrains  could  be 
attained >  the  practical  difficulty  of  designing  rigid  anchors  necessary 
to  maintain  the  fabric  under  the  high  strain  remains. 

2.4  SUMMARY 

In  these  simplified  analyses,  slip,  friction  and  geometric 
non linearities  were  not  considered. 

The  results  show  that  a  blaxlally  woven  fabric  may  fail  at  a  higher- 
impact  velocity  for  a  given  area  density,  than  another  fabrlo  of  a 
different  weave  pattern.  This  is  true  for  out-of-plane  disturbances, 
and  in-plane  disturbances. 


CHAPTER  3 


CHARACTERISTIC  THEORY 

3.0  Introduction 

In  this  chapter,  an  application  of  the  method  of  characteristics 
to  the  analysis  of  the  dynamic  response  of  densely  woven  faerie  Is 
described.  The  fabrlo  Is  analysed  as  a  network  of  the  yarns.  The 
expressions  for  the  kinetic  and  strain  energy  In  a  yarn  are  stated  and 
both  the  equation  of  motion  of  the  yarn  and  the  compatibility  equations 
for  the  yarn  crossovers  are  obtained.  A  technique  using  the  theory  of 
characteristics  is  then  employed  in  a  discrete  fashion  to  determine  the 
motion  of  each  crossover  or  node.  The  motion  of  the  fabric  Is  then 
defined  by  the  motion  of  each  node. 

The  method  has,  computationally,  advantages  over  a  finite  element 

model  with  consistent  mass  formulation  in  that,  for  the  same 

discretisation,  the  rank  of  the  largest  matrix  to  be  inverted  in  the 

2 

solution  by  the  method  of  characteristics  is  3  against  n  (where  n  is 
the  number  of  nodes  in  the  finite  element  model).  It  is  also  better 
suited  for  wave  front  predictions.  This  is  because  the  positions  of 
the  waves  fronts,  even  then  they  are  between  nodes,  could  be  determined 
directly  with  the  characteristics  equations.  Whereas  in  the  finite 
element  solution,  the  wave-front  is  not  usually  sharp,  but  smoothened 
and  spread  out. 

3*1  THEORY 

The  following  assumptions  are  made; 

1)  the  centreline  of  the  yarns  are  initially  straight. 
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II)  The  yarns  «r«  rigidly  fixed  at  the  nodes,  slip  and 
separation  of  the  yarns  are  not  allowed. 

III)  the  yarns  are  completely  flexible  and  the  energies 
associated  with  bending  or  twist  is  negligible. 


3.1.1  FTBRE/TARN  DYNAMICS 

For  a  hyperelastic  yarn  material,  there  is  a  work  or  strain  energy 
functional  V  which  is  a  function  of  only  the  extensional  Green  strain  E; 
the  dependency  of  this  functional  on  the  other  components  of  strain  is 
assumed  to  be  eery  weak  since  the  yarn  is  stiff  in  tension  but  weak  in 
flexure  and  shear. 


the  extension  Green  strain  E  Is  written  as 


£.  =  u‘  ♦  £(u'  +  v‘  +  w'  ) , 


3.1.1 


where  u,  v,  w  are  the  motions  of  the  yarn  in  the  extensional  and  two 
transverse  directions,  Fig.  3.1.  the  strain  energy  functional  is  thus 


W  =W(EL,»s  W(u',v>  w'). 


3.1.2 


The  kinetic  energy  of  the  yarn  is  simply  obtained  and  the 
Lagrangian  j^for  *  complete  yarn,  length  L,  density  and  rosette 
area  A  is  thus 


/ l^vVw1)-  Wa/*’ 


where  J  - 

is 


3.1.3 
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and  t  and  s  ara  time  and  position  along  the  yarn. 

Now  by  applying  Hamilton's  principle  to  each  yarn,  and  ignoring 
for  a  moment  the  end  conditions,  then  the  following  equations  hold  with 
the  span  of  the  yarns. 

^(l+u')vV-0(,W)wU-. o  . 

pv- ^0+uyu«- [g'v'V^jv".  p v'wV"  -o 

These  three  coupled  ron-linear  partial  differential  equations  form 
a  hyperbolic  set  with  three  pairs  of  characteristic  equations.  One  of 
these  palx*s  is  associated  with  the  propagation  of  extension  waves  and 

5 

where  T,  the  2nd  Piola-Kirohhoff  stress,  is  given  by  dW 

d£ 

and  d^f  «  dT  is  the  modulus  of  elasticity  of  the  yarn.  The  remaining 
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pairs  of  characteristics  propagate  with  the  common  velocity 

c**  Fir 

where  Is  the  propagation  velocity  of  transverse  waves. 

At  this  point,  a  set  of  parameters  is  substituted  solely  for 


propagates  with  the  velocity 


Ce.  = /(T+  )!? 


/ 


simplicity  and  ease  of  development; 
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let  p»U,  f=W  3,1,6 

*nd  Q-  ^aV^sW*  3.1*7 

Then  £  a  Q  *  $*+ 1^*)  ?  3.1.8 

and  the  equations  of  motion  become 

°>  3.i.9(.) 

fy  ^0*e)4>e‘ -  o>  3.1.9(b) 
P-  d$Q*b]Y&'-^AVV-{j$:y1+$r}Vl  *°»  3-1,,(0) 

and  e-p‘  S°»  3.1.9(d) 

( 1  -O,  3.1.9(e) 

*yp -r1  *0.  3.1.9(f) 

This  set  of  six  first  order ,  non-linear ,  partial  differential  equations 
will  be  referred  to  in  the  later  section  on  the  application  of 
characteristics . 


3.1.6 

3*1.7 

3.1.8 


3.1. 9(o) 
3.1.9(d) 
3.1.9(e) 


3.1.2  Mode  Conditions 

In  the  previous  section,  the  end  conditions  associated  with  each 
yarn  have  been  ignored.  They  will  now  be  considered  for  each  yarn 
intersecting  a  node  and  a  compatibility  condition  for  each  node  will  be 
applied. 


Firstly,  the  end  variation  for  each  yarn  can  be  shown  to  be 

sis>Sutvsv+B'*K'j 

applied  it  each  end  of  the  yarn. 


How  let  each  yarn  originate  from  the  node  under  consideration  and 
now  label  each  yarn.  The  natural  boundary  condition  then  becomes 

57  *" Sq;  +  SVc,  +  =  0 

“  Ja» o  3.1.10 


where  there  are  m  yarns  meeting  at  the  node. 


/  '• 


Now  let  the  joint  movement  be  u,  v,  w,  and  since  each  yarn  end  must 


respond  with  the  node,  the  following  compatibility  equations  result 


Oi  *  U,CO&t  +  V  SinX, 

3*1* 11(a) 

Vi  *  -  U  Sin  *1  ♦  VCO^  i 

3.1.1Kb) 

Wi  «  V, 

3.1.11(c) 

where ^  is  the  Initial  angle  orientation  to  the  frame  axis* 

Fig  3.3 

It  then  follows  that  the  constraint  equations  apply 

-  SucosXi  +  SVcm**;, 

3.1.12(a) 

Sv;  =  -  Su;  Ziifti  +  Sv  CoiKi , 

3.1.12(b) 

5wi  s  Jto, 

3.1.12(c) 

and  the  node  equations  then  become 


+  =  °>  3  1.13 

\s  l  »tol 

s^Mfi  *o9 

is.  i 

since  Jo,  JV  and  §" V  are  independent;  however,  since  they  are  also 
arbitrary  as  they  are  virtual  displacements,  it  must  follow  that  their 
multipliers  are  zero  and  thi  node  equations  then  become 

-hsuli}  =o, 

[ih'j&lll+tts.nXit  facosXi}  =0, 

fc-*  I 

If  there  is  a  point  mass  located  at  the  node  then  these  equations 
have  to  be  re-written  to  take  account  of  the  kinetic  energy  of  the  mass 
located  at  the  node*  They  now  become 
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&nXi  J  -  o  > 

H:|Aij-|{((+ec)^^i+^cosXiJ  *°» 

Mr-.z  Ai^Vi  =°- 

3.1.3  Solution  Strategy 

The  solution  to  the  fabric  dynamics  problem,  as  has  been  stated,  is 
defined  by  the  deflections  of  each  node  or  crossover.  In  woven  fabrics 
the  free  span  of  yarn  is  small  since  the  material  is  an  array  of  closely 
packed  yarns,  it  is  sensible  to  consider  the  node  equations  as  equally 
important  as  the  yarn  equations. 

Consequently,  the  yarn  equations  are  written  in  characteristics 
form,  and  more  specifically  as  backward  going  characteristics,  and 
meeting  the  node  all  at  the  same  time;  they  now  become,  in  difference 
form, 

3.  ui6(a) 

<Mpi  -(i+M/ty +M&  Aft  1  -  o  >  3.1.16(b) 

i^ifip;.  -U+oOAfj  t]:0.  3.1.16(c) 

Equation  3* 1.16 (a)  is  the  backward  going  extenslonal  wave  and  the 
difference  operator  is  applied  between  point  d  and  point  a  along  the  I 
characteristic,  Fig.  3*2  Equation  3.1.16(b)  and  3.1.16(c)  are  applied 
along  the  backward  going  transverse  characteristics,  II,  and  the 
operator  A  in  this  case,  is  the  difference  between  the  values  at  point 
d  and  point  b.  The  time  step  At  and  hence  the  point  d  i3  dictated  by 
tho  condition  that  all  of  the  characteristics  shall  originate  from  a 
point  in  the  free  span  in  the  yarn  and  not  from  beyond  the  neighbouring 


nodes. 


Given  the  point  d  or  the  time  step  At,  the  points  a  and  b,  the 
origin  of  the  relevant  characteristic  is  found  by  projecting  backwards 
from  point  d  with  the  appropriate  slopes  Cg  and  Ct;  however,  these 
slopes  do  change  since  they  depend  on  the  local  strain  and  hence  the 
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local  deflections.  An  iterative  process  is  then  used  to  account  for 
this  curvature,  the  procedure  being  as  follows:  First,  using  the 
conditions  at  the  node  (d* )  which  are  known  and  at  node  c,  an  estimate 
is  made  for  points  a  and  b.  The  values  of  strain  and  hence  propagation 
velocities  at  these  points  is  then  evaluated  and  averaged  with  those  of 
point  d.  These  average  quantities  are  then  used  to  find  a  new  estimate 
of  points  a  and  b  and  the  process  is  repeated  a  number  of  times 
(approximately  five)  to  ensure  a  reasonably  accurate  location  of  the 
points  a  and  b.  the  iteration  scheme  is  shown  in  Fig.  3.4  and  Fig.  3.5 
The  equations  3.1.16(a),  (b),  (c)  can  then  be  used  to  search  for 
the  six  unknowns  of  the  yarn,  namely  pi,  qi,  ri  and  01,  01,  *^i  at  the 
point  d.  However,  the  velocities  are  constrained  by  the  node 
compatibility  equation  to  the  global  node  velocities. 


That  Is 


Pi  « pCosKv  +  > 

t:-ps’i*0U  +  ^CosXi, 


3.1.17 

n-r. 

Thus,  at  each  node  there  are  three  velocity  unknowns  (p,  q,  r)  and  three 
deformation  gradients  (01,  01,  y  i)  for  each  connecting  yarn.  There  are 
three  backward  going  characteristic  equations  for  each  connecting  yarn 
and,  finally,  there  are  the  three  node  equations  to  complete  the 
solution.  j 

Two  approaches  for  the  incorporation  of  these  nodal  equations  are 
now  considered.  Firstly,  if  the  problem  is  solely  concerned  with  fabric 
dynamics  without  'any  impacting  masses,  the  equations  are  firstly 
differentiated  and  then  numerically  Integrated. 

-  The  differentiated  equations  are 

.2  Ai'Jg  l&feX ;  SMOwh  -fynlJO*  ft)  et 

i  •  2  Tel  —  ® 

•  o 


or  in  matrix  form 


6i  =  l9i,<t>i>fi}  • 


3.1.20 


this  matrix  equation  can  be  written  in  difference  form 

•  1  Mfej'  -x»!° 

where  S&e  superscripts  D  and  0  refer  to  the  positions.  Fig,  3.5. 


3.1.21 


Using  this  approach,  the  problem  is  then  well  posed  for  numerical 
evaluation  using  simple  matrix  algebra.  For  each  joint,  the  problem  to 
be  solved  is  the  following  matrix  inversion: 

e,  Ea  -  -  rm 

6,  o  -  —  S: 

* 

o  B,  -  -  -  f; 

I  I  -  -  - 

i 

I  I  —  —  —  , 

I  I  -  -  -  1 

where  the  matrices  A,  B  and  C  are  obtained  from  the  characteristic 

equations.  The  rank  of  the  above  system  is  3m  ♦  3  where  there  are  m 
connecting  yarns  into  that  Joint.  Each  joint  can  be  solved,  for  that 
time,  independently  of  the  neighbours,  and  a  march  through  the 
structure  will  complete  the  solution  for  one  time  step. 
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The  problem  above  can  be  reduced  by  noting  that 

AlP+Bl9;  =  £■'<■  -f or  is  1,2,  ...  m,  3.1.23 

““.±E c9;  =  P„  .  3-'-24 

Then'’5’'  Q~-  >  3.1.25 

“d  Mnc'  i.  Ei sr'CCi  -Aip  ]  =  JF  •  3-’-26 

Solving  tfieV  for  p  gives 

P  =  (i9A',A;)',[l(EiBiJC;-f  J. 

~  tel  L^< 


3.1.27 


The  problem  has  thus  been  reduced  to  that  of  finding  the  inverses 
of  matrices,  rank  3»  rather  than  matrices,  rank  3m  ♦  3. 


7.  . 


\\| 


In  the  case  where  there  are  impacting  masses,  the  differentiation 
of  the  Joint  equation  leads  to  rates  of  change  of  acceleration  and  there 
is  no  improvement  in  this  formulation;  instead  a  finite  difference 
formulation  of  the  joint  equation  is  used;  theyj:hus  beccme 

fc-fc 


formulation  of  the  joint  equation  is  used;  tneyjthus  beccme 

.  A.  J  j  j  —  fcsWX*}*  3.1.28(a) 

ffi.lyjefo Sm% * c>tM‘  3.1.250,, 

*«  J  ul  J  +  &  CofXt?  p 

<*  &  -  rc*  ^  s-Vswtic  3.i.28< 


3.1.28(c) 


—d  —  -xri  —  *  e»c  * 

and  in  this  case  the  matrix  form  is 

CtJp,-f£i^  =  DJPft+|’E«8:<  3. >.29 

A  similar  mode  of  solution  is  employed  to  that  previously  described  and, 


in  this  case, 


p  =  +f  J  •  3.1.30 

**  (S  1  4“ 

In  this  latter  approach,  improved  approximation  dW,  the  intermediate 

dE 

values  of  0,  0  and  and  the  propagation  velocities  CQ  and  are 
obtained  by  performing  a  two  pass  system  at  each  Joint.  The  first  pass 
finds  approximate  values  for  the  joint  motion  and  yarn  deformation 
gradients  at  point  .  These  are  then  used  to  improve  the  approximation 


3.1.30 


N. 
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for  the  above  and  these  Improved  values  are  then  used  to  find  better 
values  for  the  joint  movements  and  yarn  deformation  gradients. 

3.2  UNEAR  OPT-OF -PLANE  TRANSVERSE  WAVE  MODEL 


This  model  involves  only  the  out-of-plane  transverse  wave,  as  all 
in-plane  motions  were  neglected.  This  model  was  developed  to  overcome 
the  buckling  of  the  yarns  which  may  occur  when  the  equations  developed 
In  Section  3. 1  are  applied  to  the  analysis  of  a  fabric  under  low  pre- 
strain.  The  minimum  prestrain  above  which  no  yarn  of  a  fabric  will 
buckle  increases  with  the  velocity  of  Impact. 


The  model  is  applied  to  a  highly  stressed  fabric  where 


<iw,  45  O2*.*  ^W.  V2,.  or  Af.W2.  and 

dE  OE  3  2e  *  3e  a 


the  product  terms  are  small.  Equation  3.1.4(c)  then  reduces  to 

fw  -  T.w'  *  *  0,  3.2.1 


where  T  s  m  the  tension  which  is  taken  as  constant.  The  corresponding 
characteristic  equation  is:  • 


Aij/* 

3.2.2 

.  ct  •  /(T/f) 

3.2.3 

The  nodal  equillbrum  equation  is: 

f*| 

Mr  -  ny;  -O. 

3.2.4 

The  solution  is  obtained  with  the  strategy  described  in  Section 


3.1.3 
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3.3  QPASI-LINEAR  OUT-OF-PLANS  TRANSVERSE  WAVE  MODEL 

In  order  to  maintain  the  accuracy  of  the  linear  model  described  in 
Section  3.2  without  losing  much  of  its  simplicity,  the  condition  that 
the  tension  in  the  yam  remain  constant  during  motion  was  removed.  The 
relevant  equations  are: 


and  the  equilibrum  equation 
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CHAPTER  4 

FINITE  ELEMENT  MODELS 

4.0  Introduction 

The  theory  of  the  finite  element  method  has  been  described  in  many 
books  and  publications  and  would  therefore  not  be  repeated  here. 

Interested  readers  are  referred  to  either  the  book  by  Zienkiewiez 
(57)  or  the  book  by  Irons  and  Ahmed  (58). 

There  are  many  possible  approaches  to  the  formulation  of  the 
equations  of  state  of  the  finite  elements.  Some  Involved  only  one 
Independent  field  while  others  involve  two  or  more  independent  fields. 
The  variational  statements  used  in  the  work  reported  here  Involve 
displacement  as  the  only  independent  field. 

A  short  theoretical  background  to  the  equations  of  motion  is 
given. 

Four  different  models  are  described.  Each  was  developed  for  the 
investigation  of  the  significance  of  specific  fabrics  parameters  to  the 
overall  structural  performance.  The  four  models  are: 

1)  A  space-truss  model  -  with  crimp,  slip  and  shear  in  single 
layers. 

2)  A  membrane  model  for  single  layers,  as  well  as  multilayers. 

3)  A  combination  of  membrane  elements  and  three  dimensional 
elements  to  model  compression  and  slip. 


48 


4)  A  two-dimensional,  layered  membrane  model  which  allows  slip 
between  layers,  called  the  master-slave  layer  model. 


4.1  GENERAL  BACKGROUND  THEORY 


In  the  finite  element  approach  to  the  analysis  of  a  body  it3  volume 
is  divided  into  a  finite  number  of  subregions  called  elements.  Within 
each  element  (subregion)  the  unknown  variables  are  represented  by 
simple  functions  of  their  values  at  points  on  (or  within)  the  subregion 
called  nodes.  The  nodal  variables  are  calculated  in  a  manner  that 
leaves  the  corresponding  variational  statement  stationary. 

Consider  a  volume  V  in  a  space  defined  by  cartesian  coordinate 
system  (x^,  x2,  x^)*  Let  the  nodal  displacement  vector  (U1 , 

U2,  0^)  for  node  n  p 

The  general  field  displacement  U  =  ^  where  there  are 

P  such  nodes  and, 

the  matrix  i3  the  shape  function  for  the  nth  node. 

The  Cartesian  components,  £ij,  of  the  Green's  strain  tensor  are 


defined  as 


p. .  .  I  (  +  Vi*  .  Vfk  \  .  .  .  ,  .  . 

It  is  assumed  that  there  exists  a  strain  energy  density  function, 
,  referred  to  the  undeformed  state  of  the  body  which  depends  only  on 
the  components  of  strain  tensors?  ••  ••  £-). 

To  define  the  state  of  stress  in  the  body,  the  second  Piola- 
Klrchhoff  stress  tensor  is  used,  given  by 


'SiV  — 


_  df 


^  be-,K 

As  in  the  classical  theory  of  elasticity,  it  can  be  shown  that 


4.1.2 


4.1.3 


and  that 


Cb  »  -k  C.V • 


where  Ctkjl  are  tha  components  of  a  fourth-order  tensor  of  material 
properties. 


The  total  potential  energy  of  the  system  consists  of 

1)  the  strain  energy,  given  by  0  s 

2)  the  energy  due  to  the  body  forces.  Lfl*M 

is  mass  per  unit  undeformed  volume) 

3)  the  energy  due  to  surfaces  forces  (these  forces  may  be  non- 
conservative) 


«  X.O 

Applying  the  principle  of  virtual  displacement,  after  summing  up  for 
all  the  elements  of  the  body,  the  variational  expression  1st 

c*cto 

The  strain 

point  displacement  vector  as£  ■[&•)&>  •••  S^BjU  ,and  the  virtual 

nodal  displacement  SU  can  be  expressed  as  and  eqtn. 

4.1.4  as 

/•dints 

(  0  denotes  differentiations  with  respect  to  time) 

Since  the  expression  within  the  brackets  must  always  be  zero  for 
any  arbitrary  variation  of  displacement, 


SjfF-u^v/  -  X  -Su  =  0  • 


4.1.4 


field,  ”3  is  then  related  to  the  nodal 


Y /CTrIC.:  J..  .  C/iCJ..  \/  \ 


Equations  (4.1.6)  can  be  written  in  matrix  fora  as 

m 

the  column  veotor  F  is  dependent  on  the  strain  and  therefore  changes  as 
the  volume  deforms,  while  the  column  of  externally  applied  forces  R  may 
also  vary  with  time. 

The  original  mass  matrix,  [m],  obtained  by  integrations,  was 
consistent.  Thereafter,  all  the  non-diagonal  elements  in  a  column  were 
lumped  with  the  diagonal.  All  the  results  presented  later  were 
therefore  obtained  using  a  diagenalised  mass  matrix. 

# 

To  obtain  the  displacement  vector  at  a  timet^j?  ,  the  central 

2 

difference  scheme  is  used  given  by 


ut+Afc  *  -F][fl]~'+  We 

4.1.8 

utt4£  =  Ut(ut+4%  ii*) 

4.1.9 

4.2  THE  SPACE-TRUSS  MODEL 
Introduction 

Photographs  of  impacted  fabrics  show  that  yarn  strain  distribution 
in  an  uncoated  fabrio  may  vary  discontlnuously  from  yarn  to  yarn.  This 
evidence  points  to  the  possibility  that  better  results  might  be 
obtained  by  modelling  a  fabric  as  an  assembly  of  yarns,  appropriately 
arranged,  than  as  a  thin  ctontlnuum.  A  model  of  fabric  as  an  assembly, 


of  yarns  whose  relative  motions  follow  a  pattern  determined  by  the 
fabric  weave  is  described.  The  model  is  designed  to  allow  the 
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influence  of  the  initial  distribution  of  impact  load  can  also  be 
investigated. 

4.2. 1  KINEMATICS  AND  STRAINS 

The  basic  element  of  this  model  is  the  cable  element  which  may  be  a 
simple  two-node  rod  or  a  higher  order  element.  The  element  may  have  all 
the  nodes  lying  on  a  straight  line,  or  the  equations  of  its  centre  line 
may  be  described  by  a  specified  polynomial  (or  transcedental) 
expression.  By  the  appropriate  choice  of  the  shape  of  the  centre  line, 
crimp  is  modelled. 

Pig.  4.1  (b)  shows  an  arrangement  of  four,  2 -node,  straight 
elements,  to  reproduce  plain-weave  while  Fig.  4.2(b)  shows  the 
combination  of  two  3-oode  elements  to  model  the  same  weave.  Pig  4.2(c) 
shows  a  model  of  a  trlaxlal  weave,  and  Fig.  4.2(d)  shows  a  knitted 
fabrio  model.  The  nodes  numbered  (5)  and  (6)  are  actually  vertically 
afrove  one  another  but  were  separated  in  Fig.  4.2(d)  only  for 
illustration  purposes.  The  same  is  true  of  nodes  (7)  and  (8).  Only  one 
strain  is  of  interest,  that  along  the  yarn. 


Following  Ozdemir,  the  distance  along  a  yarn  from  the  origin  to  any 
point  is  expressed  as  a  function  of  the  distances  from  the  origin  to  the 
nodes.  The  values  of  the  Green  Strain  of  the  yarn  based  on  such  nodal 
length  is  reported  to  be  more  consistent  than  that  obtained  by  the  usual 
methods  (27). 

For  a  2-node  element,  the  equation  obtained  by  this  formulation 
coincides  with  those  obtained  using  the  usual  rod  formulation. 
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4.2.2  MODELLING  OF  SLIP  AT  YARN  CROSSOVERS 

The  discussion  will  be  restricted  to  a  plain-weave  fabric  vith 
orthogonal  yarn  directions. 

In  the  space-truss  model  assembled  from  3-node  elements ,  slip  wj 
allowed  at  yarn  crossovers. 

To  simplify  the  calculations,  it  is  assumed  that  the  tvo  nodes 
representing  the  crossover  points  lie  on  top  of  one  another  throughout 
the  motion  of  the  fabrio.  This  approximation  is  reasonable  in  cards 
where  yarn  pull-through  does  not  occur,  and  the  total  slip  Is  small 
relative  to  yarn  length. 

| 

At  the  beginning  of  a  time  step,  both  the  magnitude  and  direction 
of  the  frietlonal  force  opposing  slip  are  unknown.  To  overcome  this 

•  j 

difficulty  a  two-pass  step  was  adopted.  j 

i 

At  the  beginning  of  the  first  pass,  the  frictional  resistance  is 
set  equal  to  zero, and  ,  the  transverse  velocities  of  the  two  overlapping 
nodes  and  a  first  estimate  of  the  direction  of  the  relative  motion  are 
calculated.  The  frictional  force  is  put  equal  to  the  product  of  the 
normal  reaction  and  the  corresponding  coefficient  of  sliding  friction. 
In  the  second  pass,  the  frictional  force  is  applied  as  an  external  force 
each  node,  acting  along  the  slip  direction  estimated  above,  but  opposed 
to  it.  Better  estimates  are  then  calculated  for  the  magnitude  and 
direction  of  the  slip  and  hence  of  the  components  of  the  frictional 
resistance.  The  second-pass  is  repeated  a  few  times  until  convergence. 
See  Fig.  4.3 


Let  the  relevant  equations  of  notion  of  the  nodes,  assembled  as  in 


equation  (4iJ7)  be 
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♦ 

F5z  s 
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i  the  prefix  (*)  denotes 

the  Kth 

4.2.1 


The  significance  of  an  approximation  mentioned  briefly  before  will 
be  stated.  This  was  the  approximation  involving  in  putting  the  node 
(5)  and  node  (6)  of  each  element  of  the  fabric  vertically  above  one 
another  at  the  beginning  of  each  time  step,  despite  the  cummulatlve  slip 
during  the  earlier  time  steps.  This  approach  obviates  the  need  to 
determine: 


(a)  the  actual  points  where  the  orthogonal  yarns  cross. 

(b)  the  nodal  equivalents  of  the  normal  reaction  and  the 
frictional  resistance  which  act  at  the  points  which  would 
otherwise  have  been  calculated  in  (a). 

If  this  approximation  were  not  made,  then  the  value  of  the  normal 
reaction  would  also  have  had  to  be  recalculated  at  every  pass. 
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4.2.3  FABRIC  SHEAR 

Shear  deformation,  at  a  node,  Is  defined  as  the  relative  rotation 
at  that  node  between  a  yarn  and  an  orthogonal  yarn  which  overlaps  it  at 
that  node.  The  shear  resistance,  T,  is  related  to  the  relative 
rotation,  0,  and  the  normal  reactions  at  the  crossover,  Fz»  by  the 
equation: 

T  »  C,  ♦  (C2  ♦  C3.Fz)  fr  4.2.2 

where  Cj,  C2  and  are  material  and  weave  properties.  The  equation  is 
adopted  from  the  work  by  Kawabata  et  al.  and  discussed  in  greater 
details  in  the  next  seotion. 

4.3  THE  MEMBRANE  MODEL 
Introduction 

In  this  model,  the  fabric  is  treated  as  a  thin  membrane  made,  from 
densely  woven  yarns.  The  yarn  are  thin,  perfectly  flexible, 
incompressible  rods.  The  weave  pattern  is  modelled  by  the  "stereo 
model"  described  by  Kawabata  et  al.  and  shown  in  Fig.  4.4  and  Fig.  4.5. 


4.3.1 


KINEMATICS  AND  STRAINS 


The 

Exx 

Eyy 

Eyxt  Exy 


relevant  strains  are  (Exx,  Eyy,  Exy,  Eyx),  given  by: 


«  +  1  (7**\  2  [dv] 

dy  2  L\dyj  *  Vdy/ 


1  f  dy  dv^  1  f*dv 

2  (  dy  *  dxy  *  2  [Jdx  * 


p> 
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4.3.2  CONSTRAINTS  AND  GEOMETRY  OF  WEAVE 

In  order  to  satisfy  the  weave  pattern,  both  force  and  geometrical 
constraints  have  to  be  satisfied  for  each  element.  This  mathematically 
means  an  addition  of  extra  terms  to  the  variational  statement  of 
equation  4.1.4.  Fig.  4.6  shows  the  initial  geometry  of  the  centrelines 
of  the  intersecting  yarns  before  deformation  and  Fig.  4.7  show  the  same 
lines  after  deformation. 

Let  Exx  be  the  strain  in  the  plane  of  the  fabric  in  the  plane  of 
yarn  1t  and  the  Eyy  the  strain  along  the  perpendicular  direction.  These 
strains  are  calculated  using  equations  (4.3.1).  Let  e1  be  the  actual 
strain  in  yarn  1.  which  lies  in  the  xz  plane.  If  the  pitches  in  the  x 
and  y  directions  are  P1  and  ?2  respectively,  the  equations 

Lf(1  ♦  e1)2sin201  a  P^O  +  Exx)2  ♦  h,2,  4.3.2 

122C1  *  e2)2sin202  ,  p22  (1  ♦  Eyy)2  ♦  h22,  4.3.3 

arise  from  the  assumption  of  no  slip. 

Also  from  the  incompressibility  assumption 

tj(1  eos01  ♦  £2  (  |  ♦  e2)cos02 

*  (^qCosB^q  ♦t20cose2Q  4.3.4 

*  h1  *  h2  *  h10  *  h20 

Let  the  tensions  in  the  yarns  be  T1  and  T2  respectively,  then,  from 
consideration  of  equilibrum  at  the  crossover  (see  Fig 

T^cosO^  a  T2cos02  4.3.5 


/ 


a. 
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At  any  given  time,  Exx  and  Eyy ,  the  strains  in  the  plane  of  the 
fabric  can  be  calculated  using  equations  4.3.1.  then,  using  the  Newton- 
Ralphson  method,  equations  (4.3.2)  to  (4.3.5)  are  solved  iteratively  to 
obtain  e^,  e2  and  therefore  (e^)  and  T2  (e2).  The  values  of  and  T2 
can  now  be  used  for  the  integration  of  equations  (4.1.6)  to  obtain  the 
equations  of  motion. 

4.3.3  SHEAR  EQUATIONS 


.  Until  the  shear  force  is  increased  beyond  a  limiting  value, 
relative  rotations  between  the  yarn,  in  other  words,  shear,  would  not 
occur.  Beyond  this  limit,  the  shear  strain,  Exy  or  Eyx,  increases  in  a 
complex  way  which  depends  on  the  magnitude  of  shear  force  as  well  as  the 
reaction  between  the  yarns,  and  on  the  direction  of  rotation. 

The  model  was  developed  by  Kawabata  et  al  (3)  and  could  be  written 

as 

T12  *  ±  Ci  *  (C2  ♦  C^TjCOsQj)  Exy  4.3.6 

where  is  the  shear  force,  is  the  frictional  constraint,  while  C2 
and  represent  the  effect  of  elastic  behaviour  of  the  contact  area. 

4.4  THREE-DIMENSIONAL  MODEL 

Experiments  show  that  in  a  multilayer  of  fabric,  the  layers 
closest  to  an  Impacting  projectile  fail  first  even  when  the  layers  are 
identical  and  slip  is  negligible.  Such  a  sequence  of  perforation  cannot 
be  easily  predicted  if  the  layers  are  treated  as  thin  plates. 
Otherwise,  a  failure  criterion  which  involves  the  magnitude  of  the 
on  the  surface  of  each  layer  may  be  necessary.  Apart 


from  this  observation,  it  has  be  suggested  by  Marom  and  Bodner  (23)  that 
compressive  resistance  is  significant  to  the  overall  structural 
mechanics  of  layers  of  comparatively  thin  1100-H16  aluminium  sheets 
and  beams. 

The  three-dimensional  model  is  therefore  designed  to  investigate 
the  significance  of  the  yarn  compressive  strength,  and,  at  the  same  time, 
the  significance  of  inter-layer  slip. 

Under  impact,  some  layers  may  be  penetrated  and  separated  from  the 
forward,  unpenetrated  layers.  This  model  allows  such  separation  to 

occur,  there  is  also  a  check  for  when  a  previously  separated  layer  re- 

! 

establishes  contact. 

i 

I 

The  in-plane  resistance  of  a  fabric  is  provided  by  a  two- 

dimensional  membrane  element.  The  interface  between  two  layers  is 

i 

modelled  by  three-dimensional  parallelpiped  elements  which  also  provide 

.  I 

resistance  against  shear  and  compression.  Each  three-dimensional 

i 

element  has  zero  mass,  and  its  thickness  equals  half  the  sum  of  the 
thicknesses  of  the  fabrics  on  either  side  of  the  interface  it 
represents. 

t 

Pig.  4.8  shows  the  combination  of  elements  to  model  three  fabric 
layers. 

4.4.1  STRAINS  AND  STRESSES 

For  the  two  dimensional  membrane  elements,  the  relevant  strains 
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*ife  2  Hi) "  ‘fell 

2yy  »  —  ♦  iflSsf  *  JiSi)  *] 

dy  2  |\dyl  \^yj  [dyl  J 
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2  \  dy  dx/  2  I  dx  dy  dx  dy  dy  dxj 

Prom  experiments,  the  in-plane  stiffness  coefficients  of  a  fabric, 

along  and  perpendicular  to  the  weft  direction,  are  measured.  They  could 

be  expressed  in  a  relation  between  stresses  and  strains  as:-. 
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4.4.2 


During  the  arrangement  of  fabrics  into  a  multilayer  assembly,  the 
warp  direction  of  a  layer  may  be  laid  at  an  angle  0  to  that  of  the 
reference  layer.  This  angle  may  be  different  from  layer  to  layer.  It 
is  necessary  to  transform  the  constitutive  relations  so  as  to  relate  the 
force  in  each  layer  and  the  strain  in  the  reference  directions. 

The  terms  Dij  of  the  constitutive  relations  in  equations  4.4.2.  is 
then  replaced  by  Dij  (0),  which  are  given  by  the  expressions: 

1  h  h  p  p 

+  ^22  s*n  ®  +  ^®i2  +  s*n  ‘ ®  cos*® 

D^22  s  D^sin^e  ♦  D ^  cos*1©  +  C2D12  ♦  40^)  sin2®  cos2® 

®12  *  ®21  *  ®12  +  (D^i  +  Dg2  *  20^2  —  ^^33)  sin2®  cos  9 
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D22  ♦  2D33^  sin  9  cos® 

D33  *  °33  +  (D11  “  °12  ~D33^  sin^©  cos®  - 
D22  ♦  2D33^  ain®  cos^9 
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4.4.2 


STRAIN  AND  STRESSES  FOR  INTERFACE  ELEMENTS 


The  nodal  displacements  are  the  same  as  for  the  two-dimensional 
fabric  elements.  The  important  strains  are:  (Exx,  Eyy,  Ezz,  Exz,  Eyz) 
which  are  defined  as: 

■»  ■  a  *  ?[te)2  *  tel2  ♦  tel] 
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The  derived  stresses  which  are  of  interest  are  Sxx,  Syy,  Szz,  Sxz,  Syz. 
they  are  related  to  the  strains  in  equations  4.4.4  by  the  expressions  of 
the  form: 
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Where  the  stiffness  coefficients  C.^,  C22»  depend  on  both  strain  and 


strain  rates 


For  unbonded,  layered  fabrics,  C^,  C12,  C21,  Cjg,  C^1 ,  C^2 

are  made  very  small  and  put  equal  to  the  average  compressive  modulus 
of  the  fabrics  on  either  side  of  the  interface.  When  Ezz  is  positive, 
that  is  separation  occurs,  and  the  other  coefficients  are  set  to 
a  very  small  number  (lO*1*).  When  the  interlayer  is  compressed,  that  is 
Ezz  is  -  ve,  la  set  to  the  value  stated  earlier. 

The  stresses  are  integrated  through  the  volume  to  obtain  the  nodal 
load  vector  (Fx,  Fy,  Fz).  In  cases  where  the  layers  are  free  to  slip, 
the  resultant  in-plane  force,  (Fx  ♦  Fy  )*  is  adjusted,  if  necessary,  so 
that 

(Fx2  +  Fy2}*  ^  n  Fz  4.4.6 

where  is  the  coefficient  of  friction. 

In  a  bonded  interface,  the  stiffnesses  in  the  xx,  yy,  and  zz 
directions  are  set  equal  to  that  of  the  binding  resin. 

4.5  THE  MASTER-SLAVE  LAYER  METHOD 


4.5.0  Introduction 

The  model  iiscussed  in  the  last  section  involved  the  use  of  three- 
dimensional  elements  to  simulate  the  compressive  resistance  of  the 
yarns  and  the  slip  between  layers.  The  thickness  of  each  fabric  is 
however  very  small  compared  with  its  in-plane  dimensions:  0.5mm  to 
150mu  by  150mm.  Very  many  elements  are  necessary  to  keep  the  aspect 
ratio  of  the  three-dimensional  elements  low  and  to  ensure  accurate 
stress  results.  This  requirement  drastically  restricts  the  number  of 


layers  which  could  be  modelled.  A  solution  to  the  problem  Is  to  avoid 
the  use  of  three  dimensional  Interface  elements  altogether.  The 
master-slave  layer  method  Is  an  attempt  to  achieve  this  end,  and  at  the 
same  time  allow  for  slip  between  layers. 

The  computational  procedure  Is  similar  to  that  used  In  fluid- 

structure  dynamic  analysis  to  investigate  the  motion  of  two  contacting 

materials  which  have  a  relative  tangential  motion  along  their 

interface,  but  which  remain  in  contact  throughout  this  motion.  See  Ref 

44.  The  material  lying  on  one  side  of  the  plane  of  contact  is 

designated  as  the  'master'  material,  while  the  material  on  the  other 

side  is  the  'slave*  material.  Generally  the  stiffer,  or  the  denser 

• 

material,  is  chosen  as  the  'master'  material.  Both  materials  are 
divided  into  elements  as  usual,  the  nodes  of  the  elements  of  the  master 
material  lying  on  the  interface  are  given  different  numbers  from  the 
nodes  of  the  slave  elements  lying  on  the  same  surface  even  when  some 
nodes  are  initially  coincident.  The  values  of  the  variables  at  the 
nodes  of  the  slave  elements  which  lie  on  the  interface  are  fixed  by 
those  of  the  nodes  of  the  master  elements  on  which  they  lie.  In  turn, 
the  slave  nodes  are  deemed  to  exert  forces  on  these  master  element;. 

the  investigation  of  the 
dynamic  response  of  a  multilayer  systemi  many  modifications  are 
necessary.  For  example,  the  master  layer  at  any  time  is  the  layer 
directly  in  contact  with  the  projectile  at  that  time.  Since,  during 
motion,  layers  may  get  penetrated  and  separate!,  the  layer  designated  as 
the  master  layer  must  be  upgraded  at  the  beginning  of  each  time  step  to 
be  the  layer  just  below  the  last  penetrated  layer.  To  satisfy  this 
requirement,  a  check  is  carried  out  at  the  start  of  each  time  step  for 
elements  which  have  failed.  If  elements  in  the  t^op  layer  have  failed, 


In  the  application  of  such  a  scheme  t- 
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then  the  projectile  is  assumed  to  be  in  contact  with  the  next  layer. 

The  check  is  continued  until  the  projectile  is  in  contact  with  an 
unpenetrated  layer. 


Once  the  master  layer  has  been  identified,  as  above,  the  next  step 
is  to  calculate  the  acceleration  during  the  time  step.  First  the  nodal 
loads  (Vector  F  in  equation  (4JJ6))  are  determined  for  the  nodes  of  each 
layer  not  yet  penetrated.  The  external  forces  acting  at  the  nodes  of 
the  elements  of  the  master  layer,  due  to  inertia  and  strains  in  the 
other  layers,  are  determined.  The  layers  are  assumed  to  be  constrained 
against  slip  arid  the  accelerations  of  the  master  nodes  are  calculated. 
From  the  accelerations  of  the  master  nodes,  an  estimate  can  be  obtained 
of  the  acceleration  of  a  node  on  any  slave  layer.  This  is  only  possible 
because  it  was  assumed  that  slip  did  not  occur.  The  reaction  force 
between  the  layers  is  calculated  at  each  node,  and  modified,  if 
necessary,  to  satisfy  the  frictional  properties  of  the  interface.  The 
acceleration  of  each  node  (master  or  slave)  is  recalculated  with  the 
modified  values  of  reactions. 


o' 


o 


o 


4.5.1  BASIC  PRINCIPLES  :  MASTER-SLAVE  LAYERS 


6 


Consider  the  membrane  elements  T  and  S  in  Fig.  4.9.  The  element  S 
lies  on  element  T. 


Let  the  element  T  be  the  master  element  and  the  element  S  the  slave 
element.  The  velocities  of  the  nodes  A,  B,  C  and  D  of  the  element  S  are 
therefore  determined  by  their  positions  on  the  element  T.  This 
statement  is  true  as  long  as  contact  is  maintained  between  the  elements 
without  slip. 


Q 


iCV 


The  equation  of  motion  of  the  nodes  of  element  T,  following 


equation  (4.1.17)  can  be  written  as 
”i  Si  ♦  h  «  Bj  • 

The  equation  of  motion  of  the  nodes  of  element  S  Is  similarly 
written  as 

Mi  ♦  Ji  *  h  t  4*5*2 

Where  the  externally  applied  force  vector  R1  is  assumed  to  be  only 
due  to  the  reaction  of  element  T  on  element  S. 

0 

The  velocity  vector  of  A,  B,  C  or  D  is  given  by 

-i  *  ^  *i  Hj  y  4*5.3 


where  is  the  value  of  the  shape  function  of  the  Ith  node  of  element 
T,  evaluated  at  the  position  of  node  1  of  element  S. 

Using  the  approximation 


»  -fc. 

U  x  U 


it 

+  A-fc.U  v 


avid  assuming  that  slip  Is  small  we  car.  rewrite  the  expression  (4.5.3)  in 
ter ks  of  accelerations  as: 


0 

-1 


r  »i  5.  . 


4.5.4 


64 


MiXhi  §!  ♦  fi  *  Ji  •  4.5.5 

By  virtual  work  considerations  It  can  be  shown  that  the  force  at 
node  1  of  element  T  due  to  the  action  of  external  force  Ri  at  A  is  given 
by 


5l  *  -  l*i  h  •  -  NlMiZNj3j  -  MiEi  •  *»-5.6 

The  total  external  force  ftj  at  node  I  is  the  sum  of  the  equivalent 
nodal  forces ,  Hj»  due  to  the  reactions  on  elements  which  share  node  I  :  t 
a  common  node,  that  is 

j?i  *  Z  j?2  ■  -  Hi  it  -  4*5*7 

Equation  4.5.1  can  now  be  rewritten  as 

MI  5l  ♦  Jl  ♦  NI  J  Nj  Oj  ♦  *  0  7  4.5.8 

the  summation  variable  L  refers  to  the  number  of  slave  nodes  which 
lie  on  master  elements  to  which  node  I  is  common,  while  the  variable  Q 
refers  to  the  number  of  nodes,  J,  of  the  particular  master  element  on, 
or  within,  which  the  slave  node,  i,  lies. 

By  this  procedure,  the  number  of  unknowns  has  been  reduced,  and 
only  the  acceleration  vectors  of  the  master  nodes  need  to  be  calculated. 
The  acceleration  of  the  slave  nodes  are  later  evaluated  using  equation 
(4.5.4). 

It  can  be  shown  that  the  expression  4.5.8  is  correct  even  when  the 
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slava  layer  (on  which  node  1)  lies  is  not  in  direct  contact  with  the 
■aster  layer,  but  is  in  contact  with  another  slave  layer  which  is  in 
eontact  with  the  toaster  layer.  The  expression  is  correct  fcr  as  long  as 
the  velocity  vector  of  the  slave  nodes  satisfies  equation  (4.5.3),  that 
is,  for  as  long  as  there  is  no  slip  between  the  fabrics  or  the 
separation  of  one  fabric  from  the  rest. 

# 

4.5.2  APPLICATION  TO  MULTI -LAYER  ANALYSIS 


The  layers  are  identically  divided  into  elements,  so  nodes  on  the 
lower  layers  are  vertically  belew  nodes  on  the  top,  master,  layer.  Once 
impact  occurs,  slip  generally  takes  place,  and  a  node  on  a  lower  layer 
nay  no  longer  be  vertically  below  any  node  in  the  top  layer.  Therefore, 
in  order  to  calculate  the  acceleration  vector  of  a  slave  layer  nods  from 
those  of  master  layer  nodes  (using  equation  4.5.3)  the  following  steps 
have  to  be  taken: 

1.  the  identification  of  the  element  in  the  top  layer  within  which 
horizontal  projection  the  slave  layer  node  lies. 

2. '  the  calculation  of  the  values  of  the  parametric  coordinates  of 

the  point  in  the  master  layer  element  (identified  in  1)  which  is 
vertically  below  the  slave  layer  node. 

4.5.3  IDENTIFICATION  OF  MASTER  LAYEB  ELEMENT 

BELOW  A  GIVEN  SLAVE  LAYER  NODE 

The  search  scheme  described  here  is,  strictly  speaking,  only 
suitable  for  elements  bounded  by  straight  edges.  It  can  however  be 


easily  modified  for  use  with  other  shapes. 


Referring  to  Fig  4.10,  the  problem  is  to  determine  whether  or  not  a 
given  point  N  lies  on,  or  within,  an  area  bounded  by  straight  sides. 
See  Fig  4.11,  also. 

The  nodes  of  the  elements  are  numbered  in  as  anticlockwise  sense. 

The  coordinates  of  the  rectangle  1234  which  envelopes  the  master 
element  is  found.  If  the  coordinates  of  N,  (Xp,  Tp),  do  not  satisfy  the 
conditions 

Xmln  Xp  £  Xmax 
and  Tain  4  Tp  Tmax 

then  the  search  procedure  is  terminated. 

If  both  conditions  are  satisfied,  then  the  distance  from  the  node  N 
to  each  of  the  sides  of  the  element  L7LK  is  calculated. 

If  point  (Xp,  Tp)  lies  within  the  element  whose  sides  have  end 
points  (Xj,  T.j)  and  (X2,  T2)  respectively,  then 

<T1  “  V*Xp  *  (X2  '  V  -*P  ♦  <XrT2  “  W  y  0  4.5.11 

If  this  inequality  is  not  satified  by  any  side  of  the  element,  the 
point  (Xp,  Tp)  does  not  lie  within  the  element.  The  condition  given  by 
equation  (4.5.11)  follows  from  the  choice  of  node  numbering  system. 


4.5.9 

4.5.10 
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Once  it  has  been  shown  that  a  node  lies  above  an  element,  the  next 
step  is  to  calculate  the  parametric  coordinate  of  the  vertical 
projection  of  this  node  on  the  master  element.  Thi3  is  the  object  of 
the  next  section. 


4. 5.* 


CALCULATION  OF  THE  PARAMETRIC  COORDINATES 
OF  A  POINT,  GIVEN  ITS  CARTESIAN  COORDINATES 


The  cartesian  coordinate  of  a  general  point,  A(X,  Y,  Z),  on  an 
element  can  be  calculated  from  the  nodal  values  using  the  expression 


»  .  £»<">  »(n)  <  j ,  n>, 


4.5.12 


where  the  nodal  shape  functions  N'n'  are  evaluated  at  this  point.  (Note 
that  the  slave  node  and  its  vertical  projection  on  a  master  element  have 


the  same  x  and  y  coordinates). 


The  parametric  coordinates  (  ^  I’Hl^  of  Point  1  ^XI,YI^’  are 


obtained  from  the  solution  of  the  equations 


~  *  £  x(n)  NCn)  ( £  , r^) 


4.5.13 


Yj  *  £y(n)  N(n)  (J  ,  q) 

The  equations  can  respectively  be  rewritten  as 


4.5.14 


=  h-i*M  N  (5 


4.5.15 


and  F2  (^,  r])  =  Yx  -£y(n)  N(n)  (£  ,^  ) 


4.5.16 


The  solution  procedure  used  to  solve  functions  F1  and  Fg  are  based 


on  the  Newton's  method  for  non-linear  equations. 
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4.5.5  ORGANISATION  OF  COMPUTER  PROGRAM 


The  diagonalised  mass  matrix  is  stored  in  an  array  MASS. 

At  the  beginning  of  a  time  step,  the  coordinates  of  all  the  nodes 
(slave  and  master  nodes)  are  known.  The  vector  of  internal  nodal 
forces,  F,  due  to  strains  are  evaluated  by  numerical  integration  and 
.stored  In  an  array  FOR  (NNODS,  NVAR)  where,  NNODS  =  total  number  of 
nodes  in  the  active  system,  that  is,  the  nodes  of  all  the  layers  that 
have  not  been  penetrated  or  separated,  and  NVAR  *  number  of  variables  at 
each  node. 


Next  the  element  in  the  master  layer,  above  which  a  node  on  a  slave 
layer  lies,  is  identified  as  described  in  Section  (4.5.3) .  Let  the 
master  element  be  called  E.  The  contribution  of  the  slave  node  to  the 
equation  of  motion  of  each  of  the  nodes  of  element  E,  given  by  H^,  (see 
equation  4.5.6  of  Section  4.5.1)  is  then  summed.  The  coefficient  of  the 
acceleration  term  of  this  contribution  (first  term  in  bracket,  eqtn. 
4.5.8)  is  added  into  the  corresponding  position  of  an  Initially  null 
array j CONTR ( NNOLA,  NNOLA), where 

NNOLA  s  number  of  free  (unfixed)  nodes  in  a  layer. 

The  second  term  within  the  bracket  (eqtn.  4.5.8)  is  added  to  the 
corresponding  location  in  an  array  RHS  (NNOLA,  NVAR). 

Note  that  Fi  in  equation  (4.5.6)  is  the  transpose  of  the  terms 
stored  in  the  ith  column  of  array  FOR (NNODS,  NVAR). 


•-  . 


The  equations  obtained  after  s omming  up  the  contributions  of  all 
the  slave  layer  nodes  and  completing  the  summation  of  equation  (4.5.8) 


where 


LA  *  the  modified  mass  matrix, 
(equivalent  to  array  CONTR). 


no  longer  diagonal  but  sparce 


Q  *  Vector  of  the  'acceleration  of  master  layer  nodes 
(*  NNOLAxNVAR). 


vector  of  modified  generalised  forces  at  master  layer  nodes 


(a  NNOLAxNVAR). 


The  equation  is  then  solved  for  Q. 


4.5.6  NODAL  ACCELERATIONS  WITHOUT  SLIP  RESTRICTION 

This  last  part  of  the  procedure  involves  the  calculation  of  the 
accelerations  of  each  active  (both  master  and  slave  nodes)  with  the  no- 
s lip  restriction  relaxed.  The  acceleration  vectors  for  the  muter 
layer  nodes,  Q,  calculated  in  equation  4.5.17  are  used. 

Two  arrays  are  set  up:  RETOP  (NNOLA,  NVAR)  and  REBOT  (NNOLA,  NVAR) . 
the  first  array  contains  the  vector  of  forces  externally  applied  to  the 
top  surface  of  the  layer  under  consideration,  while  the  array  REBOT 
(NNOLA,  NVAR)  contains  the  forces  externally  applied  to  the  lower 
surface  of  the  array.  Initially  all  the  elements  of  each  of  the  two 
arrays  is  set  to  zero. 

The  steps  taken  in  order  to  calculate  the  acceleration  of  an 


vv 
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unfixed  node  in  any  layer,  starting  from  the  top  (master)  layer  are: 


Note:  The  master  layer  is  the  topmost  layer. 


Two  column  vectors  and^R.Jj  are  set  up.  [?t]  *  [v  v  **3  j- 

contains  the  external  reactions  at  the  node  due  to  the  effect  of  the 
layer  Just  above. 

r  r  T  PT 

a  {^xb*  Ryb*  RZb|  ia  the  extepnal  **«aetion  at  the  node  due  to  the 
reaction  between  the  layer  under  Investigation  and  the  one  below  it. 

^4 

The  total  external  force  at  node  x*0i’Wi  *  C-]  i'  "-5*18  t 


For  a  slave  layer is  calculated  using  the  expression 

[si  *  -  2[sj]»i  w) 


4.5.19 


where  Rj  are  force  vectors  read  from  the  array  REBOT  which  must  have 
been  earlier  calculated  for  the  layer  above  that  of  node  I.  Nj^  i3 

the  shape  function  of  node  J  with  respect  to  the  element  above  which 
node  I  lies.  Generally,  the  node  I  lies  within  one  element,  but  it  may 
lie  on  the  boundary  between  elements.  In  the  second  case,  node  I  could 
be  assumed  to  lie  on  any  of  the  elements  sharing  the  boundary. 


For  the  master  layer,  iRtjj  is  null  since  its  top  surface  is 


assumed  free. 


The  components  of  I  Rt  are  now  stored  in  the  I  column  of  array 


RETOP. 


71 


The  forces  at  node  I,  on  the  lower  surface  of  the  layer,  is  calcul 
ated  using  equations  4.1.17  and  4.5.3  are 


» 

For  a  slave  node  I,  the  acceleration  vector  is  determined  by 
using  equation  4.5.4  and  the  values  of  vector  9  which  were  evaluated  in 
equation  4.5.17. 

A  check  is  now  performed  to  see  if  the  reaction  vector 
satisfies  the  frictional  properties  of  the  interface  between  the  layer 
on  which  1  lies  and  the  one  above  it. 


The  directional  cosines  of  the  normal  at  I,  (1,  m,  n),  are 
calculated.  The  normal  component,  Rg  of  flttL  is  given  by 


"h 

S 

1*Rxt  * 

“•V  *  "•  V 

4.5.21 

and  the  components  of  the  in-plane  part  of  the  vector 

Rt  j  are: 

m 

Rxt 

-  1J« 

px 

• 

Ryt 

- 

*  1 

ry 

* 

t 

Rzt 

-  =•«, 

If 

the 

value  of 

the  normal  reaction 

»  V 

obtained  from  equation 

4.5.21 

is  negative, 

then  separation  has 

occurred  and 

the  program  is 

terminated. 

Else, 

the  satisfaction 

of 

the 

no-slip 

assumption  is 

checked : 


The  resultant  in-plane  reaction, 


Rr 


•/ 


2  2  2 
rx  *  ry  +  rz  ’ 
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4.5.22 


while  the  maximum  allowed  in-plane  reaction  is  yy  R^. 


^  /*»  then  no  adjustment  is  needed  to  the  vectorQttJ  ^  and 
there  was  no  slip  at  the  node  I  in  that  interval*  If,  however,  this 
condition  is  violated,  the  components  of  the  in-plane  reaction  are 
adjusted  to  the  values: 

rx/R r 

gM.ry/Rr  4.5.23 

£N .  r, /Hr 

The  corrected  reaction  at  the  node  I,  on  the  top  surface, T(,, 


becomes: 


[?']  IC 


R^.r^/f^r  ♦  l.Sjj 

Rjj.r^/Rr  ♦  m.Rjj 
|^,.rj/Rr  +  n.Rjj 


4.5.24 


The  acceleration  of  node  I  is  recalculated  with  the  vector!  I  ^  in 


equation  4.5.18  replaced  by  jR 


•  the  vector 
RETOP. 


[5t]  ic  13 


IC. 


th 


now  stored  in  the  I  row  of  the  array 


After  these  calculations  have  been  completed  for  all  the  free 
uodes  in  a  layer,  the  contents  of  array  REBOT  are  now  replaced  with 
those  of  array  RETOP  in  preparation  for  the  repetition  of  the  same 
calculations  for  the  nodes  in  the  next  layer.  This  completes  the 
solutions  for  one  time  step. 


r 


i 


* 
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CHAPTER  5 

RESULTS  ArtD  DISCUSSION 

5.0  Introduction 


In  this  chapter  the  numerical  results  obtained  using  the  models 
described  earlier  will  be  presented  and  discussed.  The  potential  use  of 
the  different  models  in  the  accurate  analysis  of  fabric  mechanics  will 
also  be  discussed  and  finally  the  results  will  be  compared  with  those  of 
other  investigations.  Of  paramount  interest  is  the  prediction  of  the 
velocity  of  impact  at  which  a  given  fabric  fails  *  for  a  given  projectile 
mass.  Many  models  have  been  proposed  for  fabric,  yarn  and  fiber 
failure,  some  of  which  are  based  on  the  strain  energy  absorbed  by  a  unit 
length  of  the  yarn  while  others  are  based  on  the  strain  history  of  the 
yarn.  In  Ref  (54)  and  (55)  models  of  load  shearing  between  the 
constituent  fibers  of  a  yarn  under  tension  were  proposed.  The  loading 
to  which  the  yarns  in  a  fabric  are  subjected  are  mainly  due  to 
extension,  although  as  suggested  in  Chapter  2,  failure  may  occur  due  to 
local  overstressing  at  a  node/  The  overstressing  may  be  due  to 

nodal  constraints  or  wave  reflection.  In  order  to  minimise  the 
complexity  of  the  analysis,  a  simple  failure  criterion  was  chosen.  This 
condition  Is  that  when  the  strain  at  a  point  in  the  fabric  exceeds  the 
maximum  value  allowed,  the  yarn  breaks.  The  magnitude,  of  course, 
varies  from  material  to  material.  A  fabric  is  penetrated  when  a  yarn  in  • 
the  region  of  impact  is  broken.  Results  were  only  obtained  for 
materials  whose  mechanical  properties  were  independent  of  strain  rates. 
This  decision  was  taken  because  earlier  reports  (22)  indicate  that  the 
significance  of  such  strain  rate  dependency  on  the  mechanical  response 
of  fabrics  is  minimal. 
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Only  two  yarn  materials  are  investigated.  They  are  Nylon  and 
Kevlar.  This  choice  was  based  on  the  fact  that  these  two  materials  are 
the  most  commonly  used  for  weaving  protective  clothing  and  sheaths. 
Different  types  of  nylon  yarns  were  investigated.  Some  had  linear 
stress  strain  relations,  while  the  complex  con3titive  relations  of 
others  are  represented  by  three  linear  sections.  Only  Kevlar  29  was 
modelled. 

Input  Data 

Projectile 

Radius  2.76  mm 

Mass  1.1  g 

The  projectile  is  assumed  to  be  rigid  and  the  effect  of  Its  shape 
was  not  modelled. 

Impact  Velocity 

Impact  takes  place  normal  to  the  plane  of  the  fabric  panel.  A 
velocity  range  of  100  m/s  to  500  m/s  was  investigated. 

Fabrics 

Size  The  in-plane  size  of  the  fabric  panel  investigated  in  each 

case  was  approximately  154  mm  by  154  mm.  This  size  was  chosen  to 

correspond  with  that  cl  the  specimen  used  in  some  fabric 

penetration  experiments  in  UMIST. 


The  thickness  of  the  fabric  was  assumed  to  be  at  0.5  mm  for 


purposes  of  numerical  work. 


Yarn  Materials: 


For  convenience ,  the  units  of  the  yarn  modulus 

$ 

Newton/strain/Tex . 


Kevlar  29 


Modulus  t  48.56 


Breaking  strain  t  3.23  per  cent 


Nylon 


Type  A 

Strain  Range  Modulus 


0.0 

-  0.05 

1.75 

0.116 

-  0.20 

4.16 

0.20 

-  0.24 

1.75 

is 


Breaking  strain  24.0  per  cent 
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Strain  Range 


Modulus 


0.0 

-  0.1 

0.629 

0.1 

-  0.1475 

1*14 

0.1475 

-  0.27 

1.43 

Breaking  strain  24.0  per  cent 


Fabric  Weight 


Fabric 

Area  Density 

Tarns  per 

cm 

Tex  of  Tarns 

Weave 

No 

units  gytnm2 

Nylon 

D305 

240 

11.4 

12.3 

103 

10.1 

Biaxial 

Nylon 

0322 

386 

9.8 

9.1 

181 

170 

Biaxial 

Kevlar 

0236 

351 

10.8 

9.2 

175 

175 

Biaxial 

Nylon 

221 

3.63  3.75 

3.77 

200  205 

201 

Triaxial 

Nylon 

483 

Biaxial 

77 


5.1  SIMPLE  MODEL  RESULTS 


Linear  Model 


The  strain  expression  used  in  the  linear  model  (appendix  1)  for 
fabrio  with  high  prestrain  is 


£  *  (1  p)  dWo  , 

ds 


5.1.1 


where  £p  is  the  prestrain!  Wo  the  transverse  deflection  at  the  point  of 
impact  and  s  the  length  along  the  yarn. 


It  can  be  shown  that  the  maximum  value  of  the  slope ,  dW,  occurs  * 

ds 

immediately  after  impact.  This  fact  can  be  observed  from  the  graph  of 
the  indentation  against  time.  Fig.  5.1. 


The  maximum  value  of  dW  *  Vjj  where  Vp  is  the  velocity  of  Impact 

ds  Cs 

and  Cs  is  the  wave  speed  along  the  yarn. 


Equation  5*1.1  can  therefore  be  written  as: 


*  (1  ♦£?)  V£ 
Cs 


5.1.2 


For  a  linear  Hooken  material 


density  and  E  is  the  modulus t  hence: 


t  Cs  i| where  is  the  yarn 


£  *  (i  ♦  £p)  vp 

/TUT 


5.1.3 
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In  other  words,  the  value  of  the  maximum  strain  developed  in  the 
yam  is  directly  proportional  to  the  impact  velocity  for  a  given 
prestrain. 


The  impact  velocity,  V^f  required  to  produce  fabric  (yarn)  failure 
at  strain  ft.  in  given  by: 

9 


£»S*CS 

(1  *;tp) 


5.1.4 


This  velocity  increases  as  the  yarn  material  modulus,  E, 

.  [' 

increases,  which  means  that  fabrics  made  of  strong  materials  are  more 
efficient . 


Fig.  5.2  shows  the  deceleration  diagrams  for  three  different 

fabric  densities  at  an  impact  velocity  of  150  m/s.  In  Fig.  5.3  the 

! 

initial  slopes  of  the  Indentation  time  curves  are  all  the  same  for  the 

i 

three  fabrics,  which  is  as  expected.  Fig.  5.4  shows  the  deceleration 

I 

diagrams  at  a  lower  prestrain,  0.002,  for  2  nylon  fabrics  and  one  of 
Kevlar.  As  expected,  the  deceleration  is  higher  for  Case  (3).  'Fig.  5.5 
shows  that  the  use  of  a  heavier  fabric,  Case  (2),  does  not  lead  to  a 
reduction  in  the  value  of  the  maximum  strain  developed  in  the  yarns, 
although  it  leads  to  a  smaller  value  for  the  maximum  indentation,  and  a 
shorter  projectile  arrest  time  than  case  (1). 

The  results  obtained  for  the  response  of  two  fabrics  of  similar 

area  densities,  Fig.  5.2  for  a  240  gm  orthogonally  woven  nylon 

—2 

fabric,  and  Fig.  5.3  for  a  221  gm  triaxially  woven  nylon  fabric, 
predict  that  the  orthogonal  weaving  pattern  provides  a  better 
translation  of  yarn  strength  to  fabric  dynamic  stiffness. 


Non-linear  Model 


The  expression  for  strain  valid  for  large  deflections  which  are 
expected  to  occur  in  fabrics  with  low  prestrains  is  given  by: 

<’*%>/('  *gf) 

i 

It  can  also  be  shown  that  the  ratio  /  jo  \  attains  its  maximum  value 

V  vvJ 

of  _£«  Just  immediately  after  time  t  a  0.  Note  that  the  magnitude  of 

the  sJope  of  the  line  from  the  origin  to  the  curve  at  time  t  is  ^o 

C  t.  S66 

3 

Fig.  5.1. 

Hence  the  maximum  strain,  £  max,  is  given, 


£~  *$)•') 

The  impact  velocity,  V^,,  above  which  failure  of  the  fabric  occurs, 
is  given  by 

Vf  =  fcp.E.m  £s)2/(1  +  £p)2  -  1)  5.1.7 


The  optimum  value  of  prestrain,  corresponding  to  the  stationary 
value  of  the  variation  of  V^.  with  respect  to  Sp,  is  given  by: 

_  =  w  -  (1  L  )/3w,  where 

r  3 


This  optimum  conditions  corresponds  to  a  high  prestrain  which 
would  he  higher  than  the  breaking  strain,  £,g,  for  most  yarn  materials 
and  is  therefore  not  practically  realisable. 

The  general  equation  for  strain  (eqtn.  5.1.6)  predicts  that: 

1  the  velocity  at  which  failure  occurs  is  independent  of  the  density 
(mass  per  area)  of  the  fabric.  This  follows  because  the  highest 
strain  occurs  immediately  after  impact,  that  is  before  part  of  the 
projectile  energy  is  absorbed  by  the  fabric  as  kinetic  energy, 

2  The  value  of  the  strain  is  independent  of  the  geometry  of  the 
weave.  This  is  true  for  any  weave  in  which  the  yarns  through  the 
point  of  impact  are  straight  (neglecting  crimp). 

3  For  equal  prestrains,  the  value  of  the  maximum  strain  is  lower  in  a 
fabric  made  out  of  the  material  of  higher  modulus  since  the  value 
of  Cs  is  higher.  Protective  clothing  are,  therefore,  more 
effective  if  made  of  high  tensile  modulus  materials. 

Pig.  5*6  shows  indentation  of  240  g/m  ,  orthogonally  woven,  fabric 
as  predicted  by  the  non-linear  model.  The  variation  of  Impact  velocity 
at  failure  with  prestrains  is  displayed  in  Fig.  5.7  for  Kevlar  and 
Nylon.  Very  low  values  of  prestrain  which  may,  somehow,  be  easily 
produced  in  cloths  are  used  in  the  calculations. 
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5.2  Characteristic  Model  results 

A  property  of  the  characteristic  solutions  described  in  Chapter  3 
la  that  it  can  only  be  succesfully  applied  to  the  analysis  of 
prestrained  fabrio.  At  low  prestrains,  buckling  occurs  and  the 
solutions  goes  unstable.  The  level  of  prestrain  required  to  prevent 
buckling  increases  with  the  inpact  velocity.  It  was  noticed  that  in 
cases  when  the  initial  velocity  of  the  transverse  wave  (which  is 
proportional  to  the  root  of  the  prestrain)  is  less  than  the  velocity  of 
impact,  compressive  strains  are  developed  at  some  points  in  the  fabrics 
and  the  solution  routine  falls.  The  results  which  are  discussed  are 
obtained  with  the  smallest  prestrains  necessary  to  prevent  buckling. 

The  effect  of  prestrain  is  indicated  by  the  differences  in  the 
graphs  in  Fig.  5.8.  The  results  obtained  by  the  finite  element  method 
(membrane  model)  almost  exactly  coincides  with  experimental  results 
(56).  Until  approximately  160^3  after  impact,  the  deflection  of  a 
fabric  with  1*  prestrain  obtained  by  the  method  of  characteristics  is 
indistinguishable  from  the  finite  element  results.  Zero  pretension 
was  assumed  in  the  finite  element  solution.  Since  fabric  failure,  when 
it  occurs,  takes  place  early  during  motion,  it  may  be  predictable  by  a 
method  of  characteristics  solution  with  small  prestrains. 

Results  obtained  by  using  this  method  was  successfully  applied  to 
predict  the  maximum  velocity  at  which  layers  of  nylon  can  be  impacted 
before  failure  occurs. 

Two  schemes  were  used  to  model  the  impact.  In  one  scheme,  the 
velocity  of  each  impacted  node  was  put  equal  to  that  of  the  orn * 
at  time  t  *  0,  and  th*  •*— -  / 


82 


i 


p 


the  deceleration  of  the  projectile.  This  corresponds  to  putting  the 
acceleration  of  those  nodes  equal  to  infinity  at  time  t  a  0. 

In  this  second  scheme,  a  period  of  acceleration  of  10^>s  was 
assumed,  and  jin  this  periodjnormal  velocity  of  each  node  in  the  impacted 
zone  was  put  equal  to 

*  Vp  sin  IT 

“2 

where  Vp  is  the  velocity  of  impact.  Beyond  this  period,  deceleration 

occurs.  The  second  scheme  was  successful,  except  in  the  high  velocity 

cases,  when  some  of  the  yarns  attached  to  an  impacted  node  failed  during 

* 

the  period  of  acceleration.  The  first  scheme  was  used  on  such  occasions 
to  determine  the  correct  residual  velocity,  and  thence,  the  energy 
absorbed  by  the  fabric  before  failure. 

Apart  from  the  high  velocity  cases  mentioned  above,  there  was  no 
difference  in  the  results  obtained  by  using  either  scheme  with  the  same 
input  data.  While  most  of  the  results  were  obtained  from  solutions 
involving  the  use  of  the  first  scheme,  others  were  obtained  uiing  the 
second.  The  solution  converges  as  the  number  of  nodes,  was  increased . 
The  results  reported  below  were  obtained  using  a  39  x  39  grid.  For  most 
fabrics,  this  means  a  representation  of  about  every  4  yarns  by  a  single 
yarn.  j 

Fig,  5-9  (a)  shows  the  sequence  of  the  deformation  of  a  386  g/m2 
nylon  fabric  on  impact  by  a  1.1  gm  projectle  at  150  m/s.  The  arrest  of 
the  projectile  took  place  in  262  ^us  during  which  the  maxJ-mm  strain  in 
the  yams  through  the  impact  zone  was  21 .75.  The  initial  fabric 


(t/  0.00001) 


) 


prestrain  was  3 $•  Is  is  evident,  Fig.  5.9  (b)  that  the  out-of-plane 
transverse  wave  reached  the  fixed  edge  of  the  panel  before  the 
projectile  came  to  rest.  This  can  also  be  deducted  from  the  velocity 
profile,  the  absence  of  a  second  'hump'  in  the  curve  of  the  strain 
history  shows  that  the  reflected  out-of-plane  transverse  wave  did  not 
arrive  back  at  the  point  of  impact  befcr~  the  projectile  was  stopped. 

i 

Fig.  5.10  shows  similar  graphs  for  a  projectile  velocity  of  160  m/s. 
The  maximum  strain  in  the  yam  for  an  impact  velocity  of  170  ms*1 ,  Fig. 
5.11,  is  quite  close  to  the  breaking  strain  of  24$  and  the  fabric  will 
be  penetrated  at  a  slightly  higher  velocity. 

I 

MOLTI -LAYERS 

In  order  to  model  the  dynamics  of  layered  fabric  panels,  by  the 

method  of  characteristics,  the  yam  and  fabric  properties  of  a  single 

j 

layer  are  multiplied  by  the  number  of  layers.  This  method  involves  the 
Implicit  assumption  that  the  deformation  of  each  fabric  in  a  panel  of 
layers  is  the  same  as  the  next.  Slip  is  presumed  not  to  occur  between 
layers. 

The  prestrain  required  to  prevent  yarn  buckling  during  motion  was 
found  to  increase  with  the  number  of  layers  in  a  panel.  This  makes  the 
analysis  of  the  results  more  complex. 

The  results  obtained  for  a  simulated  panel  of  two  layers  of 
351  g/a2  Kevlar  fabrics  are  plotted  in  Fig.  5.12  (a)  and  (b).  The 
impact  velocity  is  290  m/s.  The  two  layers  absorbed  33.0  J  before 
penetration.  This  amount  of  energy  is  close  to  the  experimental  results 
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obtained  in  UMIST  when  a  panel  of  these  two  layers  of  Kevlar  fabric, 
free  of  prestrain,  was  impacted.  With  a  prestrain  of  0.0025,  three 
layers  of  Kevlar,  modelled  as  described  above  absorbed  36. 1J  when 
Impacted  by  a  projectile  moving  at  300  ms**1.  Again  this  amount  of 
energy  falls  within  the  range  of  experimental  results  obtained  for  the 

three  layers,  under  zero  prestrain  at  impact. 

♦ 

The  results  obtained  from  the  model  of  a  panel  made  up  of  two 
—2 

layers  of  240  gm  Nylon  fabrics,  with  an  initial  prestrain  of  3t, 
impacted  at  235  ms'1,  predict  a  maximum  strain  of  26.8  per  cent  -about 
the  average  breaking  strain  for  the  size  of  the  yarns.  See  Fig.  5.13* 
For  a  panel  made  up  of  three-layers  of  fabrics,  a  maximum  strain  of  26 
per  cent  is  predicted,  with  the  same  prestrain,  when  the  Impact  velocity 
is  increased  to  245  ms”1.  Fig  5.i4  shows  the  profiles  for  the  out-of- 
plane  transverse  velocity  along  a  yarn  through  the  point  of  impact  at 
five  different  times  after  impact.  Also  shown  are  the  deformed  shapes 
of  this  yarn  at  these  same  times.  From  the  graph  of  the  out-of-plane 
transverse  velocity,  the  speed  for  the  out-of-plane  transverse  wave  can 
be  estimated,  the  wave  being  assumed  to  have  just  reached  the  point  on 
the  yarn  where  the  graph  becomes  horizontal. 

Some  results  were  also  obtained  for  multi-layer  panels  made  up  of 
higher  density  nylon  fabric.  The  maximum  strain  developed  in  a  panel  of 
three  386  gm  Nylon  fabrics,  due  to  an  impact  at  255  ms  was  18  per 
cent.  The  panel  was  under  a  uniform  prestrain  of  3  per  cent.  The 
combined  area  density  of  the  panel,  1158  gm  is  almost  half  as  much 
again  as  that  of  the  panel  in  the  last  case,  while  the  maximum  strain 
was  reduced  by  only  about  30  per  cent. 


Fig.  5.15  shows  the  shape,  at  different  times  after  impact,  of  a 
483  gm  Nylon  fabric  impacted  at  150  ms  .  An  <1  x  11  grid  was  used. 
Fig.  5.16  shows  the  velocity  variation  along  a  yarn  through  the  single 
impacted  node  and  the  shape  of  the  yarn  at  different  times.  It  could  be 
noted  that  the  out-of-plane  transverse  wave  did  not  reach  the  edge  of 
the  panel  before  the  projectile  was  stopped.  The  results  shown  in 
Fig.  5.15  are  for  the  case  in  which  the  projectile  was  fixed  to  the 
central  node  throughout  the  motion. 

Fig.  5.17  shows  the  response  of  a  966  gm  fabric,  (2x483),  to  a 
similar  impact.  The  arrest  of  the  projectile,  as  could  be  expected, 
occurred  earlier  than  in  the  last  case,  and  with  a  smaller  out-of-pane 
Indentation. 

Fig.  5.18  shows  a  typical  variation  of  boundary  force  with  time. 

QUASI-LINEAR  CHARACTERISTICS  METHOD;  RESULTS 


The  values  of  prestrain,  above  which  numerical  results  could  be 

successfully  obtained  for  Nylon  fabrics,  is  10  per  cent.  Below 

this  prestrain,  wild  velocity  fluctuations  occurred.  For  Kevlar 

fabrics,  prestrains  greater  than  0.02  per  cent  were  required.  This  fact 

makes  this  approximate  model  only  useful  for  the  analysis  of  highly 

prestrained  fabrics.  For  a  given  fabric  prestrain,  the  quasi-linear 

model  is  much  stiffer  than  the  fully  nonlinear  model.  For  example,  the 

—1 

former  predicts  that  a  projectile  fired  at  150ms  into  a  240gm  Nylon 
fabric  would  be  arrested  lOOyUs  after  impact  while  the  nonlinear  model 
predicts  that  a  much  heavier  fabric,  a  483gm  fabric  with  the  same 
prestrain  value  of  10J  would  bring  the  project  to  rest  only  140^3  after 


impact 
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Fig.  5.19  and  Fig  5.20  show  respectively  the  results  obtained 

-2 

using  the  quasi-linear  model  to  analyse  the  response  of  a  240  gm  and  a 

O  1 

386gm_<:  Nylon  fabric  to  impact  at  Isoms'"  . 


5.3  FINITE  ELEMENT  METHOD  RESULTS 

5.3.1  Introduction 

In  this  section,  some  examples  of  the  discretisations  used  in  the 
various  finite  element  codes  are  mentioned  and  examples  are  given  of  the 
results  obtained.  Comparison  is  also  made  with  a  few  experimental 
results. 


5.3.2  THE  MEMBRANE  MODEL 

Firstly,  the  results  obtained  using  constant  strain,  three-ncde, 
triangular  elements  are  described. 

Fig  5.21  (a)  and  5.21  (b)  show  the  two  patterns  of  descretisations 
Investigated.  Vith  the  pattern  in  Fig.  5.2  (a),  compressive  strains 

are  generated  in  some  elements  for  all  Impact  velocities  when  the  yarn 
directions  are  chosen  to  be  parallel  to  the  perpendicular  sides  of  the 
elements.  The  descretisation  pattern  was  therefore  abandoned  and  that 
Fig.  5.21  (b)  used  subsequently.  The  results  obtained  with  this  mesh 
was  then  compared  with  the  experimental  values  of  deformation  obtained 
in  Ref  (56).  Even  for  this  coarse  mesh  the  predicted  formation  results 
were  close  to  the  experimental.  The  need  for  a  finer  division  of  the 
panel  was  found  necessary  in  order  to  obtain  accurate  strain  values. 
The  model  was  found  to  converge  and  the  final  mesh  chosen  for  generating 
results  is  the  800  element,  441  node  mesh  shown  in  e 
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The  membrane  model  involving  constant  strain  triangluar  elements 
was  used  to  investigate  the  significance  of: 

1  Initial  distribution  of  impact  load 

2  Shear  resistance  and  jamming 

3  Multilayer  panels 


It  was  found  that  for  some  values  of  input  data ,  there  were  wild 
oscillations  in  the  value  of  strain  at  the  point  of  impact.  These 
oscillations  were  suppressed  by  the  addition  of  a  damping  factor  kMO  to 
the  assembled  equation  of  motion.  The  equation  (equation  4117)  then 
becomes 


5.3.1 


where  k  Is  a  fraction  chosen  to  lie  in  the  range  .01  -  .1. 

1  Two  starting  procedures  were  used.  In  the  first,  a  single  node 
which  was  considered  to  be  directly  hit  by  the  projectile  was  given  an 
initial  displacement.  In  the  second  procedure,  the  impact  was  shared 
between  5  nodes  each_ of  which  is  given  the  velocity  of  the  projectile. 


2  The  fabric  stiffness  against!\  shear,  as  explained  earlier,  is 
assumed  to  be  mostly  due  to  torsional!  resistance  against  the  rotation  of 
one  yarn  against  another  at  their  crossover.  By  changing  the  weft,  or 
warp,  pitch  it  was  possible  to  change  the  number  of  crossovers  in  a 
given  area  and  therefore  the  shear  stiffness  of  the  fabric.  Using  the 
351gm  ^  Kevlar  fabric,  and  also  the  386  gm“^  Nylon  fabric,  as  test 
fabrics,  both  the  pitch  and  the  tex  of  the  yarns  were  halved  (area 


density  maintained  constant).  The  results  did  not  show  any  appreciable 

change.  This  was  inevitable  as  the  shear  3tiffnes3  of  the  fabric  was 

very  low  compared  with  the  stiffness  in  tension  in  either  the  weft  or 

warp  directions.  In  fact  there  was  hardly  any  change  in  the  strain 

history  at  the  point  of  impact  when  the  shear  stiffness  was  removed. 

The  value  of  the  parameters  and  Cj  in  eqtn.  4.3.6  are  respectively 
o 

3.5  x  10  and  0.0002.  These  values  are  taken  in  Ref.  (3). 

For  values  of  shear  strains  less  than  -  0.174  rad  (10?),  the  value 
of  the  parameter  in  the  Kawabata  equations  relating  shear  torque  and 
shear  strain  was  set  to  1.943  x  10“^gm.  Above  this  value  of  shear 

4 

strain,  C£  was  increased  to  I.Ogm,  that  is  by  a  factor  of  about  5  x  10  . 
By  this  means,  it  was  possible  to  prevent  shear  strain  from  exceeding 
0.20.  However  the  system  as  a  whole  became  much  3tiffer  than  the 
physical  system.  This  method  of  modelling  jamming  was  therefore 

abandoned . 

3  To  model  the  response,  to  impact,  of  a  panel  made  of  a  number  of 
fabrics,  the  panel  is  replaced  by  a  single  layer  but  the  tex  and  the 
modulus  of  each  yarn  in  this  model  layer  is  put  equal  to  that  of  the 
yarn  of  a  single  fabric  multiplied  by  the  number  of  fabrics  in  the 
panel.  Of  course,  this  approach  is  only  possible  when  the  fabrics  are 
all  identical,  and  laid  one  over  another  with  the  weft  yarns  parallel. 

Fig.  5.23  (a)  shows  deceleration  of  a  projectile  fired  into  a 
483  gm“  Nylon  Fabric  at  150  ms“  .  The  experimental  results  are  from 
Ref.  (56).  There  was  also  remarkable  closeness  between  the 
experimental  and  the  finite  element  model  results  for  the  transverse 
deformation  at  the  point  of  impact.  In  the  experimental  3et  up,  the  top 
edge  of  the  fabric  is  fixed  while  the  lower  edge  is  loaded  by  hanging 
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weights.  The  prestrain  produced  in  the  yarn  is  estimated  to  equal 
0.0022.  In  the  finite  element  model,  the  prestrain  is  put  equal  to 
zero. 


At  the  same  impact  velocity  of  150  ms”1 ,  the  projectile  was  stopped 
by  a  351  gm  Kevlar  fabric  in  170^* s.  The  maximum  strain  developed  at 
the  impact  point  was  2.7  percent  135 yUs  after  impact.  The  maximum 
transverse  deflection  of  the  fabric  was  1.6  cm.  It  is  instructive  to 
compare  the  response  of  this  Kevlar  fabric  with  that  of  a  Nylon  fabric 
of  similar  density,  a  386  gm  fabric.  The  Nylon  fabric  brought  the 
projectile  to  rest  in  353^  s.  The  transverse  deflection  at  the  impacted 
node  was  2.97  cm.  The  maximum  strain  in  the  fabric  was  0.17  which  is 
about  70  percent  of  the  breaking  strain  of  the  yarns.  See  Fig.  5.24. 
This  maximum  strain  is  a  smaller  percentage  of  the  breaking  strain  of 
the  yarn  than  in  the  previous  case  where  maximum  strain  attained  was 
about  84  per  cent  of  the  breaking  strain  of  the  Kevlar  29  yarns.  It  can 
be  noticed  from  the  shape  of  a  quarter  of  the  fabric,  displayed  at 
different  times,  that  the  faster  in-plane  wave  was  reflected  from  the 
boundary  before  the  projectile  was  stopped.  Experimental  results  show 
that  the  shape  of  the  transverse  indentation  is  rhomboidal,  this  is 
reasonably  in  agreement  with  the  results,  especially  in  the  first 
200 JK  a  after  impact. 

Fig.  5.25  (a)  show  the  deformed  shape  of  two  layers  of  240  gm"2 
Nylon  fabric  impacted  at  200  ms"1.  Again,  it  can  be  noticed  from  the 
distortion  of  the  grille  drawn  on  the  fabric  that  the  in-plane  waves 
have  been  reflected  from  the  boundary  before  this  time.  The  strain  at 
the  pcint  of  impact,  Fig.  5.25  (b)  shows  the  second  'hump*  when  this 
reflected  wave  arrives  back  at  the  impact  point. 
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The  model  predicts  that  a  panel  of  two  layers  of  351  gm  Kevlar 
fabric  will  fail  if  impacted  at  270  ms"1.  This  is  an  underestimation, 
most  likely,  of  the  impact  velocity  at  which  failure  occurs,  as  a  single 
layer  is  penetrated  only  above  250  as"1.  However,  the  energy  which  is 
predicted  to  be  absorbed  before  failure,  34.7  J  is  close  to  the 
experimental  values  for  the  energy  absorbed  by  the  2-layer  panel  at  high 
velocities. 

Finally,  a  comparison  is  made  between  the  results  obtained  with 

the  membrane  model  Involving  constant  strain,  triangular  elements  and 

another  model  involving  8  mode  quadratic  elements.  As  will  be  shown 

later,  the  results  obtained  using  the  quadratic  elements  are  generally 

very  close  to  the  experimental  results.  Fig.  5.26  (a)  shows  the 

variation  of  normal  strain  with  time,  at  the  single  impacted  node,  using 

the  two  different  elements,  for  a  351  gm  Kevlar  fabric  impacted  by  a 

1.003  g  projectile  at  350  ms"1.  The  calculations  involving  rectangular 

. 

elements  were  stopped  immediately  the  yarn  exceeded  its  breaking  strain 
of  3.2  per  cent,  while  that  involving  the  triangular  elements  were 
continued.  There  was  appreciable  difference  between  the  variation  of 
strain  in  the  two  models.  In  Fig.  5.26  (b)  the  two  set  of  results  for 
projectile  displacement  were  almost  coincident.  These  results  typify 
the  general  behaviour  of  the  membrane  model  with  triangular  elements. 
While  it  predicts  the  central  deflection  of  the  fabric  with  reasonable 
accuracy,  the  strain  values  are  poor.  The  poor  strain  values  result 
from  the  requirement  that  the  strain  does  not  vary  within  an  element. 

A  second  disadvantage  of  this  model  is  the  spurious  oscillations 
which  have  to  be  supressed  by  the  addition  of  damping,  -'.-.s  advantage  is 


that,  relative  other  codes,  particularly  the  program  based  on  the 
method  of  characteristics,  it  uses  very  little  computer  running  time. 


5.3.3  CABLE  ELEMENTS  t  RESULTS 


Introduction 


The  results  obtained  using  either  the  2-node  element  or  the  3-node 
element  showed  oscillations  which  were  damped  as  described  in  the 
seotlon  on  the  membrane  model.  Only  few  results  are  therefore 
presented. 

The  2-node  elements  were  used  to  model  different  fabrics  of  size 
152  mm  by  152  mm.  A  standard  39  x  39  orthogonal  representation  was 
used . 


A  circular  piece  of  fabric  152  mm  diameter  was  modelled  using  the 
3-node  cable  elements. 

Single-point  and  multiple-point  Impacts  were  simulated. 

The  simulations  of  crimp  effect  was  found  difficult  as  the 
numerical  oscillations  increased  rapidly  when  the  magnitude  of  the 
crimp  is  Increased  from  zero. 

An  attempt  to  allow  yarn  slip  when  using  2-node  elements  was 
abandoned  due  to  too  numerical  instability. 

Fig.  5.27  shows  the  strain  history  at  the  centre  of  the  240  gm“2 
fabric  impacted  at  150  ms”*.  The  projectile  came  to  rest  305^3  after 
impact,  showing  that  the  model  is  stiff er  compared  with  the  membrane 


impact  of  the  two  layers  of  a  152  too  x  152  on  square  piece  of  Kevlar  29 
fabric.  In  either  model,  the  fabric  was  under  a  prestrain  of  0.005. 
While  the  values  for  the  central  out-of-plane  deflection  predicted  by 
the  two  methods  are  reasonably  close,  Fig.  5.28  (b),  there  are  large 
differences  in  the  predicted  strains.  It  however  has  to  be  remembered 
that  the  yarn  strain  plotted  in  Fig.  5.28  (a)  for  the  2-node  model  is 
for  the  maximum  strain  in  the  yarn  at  any  given  time.  This  maximum 
value  may  occur  in  the  weft  yarn  through  the  centre  at  one  Instant  and 
occur  in  the  warp  yarn  at  the  next  instant.  However  the  fact  that  no 
crimp  was  included  in  the  characteristic  model  ensures  that  the  equal 
strain  occur  in  the  four  yax-ns  connected  to  the  impacted  node. 

• 

The  results  shown  in  Fig.  5.29  (a)  and  Fig.  5.29  (b)  are  obtained 
using  the  3-node  cable  element.  Only  a  quarter  of  the  fabric  is 
modelled,  and  the  model  has  1592  nodes.  The  oscillation  of  the  strain 
values  could  be  noticed  in  both  cases.  Slip  was  allowed  in  the  first 
fabric,  and  a  coefficient  of  frictions  of  0.5  was  used.  The  percentage 
of  an  incident  wave  which  is  reflected  from  a  node  is  maximum  when  the 
yarns  are  rigidly  connected.  The  percentage  which  is  reflected  for  any 
given  value  of  friction  coefficient  lies  between  this  maximum  value  and 
zero.  This  is  why  the  strain  at  the  point  of  Impact  is  higher  when 
rigid  crossover  connections  are  assumed.  Tam  failure  and  fabric 
penetration  is  therefore  likely  to  occur  at  a  lower  impact  velocity  when 
the  yarns  are  rigidly  connected. 

On  the  whole,  the  rod  elements  are  stiff er  compared  with  either  the 
membrane  (triangular  element)  model  or  with  experimental  results. 


5.3.4  THREE-DIMENSIONAL  ELEMENTS:  RESULTS 


Fig.  5.30  shows  the  descretisation  used  in  a  model  involving  the 
use  of  3-dimensional,  8-node  parallelopiped  interface  elements  and  4- 
node  membrane  elements.  Only  a  quarter  of  the  fabric  is  investigated 
because  of  the  symmetry  of  the  problem.  Fig.  5.31  shows  the 
discretisation  used  in  another  model  in  which  the  4-node  membrane 
elements  of  the  former  model  were  replaced  by  9-node  rectangular 
elements.  The  details  of  the  interface  elements  in  this  second  model 
are  shown  in  Fig.  5.32.  There  are  four  interface  elements  covering  the 
surface  of  each  membrane  element  such  that  the  nodes  of  the  membrane 
element  coincide  with  those  of  the  interface  elements.  During  motion, 
Interpenetration  of  the  curved  membrane  surface  by  the  straight  edges 
of  the  Interface  element  will  occur  away  from  the  nodes.  This 
approximation  should  not  affect  the  results  sugnif leant ly  as  the 
interface  elements  are  only  introduced  to  produce  a  resistance  and  do 
not  correspond  to  any  physical  material.  A  single  point  integration 
point  is  used  for  the  3-dimensional  elements;  which  means  the 
generalised  reaction  forces  on  the  nodes  on  the  top  surface  are  equal 
and  opposite  to  those  of  the  nodes  on  the  lower  surface.  2x2  Gaussian 
integration  scheme  is  used  for  the  membrane  elements. 

It  was  not  possible  to  successfully  obtain  results  when  any  of  the 
three  dimensional  models  are  used  to  model  a  multilayer  assembly  for 
impact  velocities  higher  than  160  ms”\  even  in  this  range  the 
numerical  results  obtained  are  not  reliable.  The  numerical 
difficulties  arise  from  the  critical  nature  of  the  compressive 
resistance  of  the  interface  elements  to  a  successful  code 
implementation.  Experimental  data  of  the  relationship  between  the  yarn 
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not  available,  and  therefore  arbitrary  values  were  used.  If  the  assumed 
compressive  modulus  is  too  high,  the  projectile  is  prematurely 
arrested.  If  the  modulus  is  too  low,  on  the  other  hand,  the  downstream 
layers  do  not  provide  any  resistance,  and  some  interface  elements  will 
flip  over,  that  is,  their  depths  become  negative.  However,  the 
difficulty  may  not  totally  stem  from  the  non-availability  of 
experimental  data. 

e 

A  solution  to  this  problem  would  be  to  make  the  marching  step 
iterative  and  treat  the  nodal  forces  due  to  the  compression  of  the 
interface  elements  as  a  sort  of  ' reaction *  forces.  At  the  beginning  of 
the  step,  a  first  estimate  is  made  of  the  nodal  acceleration  at  the 
beginning  of  that  step  and  the  new  nodal  position  vector  and  nodal  force 
vector  are  calculated.  Generally  the  forces  at  the  nodes  would  not 
balance.  The  requirement  that  the  upward  reaction  on  a  layer  above  an 
Interlayer  element  be  equal  to  the  downward  reaction  vector  exerted  by 
the  same  Interlayer  element  only  the  lower  layer  would  not  be  satisfied, 
in  other  words.  The  nodal  acceleration  vector  is  then  adjusted  until 
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the  magnitude  of  the  difference  between  the  vectors  falls  below  a  fixed 
limit  This  modification  of  the  code  was  made  unnecessary  by  the 
success  of  the  master-slave  method,  described  next,  is  modelling  the 
mechanics  of  multilayer  panels. 

5.3.5  MASTER-SLAVE  LAYER  METHOD:  RESULTS 

Initialy,  the  method  was  tested  in  a  model  involving  4-node 
rectangular  elements.  Each  layer  was  divided  identically  and  the 
central  node  in  a  layer  is  given  the  same  displacement  to  stimulate 


impact.  The  method  was  found  to  converge  with  an  increase  in  the  number 
of  elements  into  which  a  layer  is  divided.  A  test  run  was  carried  out 

with  a  panel  of  two  layers.  One  layer  is  a  240gm  Nylon  fabric  while  the 
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second  is  a  351  gm  Nylon  fabric.  The  panel  was  impacted  at  500  ms  . 

In  the  first  run,  the  lighter  fabric  was  put  in  front  (nearer  to  the 

projectile)  and  in  the  second,  the  position  of  the  fabrics  were 

reversed.  In  each  run,  the  projectile  penetrated  the  panel,  24yt  s  after 

impact,  with  a  residual  velocity  of  457  ms”11.  The  coincidence  of  the 

results  of  the  two  min s  confirms  the  stability  of  the  procedure. 

The  Influence  of  the  frictional  resistance  between  the  layers  was 
cheeked  for  values  of  coefficients  of  dynamic  friction  between 
0.0  -  0.03.  Frictional  resistance  against  interlay  slip,  or  the 
surface  treatment  of  the  fabrcs,  did  not  significantly  affect  the 
overall  structural  behaviour  of  the  panel  for  this  range  of  frictional 
coefficient. 

For  a  single  layer  of  fabric,  it  was  found,  that  the  membrane  model 
based  on  4-node  element  always  predicts  that  fabric  failure  occurs 
earlier  than  experimentally  observed.  The  difference  is,  partly, 
because  the  linear  nature  of  the  variations  of  variables  within  the 
4-node  element  prevents  an  adequate  variation  of  strains.  The  assumed 
strain  variation  pattern  which  results  from  the  differentations  of  the 
displacements  cannot  accurately  model  the  rapid  change  of  strains  which 
occurs  in  the  fabric;  especially  in  the  area  close  to  the  impact  point. 
The  model  therefore  predicts  higher  yarn  strains.  The  result  is  that 
the  breaking  strain  is  reached  too  early  during  motion  and  failure 
occurs  at  low  impact  velocities,  for  example,  a  projectile  striking  a 
351  am  Kevlar  fabric  at  500  ms  is  predicted  to  have  penetrated  with 


To  improve  the  prediction  of  fabric  strains,  8-node  elements  were 
subsequently  used  to  obtain  results  comparable  to  the  available 
experimental  results.  Nodal  strains  were  calculated  from  strains  at 
integration  points  by  the  use  of  a  least  square  smoothing  routine  (59). 

Still  using  4  -node  elements,  the  ability  of  the  procedure  to 
model  the  separation  of  layers  during  motions  was  tested.  A  two-layer 
panel  of  a  386  gm"2  Nylon  fabric  and  a  351  gm"2  Kevlar  fabric  was 
modelled.  The  resistance  to  Interlayer  slip  was  put  at  zero.  In  the 
first  run,  the  Kevlar  fabric  was  put  nearer  to  the  impacting  projectile. 
The  Kevlar  fabric  was  penetrated  10^4  s  after  the  start  of  motion  when 
the  projectile  speed  has  been  reduced  to  492  ms-1 .  For  the  remaining 
part  of  the  motion,  the  projectile  was  in  contact  with  the  Nylon  fabric 
which  then  moved  away  from  the  penetrated  fabric.  The  second  layer  was 
penetrated  I4^u  s  later  when  the  breaking  strain  of  nylon  was  reached. 
The  residual  projectile  velocity  was  459  ms-1.  In  the  next  run  the 
nylon  fabric  was  put  nearer  to  the  projectile  so  that  separation  does 
not  occur  throughout  the  motion,  although  the  Kevlar  fabric  failed  when 
the  strain  at  the  centre  of  the  panel  exceeded  0.0323. 

Fig.  5.33  shows  the  deceleration  curve  and  the  change  of  vertical 
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displacement  at  the  impacted  node  for  a  240  gm  Nylon  fabric  impacted 
at  150  ms”1.  The  8-node  rectangular  elements  were  used.  Although  there 
are  no  comparable  experimental  results,  the  numerical  results  seen 
reasonable. 

The  strain  values  were  improved  by  the  use  of  8-node  rectangular 
elements.  The  residual  velocity  with  which  a  1.0  g  projectile  with  an 
Initial  velocity  of  500  ms”1  leaves  a  351  gm”2  Kevlar  fabric  was 


predicted  to  be  472  os"1.  The  energy  absorbed  by  the  fabric,  13.6  J, 
compares  well  with  the  value  of  14.6  J  which  was  obtained 
experimentally  in  UMIST  for  an  impact  velocity  of  512  ms"1.  The 
variation  of  the  energy  absorbed  with  velocity  is  very  weak  at  such  high 
velocities. 

The  posibility  of  modelling  multilayer  panels  as  single  layer 

panel  was  also  investigated.  This  was  achieved  by  multiplying  the  tex 

of  the  yarns  and  the  fabric  density  of  a  single  layer  by  the  number  of 

layers  in  panel.  This  implicity  involves  the  assumption  that  each  layer 

deforms  identically  to  the  other,  which  means  that  inter-layer  slip  is 

insignificant.  This  approximation  is  reasonable  at  higi  velocities. 

Experimental  results  obtained  in  UMIST  (Textile  Technology)  (61)  Showed 
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that  while  9  layers  of  240  gm  Nylon  fabric  are  totally  penetrated  by  a 

projectile  travelling  at  405  ms"1,  10  layers  are  not  penetrated  even 
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when  the  projectile  velocity  is  increased  to  407  ms  .  With  the 

approximate  multilayer  model  Just  explained,  it  was  predicted  that 

| 

whereas  10  layers  of  the  fabric  are  totally  penetrated  at  a  velocity  of 

i 

420  ms"1,  while  the  projectile  is  stopped  when  the  velocity  is  reduced 
to  400  ms"1 . 

Fig  5.34  shows  thu  variation  of  the  energy  absorbed  before 
penetration,  for  increasing  numbers  of  351  gm  Kevlar  layers.  Both 
the  experimental  results  and  the  numerical  predictions  are  shown. 

In  general  it  is  observed  that  this  approximate  method  of 
modelling  multilayer  panel  produces  results  which  are  very  close  to  the 
experimental  results.  However  the  difference  between  the  numerical  and 
experimental  results  increases  as  the  residual  energy  approaches  zero. 


In  thl3  region,  the  size  of  fabric  is  Just  enough  to  absorb  all  the 
initial  energy  of  the  projectile.  The  contact  time  between  the 
projectile  and  the  fabrics,  especially  the  last  layers,  is  then 
relatively  high  and  some  factors,  neglected  in  the  approximation, 
become  significant. 

The  master-slave  layer  method,  with  allowance  for  interlayer  slip, 
that  is  the  system  described  in  section  4.5,  was  successfully 
Implemented  using  8-node  rectangluar  elements.  However,  solution  time 
was  found  to  be  multiples  of  that  required  for  the  same  problem  when 
slipping  was  suppressed.  Host  of  the  extra  time  was  taken  in  the  search 
subroutine  and  in  the  calculation  of  the  normal  components  of  the 
reaction  between  the  layers.  Also,  as  explained  in  section  4.5,  the 
modified  mass  matrix  (eqtn.  4.5.17)  is  then  no  longer  diagonal  and  has 
to  be  inverted.  The  only  condition  in  which  the  original  method  need  be 
used  Instead  of  this  more  efficient  approximation  is  when  the  effect  of 
fabric  slip  is  under  investigation  or  when  panel  is  assembled  from 
fabrics  of  different  materials.  Under  these  conditions,  the  panel  may 
not  be  replaced  by  a  single,  equivalent,  fabric. 


Because  of  the  storage  and  execution  time  advantages  of  the 
approximate  method,  it  is  to  be  preferred  to  the  complete  master-slave 
layer  method,  except  under  special  conditions. 

Experimental  results  (37)  show  that  it  is  more  effective,  in  a  two 
layer  panel,  to  put  the  denser  layer  upstream  instead  of  downstream, 
the  difference  in  effectiveness  cannot  be  predicted  by  this  model.  This 
is  possibly  linked  to  the  fact  that  influence  of  localised  yarn 
compression  on  yarn  rupture  is  neglected. 


the  whole,  this  model  can  predict,  with  reasonable  accuracy, 
the  mechanics  of  the  response  of  multilayer  fabrics  to  impact. 
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CHAPTER  6 

CONCLUSION 

6.0  Introduction 

The  accuracy  and  usefulness  of  the  different  models  in  predicting 
the  V^q,  which  is  defined  as  the  speed  at  which  50f  of  the  projectiles 
fail  to  penetrate  a  test  panel,  is  discussed  in  this  chapter.  The 
computer  times  required  to  obtain  numerical  methods  by  these  methods 
are  compared . 

6.1  The  Simple  Variational  Method 

Although  the  computer  units  used  in  solving  ths  equation  is  small, 
results  are  only  obtainable  for  fabrics  with  high  prestrains.  While  the 
method  is  useful  for  comparing  different  fabric  weaves,  it  is  not 
suitable  for  obtaining  the  VgQ  of  a  panel. 

6.2  Nodal  Impedance  Analysis 

This  analysis  may  be  useful  in  predicting  the  relative  advantages 
of  different  fabric  geometries,  as  far  as  their  response  to  impact  is 
concerned.  However,  the  suggestion  by  two  earlier  investigations  is 
that  crossover  effects  may  not  significantly  magnify  strains  as 
considerable  slipping  occurs  (61),  (62). 

6.3  Method  of  Characteristics 

This  solution  requires  the  maximum  computer  time,  for  a  given 
fabric  and  projectile,  of  all  the  approaches  and  methods  used  in  thi3 
investigation.  Apart  from  this  disadvantage,  a  high  prestrain  which 
increases  with  the  impact  velocity  is  required  for  it3  successful 


implementation.  The  use  of  this  method  should  be  restricted  to 
providing  a  basis  for  comparing  different  fabrics  as  it  is  too  expensive 
for  regular  usage. 

6. 4  Finite  Element  Solutions 

The  various  models  are  treated  separately. 

(a)  Rod  elements  -  (2  node  and  2  node) 

These  produced  well-conditioned  results,  with  no  oscillations  for 
some  fabrics.  But  oscillations  occur,  especially  for  Kevlar  fabrics. ' 
This  model  provides  a  means  of  investigating  the  effect  of  rigidity  of 
yarn  connections  at  crossovers.  However  methods  of  supressing  the 
strain  oscillations  have  to  be  introduced  before  the  model  could  be  used 
in  production.  The  model  based  on  2-node  rod  elements  makes  the  most 
efficient  use  of  computer  time,  compared  with  all  the  other  models. 
They  also  provide  an  approximate  approach  to  assessing  the  effect  of 
crimp. 

(b)  The  3  node  (triangular)  membrane  Element 

This  is  an  accurate  model  for  investigating  the  mechanics  of  a 
single  fabric  layer.  While  the  results  are  well-behaved  for  most  input 
data,  oscillations  have  to  be  supressed  for  some  input  data.  The  model 
is  quite  efficient  in  computer  time  usage.  It  is  useful  for  different 
weaves  of  fabric  and  also  for  knitted  fabric. 

(c)  The  3-Dimensional  Model 

Further  work  is  required  to  render  this  model  usable  for  routine 
design.  The  solution  might  be  to  replace  the  two-dimensional  membrane 
elements  by  20-node  parallelopiped  and  remove  the  interface  elements. 

(d)  The  Master-slave  layer  method 

This  model  involving  4-node  rectangular  elements  produces  poor 
predictions  of  values.  However  it  is  useful  as  an  inexpensive  way 
of  studvinff  hh»  rtf*  •  n*» 
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Better  results  are  obtainable  by  using  8-node  rectangular 
elements.  However  the  running  times  are  high,  compared  with  either 
that  involving  3-node  triangular  elements  or  4-node  triangular 
elements.  This  code  is  useful  for  routine  calculations. 


6.5  Factors  affecting  overall  behaviour  of  fabric 

Ally  a  description  rather  than  a  quantitative  comparison  of  the 
significance  of  the  various  yarn  and  fabric  properties  would  be  given. 

In-Plane  Shear 

Shear  stiffeners  seem  to  have  little  significance  on  fabric 
performance  until  its  value  is  high  enough  to  reduce  the  ability  of  the 
yarns  to  absorb  the  projectile  kinetic  energy  as  longitudinal  strain 
energy.  This  will  happen  when  a  very  stiff  coating  is  applied  to  a 
fabric.  Otherwise  the  change  in  fabric  shear  stiffness  obtained  by 
washing  may  not  be  significant. 

Crimp  interchange 

The  demonetisation  has  not  been  fine  enough  for  definite 
conclusions  on  the  significance  of  crimp  interchange.  Its  effect  on 
fabric  performance  may  vary  with  the  impact  velocity. 

Friction 

Increasing  the  friction  of  yarn  crossovers  will  increase  the 
percentage  of  reflections  at  nodes,  and  will  therefore  reduce  fabric 
effectiveness.  Friction  between  layers  does  not  seem  to  be  very 
significant  close  to  the  VgQ  since  penetration  occurs  before 
appreciable  movement  of  the  panel  can  take  place.  The  influence  of 
friction  —  generated  heat  on  modulus  and  breaking  strain  was  not 
evaluated . 

Size  of  yarn  and  pitch 

The  tex  of  a  yarn  seems  to  affect  the  value  of  its  breaking  strain 
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whUe  its  crimp  varies  with  pitch,  for  a  given  tex.  Apart  from  these 
indirect  effects,  neither  the  text  of  a  yarn,  nor  the  pitch  in  any 
direction,  seem  to  significantly  affect  fabric  performance. 

Bending  stiffness  of  yarn 

Although  some  earlier  investigations  conclude  that  yarn  bending 
cannot  be  neglected  if  the  mechanics  of  fabric  behaviour  at  low  strains 
is  to  be  accurately  modelled,  the  results  obtained  show  that  the 
Influence  of  bending  may  be  of  second  order  to  the  overall  fabric 


effectiveness 
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6.6  Future  Work 

(1)  The  condition  for  yarn  failure  which  was  adopted  here  is  the 
simplest  possible.  Since  the  order  of  arrangements  of  two  fabrics  has 
been  found  to  affect  f-.eir  efficiency  in  absorbing  impact  energy,  it  is 
likely  that  yam  compression  affects  the  longitudinal  strain  at  which  a 
yarn  ruptures.  The  relationship  should  be  studied  and  incorporated 
into  the  failure  condition  U3ed  in  any  numerical  method  for  predicting 
V5Q  values. 

(2)  The  attempts  to  model  a  fabric  with  yarns  which  can  slip  during 
motion  has  not  been  very  successful  due  to  numerical  oscillations  of 
yarn  strain.  The  development  of  an  alternative  method  should  be 
investigated. 
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APPENDIX  1 


Wavefronts  in  Triaxial  Weaves 


Prediction  of  Wave-Front 

Leech  and  Mansell  (31)  have  shown  that  the  wave  front  in  a 
transversely  impacted  orthogonally  woven  fabrics  is  rhombiocal. 
Below,  it  is  shown  that  for  a  plain  triaxial  weave  the 
propagation  front  is  hexagonal.  The  repeated  cell  of  this  fabric  is 
shown  in  Fig.  A1.1(a). 

The  equation  of  the  sides  of  the  hexagonal  are 


Al.la 


11. 1b 


11.1c 


where  S^,  s2  and  are  the  coordinate  in  the  directions  of  the 
yarns  (Fig.  .11.1b)  and  C2  and  are  the  transverse  wave  velocities 


in  the  respective  directions. 


1  similarity  type  of  solution  W(t1s1,s2,s3)  a  Wo(t)W  (  )  is 

S.  S- 

introduced  where  a  g— ^  ^  for  the  ragion  where  0, 

Sg  ^  0,  0  where  Wo  is  the  displacement  at  the  point  of  impact. 


See  Fig  A1.2. 


The  coordinate  is  perpendicular  to  the  wavefront  and 
the  impact  and  ^  a  1  at  the  wavefront  Fig.  A1.3. 

The  excited  area  A  s  t2C,C2  ^  2  and  the  incremental 

dA  a  t2C^Cz  ^  • 

T,  the  kinetic  energy  of*  the  material  in  area  dA  behind 
front,  at  any  time  t,  is  given  by, 

[?“<  fi  +  h  +  h*  (pO  dA* 

g  ^  are  the  mass  fractions  of  the  yam  families 
3).  The  total  kinetic  energy  for  the  whole  deformed  is 

T  a  l/|  (0^3  ♦  C2Cy  ♦  C2C3)t2  ^Wo2I1  ♦  2  ^  ♦  Wof 

t 

where  Ij  *  J  Htj  drt  ,  and  I2  a  ^W1  d  7 

Hote  (•)  denotes  differentiation  with  respect  to  ^ 
differentiations  with  respect  to  time. 

The  expression  for  the  strain  at  one  point  in  the  yam  s 

l  *  ^  (i  ♦  £s)  / 1  -  1  J 


T  a 

Pi’  h’ 


*\  s  o  at 

area, 

the  wave 

A1.2 

(1,  2  and 

|  ,  A1.3 

and  ( * ) 

,  is 


where  £  is  the  prestrair.  and 
9 

that  point. 


3  W 
*  3 


is  the  slope  along  the  yam  at 


Fop  a  linearly  elastic  material,  the  resulting  strain  energy  = 

£5  J  v.  [<’♦  /  '♦ 

A 

Sunned  up  for  all  the  familiarty  strands  in  the  area.  E  is  the  Young's 

modulus  of  the  sth  strand  and  h  is  its  average  thickness. 

s 


Assuming  highly  stressed  triaxial  fabrics  with  small  identation, 
the  system  is  now  linear  '  « 1 )  and  the  strain  energy 
becomes , 


U3f  [p^t*4 of)  KWi> 


A1.5a 


or 


—  i.  ■ 

1  *  2 

The  initial  kinetic  energy  of  the  projectile  is  j  MpWo  .  When  this 
is  added  to  the  total  energy  of  the  fabric  and  Hamilton's  principle  is 
applied,  the  result  is  the  equation  of  {notion: 

Jmp  ♦  ^Bm(C1C3+C2C1+C3C2)t2I1  Wo  ♦  2^3! mtC^+CjCj+C^Jt^Wo 
j3m(C1C2+C1C3+C2C3)(I1  -  I2)^J  Wo  =  0  , 


where  m  =  ^7  ^  C\,  ^  q-j 

ca  c}  ' 


A1.6 


The  equations  can  be  nondimensionalised  by  introducing  the  factor 
£  a  o(  t  and  \  t  to  get 


(1  ♦  *£Z)  Wo"  +  2  1  Wo* 


\I3  «■  (I,  -I2) 


Wo  =  0  A1.7 


2  /3m(C1C2+C2C3+ClC3) 

where  04  *  - T2~TTp -  and 


A  • 

^  (£*£*.  *  ^1^3  4  C»C*) 

Equation  3.7  is  very  similar  to  the  expression  (I  ♦  •£  )Wo"  +  21 
•  Wo'  +  S  Wo  s  0  obtained  for  small  indentation  of  an  orthogonally  woven 
fabric  (31)  where 


*3  *  h  -  T2 


The  only  parameter  which  affects  the  response  of  the  orthogonal 

fabric  for  a  nondimens ional  velocity  of  impact,  is  S  which  is 

1 

dependent  only  on  the  grouping  of  integrals  while  for  the  trlaxial 


fabric  X  also  depends  Ion  the  ratio  of  the  yarn  stresses. 


Large  Deformations 

The  equations  of  motion  for  large  deformations  of  the  triaxial 
fabric  is  obtained  by  the  application  of  Hamilton's  principle  with  the 
expressions  A1.3  and  A1.4  incorporated  and  is  as  follows 


i 


'  L  vNC n'o/c.z)"  J 


+jc,Klt£jj U - 1  7  Wirfjf.  __  _j _ 7  "'8 

1  *  'fcOWidv  (  r7)+(w./c,tfjl 


where  k 


1 


C1C3  +  C1C2 
0^2  +  C^Cg  +  C2C33 


cic2  *  C2C3 

^  C2C3  +  c3c1 


C2C3  ♦  ,  and  Pi*  p2*  ?3  are  the 


C^Cg  +  C2C3  + 


mass  fraction  of  the  3,,  s0  and  3,  3trand3  respectively 


APPENDIX  2 


Theory  of  the  Method  of  Characteristics 


Theory 

Consider  a  single  filament,  Fig.  A2.1,  with  displacement 
functions,  u,  v,  w  relative  to  its  initial  constraint  position.  The 
displacements  are  functions, of  x,  the  running  coordinate  along  the 
undisplaced  length  of  the  filament,  and  t  the  time. 

To  allow  for  large  deformations,  the  strain,  £  ,  is  defined  as 


£  +  w,;  j1-  1. 


(Note:  h,^  denotes  the  differentiation  of  h  with  respect  to  x). 

The  kinetic  co-energy  and  the  strain  energy  can  be  written  as: 


i  Lf  p*(0’t2  ♦  v»t2  +  V’t 


)  dX  , 


where  A  is  the  cross-sectional  area  of  the  yarn  and  E  its  Young's 
modulus,  and  ^  its  density. 

Applying  Hamilton's  principles,  the  differential  equations  of  the 
motion  of  the  string  are: 


”’tt  *  f  <1  )2  ££0  *  £)0Vx 


+  0  +  0,-  )  £ 


A2.3(a) 


v.tt  -  - - 2 

(1+  £.  T 

l" 

♦  t  )V,^ 

+  V|  A 

i.x) 

=  0 

A2.3(b) 

W,tt  -  - - 

(1+i)2 

L" 

+  l  )W,** 

+  w,x 

£.xj 

s  0 

A2.3(c) 

These  equations  are  conveniently  reduced  to  a  system  of  six 
quasilinear  hyperbolic  equations  by  the  introduction  of  the  following 
terms: 


p  »  n.t  .  «  *  V It  .  r  =  W,t  ) 

9  f  •«.  X  • 

Further  at  a  point  on  the  filament,  the  vectors  Z,  W,  and  £  are 
defined  as 


t 

Is 


These  6  equations  are  now  written  in  matrix  form  as 

Z,.  ♦  AZ.  >  0,  A2.5 

where  A  is  a  6  x  6  matrix. 

Equations  A2.3  are  the  equations  of  propagations  of  signals  or 
disturbances  along  the  filament. 

Equation  A2.3(a)  models  longitudinal  wave  motion  while  A2.3(b)  and 
A2.3(c)  are  the  equations  of  propagation  of  waves  in  the  two  transverse 


directions 


Characteristic  theory 

is  now  applied  to  equation  A2.5 

characteristic  equations  of  the  single  yarn.  These  are 

along 

dx  - 

dt  =  +  e 

I* 

D1  (W)  dZ 

=  o? 

dx  _ 

/«/. 

along 

dt  a  “  Ce 

M  P 

I" 

02  (W)  dZ 

a  0  , 

along 

dx  „ 

dt  s  *  °T 

=  /e£  /p(l  ) 

II* 

U3  (W)  dZ 

a  0  , 

<£>  dZ 

s  0 

along 

a  c 

dt  8  CT 

t/pd  ♦  £  ) 

II" 

05  (W).dZ 

3  0  , 

A2.6(a) 


A2.6(b) 


A2.6(c) 


A2.6(d) 


A2.6(e) 


°6  *  C.  A2.6(f) 

In  each  of  the  equations,  dZ  ia  the  change  in  the  vector  Z  along 
the  characteristic  curve  see  Fig  A2..2,  while 

tJj,  i  *  ^ . 6  are  the  6  left  eigen  vectors  of  the  matrix  A 

(equationAl..  12) . 

It  is  to  be  noticed  that  the  velocities  of  propagation  of  the 
transverse  waves  C^,  C,y  are  dependent  on  the  strain  and  so  vary  with  time. 

Also  C@  ^  C_  for  all  values  of  of  interest. 


The  eigenvectors  of  the  matrix  A,  equation  A2.5  are: 

°1  =  °2  *  ♦  0)»  0*  f  *  Ae(1  ♦  0),  Xe0tXe^J 

°3  =  °5  =  £  0»  (1  +  0)>  °»  At2»  AfcO  ♦  ®)»  Oj  T 

s  Og  s  •  *  (1  *  0) t  0*  Atl|/  »  ~  At(1  ♦  0)  f  1 

# 

The  equations  A2.6  are  now  rewritten  as 

(1  ♦  8)  A  Pj  *1  Tf'ii  ri  ♦  c  <1  ♦  8)  *  c  0A0 

c.l'iM'  a  0  (I  -  extensional) 

0  •  A.  P±  -  (1  +  0)  Aq.+Ct0  Act  -  Ct(1  ♦  Q)Atyi  »  0 
(II“  characteristies-inplane  transverse). 

•  Avi  -0+0)  +  ct\jf  Aa1  -  cT(i  ♦  0)  A  Y’i  *  ( 

(H“  characteristics-out-of-plane  transverse) 


(d)  Special  triaxial 
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Percentage  of  in  plane  signal  reflected  for  different  yam  arrangements 
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m  *  m+1 
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^  fT 
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Fig.  3.5 
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Fig.  4.3  Direction  of  frictioncl  force  (Kfh  Estimate) 
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Fig. 4.1 1  DiagrGmatic  representation  of  search  procedure. 
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Fig.  5.2.  Linear  model :  Deceleration  graph:  Different  area  densities. 
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7.  Variation  of  strain  with  impact  velocity-  2 materials 
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Transverse  deformation  against  time 


rig.  5.9  Results  obtained  by  method  of  characteristics 


Displacement  (m)xl0~2  Velocity  (ms^JxlO 


si 

a 

I 


Data :  Trilinear  nylon.  «p  *  0-03 ,  Vp= 235  ms"1 
3  layers.  Area  density  =  3x240 gm"2 
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Fig.  5.13  Typical  characteristic  solution- 
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Data:  Trilinear  nylon. €P=003,  Vp  =  245ms 
3  layers.  Area  density  =  3x240 gm“2 

NYLON 


25 

20 

15 


5 


0  5  10  15  20 

Projectile  arrest  in  193/i.s 


0  5  10  15  20 


Profile  displayed  after: 
43^xs,  8l^.s,  I20^,s,  I6l^xs,  I93yxs 

Fig.  5.14  Typical  characteristic  solution 


Material:  Nylon-Linear  elastic  (240gm  ,103  tex)  single 
layer  biaxial  weave 
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Fig  5.20 


Response  of  two  weights  of  fabric  to  impact  at 
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Fig.  5.27  Typical  results  obtained  using  2-riode  rod  elements 


Fig. 5.29.  Strain  history  at  impact  point  for  two  crossover 


o  Simplified  master-slave  layer  model 
(8-node  model) 

□  Experimental  results 
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Fig.A1.1b  Hexagonal  wave  front  in  triaxially  woven  fabric 


